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Film-Lubrication between Spherical Surfaces: with an Application to the Theory of 
the Four-Ball Lubricant Testing Instrument 
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(Received July 19, 1945) 


The paper is concerned mainly with the flow of viscous 
liquid between the outer surfaces of two equal spheres 
placed in an infinite body of the liquid so that the “mini- 
mum separation”’ of the surfaces (the distance apart at the 
points of closest approach) is small compared with the 
radius; one sphere is held at rest, the other rotates at con- 
stant speed about an axis which does not pass through the 
center of the fixed sphere. The assumption is made that the 
surfaces are separated by a continuous film of liquid to 
which the classical theory of hydrodynamical lubrication 
can be applied—that is, ‘‘boundary lubrication” and actual 
contact between the surfaces are excluded; the partial 
differential equation (the Reynolds equation) for the 
pressure distribution in this film is set up, and it is shown 
that with a suitable choice of boundary conditions a solu- 
tion can be found in very simple form. From this solution 
are deduced the total force exerted by the liquid on the 
moving sphere and the torque, about the axis of rotation, 
of the forces acting on the fixed sphere. The results are 
applied to the particular case of the four-ball lubricant 
testing instrument, in which a sphere is made to rotate 
about a vertical axis, under axial load, in the central space 
formed by a set of three equal spheres held stationary so as 
to touch one another with their centers in a horizontal 
plane; the whole set of spheres is immersed in the lubricant 
to be tested. The relations finally obtained express the load 


on the moving sphere and the torque on the set of fixed 
spheres in terms of the radius of the spheres, the velocity of 
rotation, and the viscosity of the liquid; they involve as a 
parameter the minimum separation of the moving sphere 
from the fixed spheres, and elimination of this gives the 
relation between the measurable quantities, load, and 
torque, which could provide a basis for the use of the 
instrument as a viscometer. It is found that the torque is 
a very slowly varying function of the load, being approxi- 
mately proportional to the logarithm of the load. Further, 
if it is assumed that the hydrodynamical theory of lubrica- 
tion is valid only for films whose thickness is greater than 
some agreed minimum, then the expression for the load on 
the moving sphere, into which the minimum separation 
enters, provides a criterion for the range of conditions (of 
which the load is much the most important) within which 
the theory does in fact apply to the instrument. Numerical 
examples are considered, and it is shown that, with 
allowance for wide departures from the standard size and 
speed of operation of the instrument, the maximum load 
which can be supported by the moving sphere in the 
hydrodynamical regime cannot be expected to exceed a few 
hundred grams weight; for loads of greater order the 
standard continuous-film theory cannot hold, and the 
forces on the spheres are no longer determined by the 
viscosity of the liquid. 





INTRODUCTION 


- the last few years a great deal of attention 
has been given to the properties of lubricants 
under conditions of ‘‘extreme pressure,’’ that is, 


under pressures so great (of the order of 10° 
kg wt./cm?) that it is not to be expected that the 
lubricated surfaces are separated by the con- 
tinuous film which characterizes hydrodynamical 
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Fic. 1. (A) Elevation and (B) plan of the 4-ball instru- 
ment. Ci, C2, C3 are the centers of the stationary spheres 
and C, that of the rotating sphere. 


or “thick film” lubrication. The literature of the 
subject has many references to the Boerlage 
four-ball testing instrument.! 

The instrument seems to have become a stand- 
ard piece of apparatus in this field, but at the 
same time it seems that the interpretation of the 
measurements made with it is not always certain. 
There is an obvious question which must present 
itself in any general discussion of the theory of 
the apparatus, and which should be amenable to 
theoretical investigation; this is, can the instru- 
ment behave in any circumstances according to 
the laws of hydrodynamical lubrication, and if it 
can, what is the range of operating conditions 
over which this is possible? There seems to be no 





'G. D. Boerlage, Engineering 136, 46 (1933). O. Beeck, 
J. App. Phys. 12, 512 (1941). D. Clayton, Inst. Mech. Eng., 
Discussion on Lubrication 2, 282 (1937). R. Schnurmann, 
Engineering 149, 567 (1940). 
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published investigation of this point; the present 
paper attempts a treatment based only on 
classical hydrodynamical theory. 

We may recall here that the four-ball instru- 
ment consists essentially (see Fig. 1) of a set of 
three equal spheres held rigidly in a fixed cup so 
as to touch one another with their centers in a 
horizontal plane, with a fourth sphere rotated 
under load about a vertical axis in the central 
space of the fixed spheres; the whole is immersed 
in the lubricant to be tested. Details of some 
methods of use with great loads (which are not 
really relevant to our present purpose) are given 
in other papers.' We base this investigation on 
two fundamental assumptions; first that the 
spheres are not deformed by the load, and 
second, that the flow of liquid around and be- 
tween the spheres obeys the classical equations 
of viscous flow, in particular that a continuous 
film in the hydrodynamical sense (as distinct 
from the ‘boundary lubrication” sense) is formed 
between the moving sphere and each fixed sphere. 

We attack the problem indirectly by first 
developing the relation between the thickness of 
this film at its thinnest part (the “minimum film 
thickness,’’ say) and the load applied to the axis. 
Suppose we find a relation W=W(ho) between 
the load W and this minimum thickness h»— 
involving of course all the parameters of the 
problem—the radius of the spheres, the angular 
velocity of the rotating sphere, and the viscosity 
of the liquid. Suppose also that we have reason 
to believe that a film of thickness less than some 
value d does not behave as a lubricating film in 
the hydrodynamical sense; then we may deduce 
that the assumption of hydrodynamical lubrica- 
tion ceases to hold for loads exceeding the value 
W(d) corresponding to the thickness 4y=d. Con- 
versely we may deduce that for loads less than 
this value there will be in fact lubrication of the 
hydrodynamical type. We have thus a criterion 
for the range of application of the classical theory. 

A complete and rigorous investigation of the 
problem would involve the solving of the full 
equations for the flow around the spheres, and 
with a geometry so complicated no success could 
be expected from such an undertaking: indeed 
it was shown by Lamb? that even the apparently 


*H. Lamb, Hydrodynamics (Cambridge University 
Press, London, 1932), p. 589. 
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simple problem of a sphere rotating in an infinite 
body of liquid is intractable if the velocity ex- 
ceeds a very small value. To bring the problem 
within the scope of reasonably simple analysis 
we make the following further assumptions: 

A. The flow in the immediate vicinity of the 
moving sphere and any one of the fixed spheres 
is unaffected by the presence of the other 
spheres. 

B. The film formed between the moving sphere 
and any fixed sphere is sufficiently thin for the 
standard assumptions of classical film-theory to 
hold ; that is, 


(a) the inertia terms in the equations of motion can be 
neglected, 

(b) the pressure can be regarded as constant across the 
film, and 

(c) rates of change of velocity measured across the film 
are much greater than those measured in any other 
direction. 


There is also the over-riding assumption, implied 
in the use of the hydrodynamical equations, that 
the only relevant physical property of the lubri- 
cant is its viscosity; and we shall take this as 
constant throughout the film. 

With these assumptions the problem is reduced 
to the solving of the Reynolds equation for the 
flow between one pair of spherical surfaces, that 
is, the partial differential equation for the 
pressure distribution in the film. The general 
outline of the investigation is then as follows: 
First ($§1, 2) the Reynolds equation is set up 
and the boundary conditions are obtained. It is 
then shown that with a suitable choice of 
boundary conditions there is a very simple 
solution which is computed as a numerical 
function. From this we deduce (§4) the total 
load supported by the film; and a further 
analysis (§5) leads to the expression for the 
viscous-frictional torque exerted on the station- 
ary sphere. Finally we apply these results to the 
particular case of the four-ball instrument and 
derive the criterion for the range of working 
within which we may suppose hydrodynamical 
lubrication to occur. 

The final conclusion is that the usual condi- 
tions of heavy loads are far outside this range, 
as is to be expected, so that in these circum- 
stances the relevant property of the lubricant is 
not its viscosity. The results have other points 
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Fic. 2. Geometry of the 2-sphere system considered. 


of interest, however ; for example, the instrument 
could be used as a viscometer, for which purpose 
it would have the advantage of permitting wide 
variations in the working conditions, including 
the rate of shear; the relations obtained here 
provide a basis for the calibration of the instru- 
ment for that purpose.’ It was partly to this end 
that the determination of the torque was made. 
Also there is the point of mathematical interest 
that the investigation leads to an exact solution 
of the Reynolds equation, which is very seldom 
possible; it has in fact been stated‘ that the 
equation cannot be solved in finite terms if the 
boundaries of the flow have curvature in two 
directions. 


1. THE PRESSURE EQUATION 
(REYNOLDS EQUATION) 


In Fig. 2, Ci, Cz are the centers of two spheres 
of radius a; the first sphere rotates with angular 
velocity w about the line MC, (which in the four- 
ball machine is vertical) and the second sphere 
is at rest; the inclination of MC; to CC. is 
taken as an arbitrary angle a whose value in the 
final application will be deduced from the 
geometry of the instrument. 

* It would be necessary to make a series of calibrations 
with liquids of widely differing viscosities to allow for 
possible changes in the dynamical pattern of the flow. 


4See, for example, R. O. Boswell, Film Lubrication 
(Longmans Green and Company, New York, 1928). 
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The coordinate system Oxyz is chosen so that 
(a) O is the mid-point of CiC2, (b) Oz coincides 
with the line of centers C,C2, and (c) the axis 
MC, lies in the plane zx (in the figure, M lies 
on Ox). 

Let p be the pressure and v the fluid velocity 
at any point; let the thickness of the film at 
any section be 2h, and let the minimum thick- 
ness be 2/9. Thus in the figure NP; = NP.=h, 
OA, =OA2=hoand P,P: have coordinates (x, y, 2) 
and (x, y, —h), respectively. 

By assumption A, p is a function of x and y 
only, and the whole treatment of the problem 
depends on the partial differential equation in x 
and y satisfied by p. The derivation of this 
equation is a standard piece of analysis, but as 
we shall need in a later section some of the 
intermediate equations we may conveniently in- 
dicate here the main steps in the process. 

If the fluid velocity v has components 4, v, w, 
then the assumption B permits us to reduce the 
general equations of motion to 


dp du dp daw dp 
=~, —=p—, —=0. (1) 
ax dz* ay 02” «Oz 


Let the velocity of the general point P; on the 
surface of the moving sphere have components 
U, V, W. Every point on the second sphere has 


zero velocity, so we can write the boundary 
conditions for these equations 


at z=h u=U v=V w=W;: 
at s=-—h u=0 v=0 w=0. 


Hence by integration with respect to z 


sth 2°—h® ap 








“= U+ —, 
2h 2u ox 

(2) 
sth 2°—h* dp 
v= V+ —, 
2h 2u «Oy 


Now the flow must satisfy the continuity equa- 
tidn div v=0 
Ou dv dw 


i.e., —+—+—=0. (3) 
Ox Oy a2 


If we substitute for u and v the values given by 


(2) and then integrate the resulting equation 
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for w between the limits —h and h so that 


‘(du dv * dw 
f (~+—)as- -f —dz= —W, 
—p\dx dy _p OZ 
we shall have eliminated the velocities u, v, w and 
derived an equation for p alone. We obtain 


0 
(ot) (0% 
oy oy 

OV 


0 
=}uW+ tur(—+"_), (4) 


Ox oy 


which is the Reynolds equation, in rather more 
general form than is usually quoted. 

So far the argument has been quite general 
and Eq. (4) holds for any system of the type we 
are considering; we must now obtain the ana- 
lytical expressions for h, U, V, W which relate to 
our particular problem. 

The expression for the film-thickness / is 
easily obtained, for if (x, y, 4) is any point on 
the upper (moving) sphere, the distance from 
the center of the sphere is a, and so 


x+y"?+(a+ho—h)? =a’, (5) 
which defines h. 
It is convenient to use the variable 


a+h—h=Z. (Sa) 
Then 
e+y+Z*=a’, (Sb) 


and the derivatives of 4, which we shall need, are 


oh OZ <x 

a ees (5c) 
Ox ox Z 

oh 0Z 

pO (5d) 
oy dy Z 


The velocities U, V, W are most readily de- 
rived by a vector method. Let i, j, k be the usual 
unit vectors along the coordinate axes; then 
(Fig. 2) the vector angular velocity of the 
sphere is 


® = (sin ai+cos ak), 
and if r is the vector C,P; the velocity of P; is 


wAr. Now 
r=xi+yj—Zk, 
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sO 
UVi+ Vj+ Wk=.(sin ai+cos ak) A (xi+yj—Zk), 


and equating the coefficients of i, j, k on the two 
sides of this equation gives the required ex- 
pressions : 

U= —wy cos a, 

V=(x cos a+Z sin a), (6) 

W=wy sin a, 


from which we have 
aV wy 


aU 
—=0, —=-~— sina. 
Ox oy 


Substitution of these values in (4) gives as the 
fundamental equation for the problem 
ap ep 3/7 ap ap 
ermine Sp 
dx® dy? 2hX\ dx 





“ay 


3yuw sin a h 
-——(1 -—)y. (7) 
2h* z 


We can at once simplify the form of this equation 
by changing from Cartesian coordinates x, y to 


polars r, 6 defined by 
x=rcos6, y=rsin @. 


The transformed equation is 


ap 1 3r\dap 10°) 
—+(-+ \o+5 
or? r Zh/or_ r* de 


3uw sin a h 
(1 -—)r sin @, (8) 
2h* 








4 


where 
h=at+h—-Z, Z2=a-r’, 


and a most important point is that # and Z are 
functions of r only. 

In this as in most problems it is of great 
advantage to work in terms of dimensionless 
variables; a little algebraic manipulation shows 
that if we define the dimensionless quantities 
p, II, so by the relations 


p=r/(2aho)', T=p/3uw sin a(2a/ho)', so=ho/a, 
the equation can be written 


1 d°7Il 


é7I1 dll 
L —=Msin@, 
g2 
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where L, M are functions of p only, given by 





Z 1 3p 
=—+————_(1 + $o— 25 op’), 9 
‘ (12500)! 0 op”) (9a) 
1 1+so+(1—2sop?)! 3 
M=p—- | 
2 1+5o0+p? 
[2 —(1+50)(1—2sop?)*]. (9b) 


If finally we regard so=ho/a as a small quantity. 
and in (9a, b) neglect so in comparison with unity, 
we obtain the simpler forms 





1 6p 
L=-+ : (9c) 
p 1+ 9? 
p 
M =———_.. (9d) 
(1+p*)? 


In this last case we can write the equation in 
the form 





p 1+ */ dp 
1 onl p 


p> 06° (1+ ’)? 


re 0 (- 6p \OIl 
Op? 





sin @, (10) 


which is remarkable in that it is completely free 
from all the parameters of the problem: all these 
have been absorbed into the dimensionless vari- 
ables II, p. This single equation therefore covers 
all possible values of a, u, a, w; but in dealing 
with this form we must remember that the 
original Eq. (9) has a singularity at r=a (where 
Z=0) so the range of p must be restricted to 
correspond to a range of r stopping well short of 
this point. This is, of course, quite in accordance 
with physical ideas; for at all points on the circle 
r=a the tangent planes to the two spheres are 
parallel to the z axis (to put it less formally, the 
surfaces are diverging infinitely rapidly) and well 
before this the film will have become thicker than 
the basic assumptions will permit. 


2. BOUNDARY CONDITIONS 


Before we can make any use of the pressure 
equation we must ascribe to p a set of boundary 
conditions, and it is here that the real difficulty 
of the problem is met. We are dealing with an 
elliptic equation for p in the two independent 
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Fic. 3. Geometry of the complete journal bearing. 
S, O are the centers of the journal and the bearing, 
respectively. 


variables x, y (or the equivalent variables r, @) 
and general theory tells us that to obtain a unique 
solution we need to know the values of p at all 
points of a closed curve in the xy plane. The 
possibility of finding a solution in reasonably 
simple form is intimately bound up with the 
choice of this curve and of the set of values of p 
along the curve. 

The difficulty arises because the film has no 
definite boundaries defined a priori by the 
geometry of the apparatus, but merges gradually 
into the main body of the fluid where the as- 
sumptions A and B are not valid. This is not to 
say that the problem is not defined physically : 
strictly the boundaries are the free surface of the 
liquid, the inner surface of the cup, and the four 
spherical surfaces; and the problem is defined by 
the obvious conditions of pressure and velocity at 
these surfaces. But in view of the small chance of 
finding a strict solution we must look for an 
alternative set of conditions for our simplified 
system, relying on general physical principles in 


making our choice. 


‘An obvious condition is that in the main body 
of the liquid, and in fact everywhere except in the 
immediate vicinity of the point of closest ap- 
proach of the two spheres, the pressure must be 
practically atmospheric; if, therefore, we choose 
the atmospheric value as the zero of pressure the 
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solution for p must reduce to very small (nu- 
merical) values everywhere except in this region, 
that is, for all except small values of r. This can be 
idealized to the statement 


p=0 for r>ro, 


where ro is some value considerably less than a, 
and suggests that for the boundary condition we 
might take 


p=0 for r=ro, O<0<2z, 


i.e., the vanishing of p on a circle of radius ro. 
This however is not a suitable condition ; for the 
equation for p has what we may call “skew 
symmetry” about the origin: if in (7) we replace 
x,y by —x, —y we merely change the sign of the 
right side. We get the same effect by replacing @ 
by 6+ in (8) or (9). If therefore the set of 
boundary values has the same property, the 
whole solution will itself have this property ; that 
is, in polars, if in the region 0<@<~7 the solution 
is p=p(r, 0) then in the region r<@<27 the 
solution will be p=p(r, 6+7) = —p(r, 0). Since 
the suggested set of values has in fact this 
property, the use of this condition will lead to a 
solution in which positive and negative values of 
equal magnitude appear in pairs; the immediate 
conclusion, that the resultant load supported by 
the film is zero, is unacceptable, and not less so is 
the implication that in half the film the liquid is 
in a state of negative pressure. 

At this stage we can profit from a glance at the 
theory of the complete journal bearing, the 
geometry of which is shown in Fig. 3. The 
classical treatment of Sommerfeld’ deals with a 
bearing of infinite length (so that no flow in the 
axial direction need be considered) and effectively 
solves the Reynolds’ equation corresponding to 
the flow in the annular space between the two 
cylindrical surfaces, with the condition that the 
pressure and its derivatives are continuous and 
periodic functions of the angular coordinate ¢; 
the striking feature of the solution is that it gives 
positive pressures in the region APB(O<¢<r) 
and equal negative pressures in the region 
AQB(x << 2r). This prediction of negative pres- 
sures has led to much criticism of Sommerfeld’s 
work ; and it seems to be generally agreed that his 
premises depart seriously from the truth in that 


5 A. Sommerfeld, Zeits. f. math. Physik 50, 97-155 (1904). 
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they make no provision for the replacement of the 
lubricant which is forced out of the ends of any 
real bearing, however long; in practice a bearing 
is connected to a source of the lubricant so that 
the leakage is continuously replaced. A rigorous 
treatment incorporating this feature would re- 
quire a strict solution of the Reynolds’ equation 
taking into account both radial and axial flow, 
and this seems to present insuperable mathe- 
matical difficulties; but an approximate solution 
developed by Morgan and Muskat® suggests that 
the introduction of a source does in fact suppress 
the negative pressures. A less sound but much 
simpler treatment has been given by Giimbel ;’ 
this writer, having assumed the connection to a 
source at atmospheric pressure, say at the point 
of maximum film-thickness ¢=0, then neglects 
the axial flow and suggests three possible em- 
pirical conditions to determine p: 

1. The pressure falls to zero at the section of 
minimum thickness and remains zero until the 
position of the source is reached (traveling round 
the bearing in the direction of rotation of the 
axle), i.e., p=0 when ¢=7, and for r<@< 27. 

1]. The pressure has its maximum value at the 
section of minimum thickness and is zero between 
that point and the position of the source, i.e., 
dp/do=0 when ¢=7, and p=0 for r<@<2r. 

Ill. p and dp/d@ vanish simultaneously (the 
position at which this occurs being determined by 
the two conditions) and ~=0 between this point 
and the position of the source. 

What Giimbel assumes is in effect that the 
lubricating film breaks from the moving surface 
either at or just beyond the narrowest section; 
the first two suggested conditions imply a dis- 
continuity in either p or dp/d@; the third is a 
little less drastic but no less arbitrary in essence. 
Nevertheless Giimbel, using these as working 
hypotheses, finds that the numerical solutions 
corresponding to the different possibilities do not 
differ seriously amongst themselves and are in 
reasonable agreement with experimental ob- 
servations. He therefore concludes that his treat- 
ment, if not well founded logically, does at least 
give useful numerical results. A further applica- 
tion of what is effectively Giimbel’s third condi- 


®° F. Morgan and M. Muskat, J. App. Phys. 9, 393 (1938). 
7 L. Giimbel, Zeits. f. tech. Physik 3, 94 (March, 1922). 
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tion has been made by Peppler® in a study of the 
lubrication between the outer surfaces of a pair of 
cylinders—e.g., a pair of rollers—and he too finds 
good agreement between theory and observation. 
If we compare the two problems—the flow 
between the spheres and that in the cylindrical 
bearing—we see that they have one feature in 
common, namely that they both deal with the 
flow in a channel of which a large part is diver- 
gent in the direction of flow; and as it is a general 
result of the classical theory of viscous flow that 
negative pressures are associated with divergent 
flow, we have here the fundamental reason for 
our difficulty. The essential feature of Morgan 
and Muskat’s solution is the introduction of the 
source, and clearly the conditions of the present 
problem do not permit this; hence we are led to 
try the application of Giimbel’s empirical treat- 
ment. Here the main flow is along the y axis and 
we may suppose that the film breaks from the 
moving surface along the x axis (since this is the 
position of narrowest section); in polars, along 
the lines 0=z, 0=2z. Further, relying upon 
Giimbel’s results, we may choose from the condi- 
tions I, II, III the one which gives the simplest 
solution, and we find that this is condition I. We 
therefore take the boundary conditions as— 


p=0 for 6=2r; O<r<nro, 


6=7, 
r=To, 


(11a) 


where ro is yet to be determined. Thus we investi- 
gate the solution of the pressure equation over a 
semi-circular region of the xy plane, and effec- 
tively neglect any contribution to the pressure 
from the flow in the diverging region. 


3. SOLUTION OF THE PRESSURE EQUATION 


The actual solution of the Reynolds equation 
(9) with the boundary conditions (11a) is easily 
seen to be 


Il =g sin 8, (12) 


where g is a function of p only defined by the 
differential equation 


dg dg 
—+L——-—z=M, (13) 
dp? dp jp 
with the boundary conditions 
g =0 when p =(, p= po. 


8’ K. O. Peppler, V.D.I. Forschungs. 391 (1938). 
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TABLE I. Solutions for g.* 














i Gi X10° G: X10° g X10* Gi X10° G: X10° 2 X108 

0 0 0 0 0 0 0 
0.1 0 10 —10 0 10 —10 

2 1 19 —18 1 19 —18 

3 3 28 —24 3 28 —25 

4 6 36 —29 6 36 — 30 

5 11 43 — 31 11 43 — 32 

6 16 48 —3i1 16 48 — 32 

7 22 53 — 30 22 53 —31 

x 27 57 —29 27 57 — 30 

9 32 61 —28 33 61 —28 
1.0 38 63 —24 38 63 —25 
1.1 43 66 —23 43 66 —23 
1.2 47 68 —20 48 68 —20 
1.3 51 70 —18 52 70 —18 
1.4 55 71 —15 56 71 —15 
1.5 58 73 —14 59 73 —14 
1.6 60 74 —13 62 74 —12 
1.7 63 75 —10 64 75 —11 
1.8 65 77 —9 67 76 —9 
1.9 67 77 —8 69 77 —8 
2.0 69 78 —7 71 77 —7 
2.0 69 78 —7 71 77 —7 
2.2 72 80 —6 74 80 —6 
2.4 75 82 —5 76 82 —6 
2.6 77 83 —4 79 83 —4 
2.8 79 84 —3 81 84 —3 
3.0 80 85 —3 82 85 —3 
3.2 82 87 —2 85 87 —2 
3.6 84 88 —2 86 88 —2 
3.8 85 89 —?2 88 89 —1 
4.0 86 90 —2 89 90 —1 
4.2 87 91 —1 90 91 —1 
4.4 88 91 —1 90 91 —1 
4.6 89 92 —1 91 92 —1 
4.8 90 92 0 92 92 0 
5.0 91 93 0 93 93 0 


| 


_ *g is defined by g =Gi+AG: with d chosen so that g =0 when p =5.0. 
rhe three-decimal values given above have been obtained by rounding 
off the original four-decimal values. 


For this function satisfies both the partial differ- 
ential equation and the boundary conditions, and 
is therefore a possible solution. The general 
uniqueness theorem then tells us that this is the 
only solution. 

All that remains now is to find the solution of 
Eq. (13), and although there seems to be no 
simple closed analytical form for the function 
g(p) a numerical solution is very easily con- 
structed. 

This function can be expressed in the form 


g =Gi+ G2, 
where G, is any solution of the equation 
1 
Gi" +LGy——Gi=M (13a) 
p 


(dashes denoting differentiation with respect to p), 
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and G: is a suitably chosen complementary func- 
tion, that is a solution of 


G2!'+ LG,’ ill =0, (13b) 
p? 
and \ is a numerical constant to be chosen later. 
With the forms of L, M given by (9a) and (9b) 
it is readily seen that: 
(a) The only solutions for G; and G: which are 
finite at p=0 must vanish at this point ; near p=0 
they have the forms 


Gi=be[1—Fo?+-- +], 
G2=p[1—{p?+---]. 


(6) Neither G; nor Gz can vanish for any other 
value of p. Therefore if g=G:+AG_2 we have at 
once g=0 for p=0, and we can find \ so that g=0 
for p=po [in fact, \= —Gi(po)/Ge(po) ]. The 
problem is now reduced to the integration of 
Eqs. (13a), (13b) with the conditions 


Gi=0, Gi =0 at p=0, 
G2=0, Ge2’=1 at p=0, 


which can be carried out by standard technique. 

Since there is a choice of two forms for the 
functions L, M two distinct integrations were 
made, one with the approximate forms of (9c), 
(9d) and one with the ‘‘exact”’ forms of (9a), (9b) 
taking hy/a=1/200—this should rank as rather 
a large value, as we shall see later. These numer- 
ical results are given in Table I, and it is at once 
clear that the error of the approximation is 
negligible. But the most important property of 
the functions G; and G2 revealed by the numerical 
solutions is that they quickly attain virtually 
constant values; more particularly, their values 
remain sensibly constant for p> 5. It follows that 
for any given value of \ the function g=G,+AGz 
is itself sensibly constant for p>5. Therefore if 
we take for po any value not less than 5, we shall 
have, with \= —Gi(po) /G2(po), g=90 for p= po. In 
other words, although the quantity po is not 
formally defined_its choice need cause no diffi- 
culty. The value chosen for the solution of 
Table I was po=5.0. 

A little reflection shows that this lack of 
sensitivity to the exact position of the boundary 
is to be expected on physical grounds. For in the 
immediate vicinity of the region of closest ap- 
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proach of the two spheres—which is the only 
region we consider—the pressure must fall off 
very rapidly with increasing distance from the 
actual point of closest approach, and therefore so 
long as we ensure that the boundary is in the 
region where the pressure has fallen to a small 
value, the exact position cannot have any great 
effect on the solution. 


4. TOTAL LOAD SUPPORTED BY THE FILM 


Since the pressure p does not vary through the 
thickness of the film the total load supported is 
simply the integral of » over the xy plane. 

Thus, 


load P= J f pdxdy = f f prdrdé 
® 0 
, Qa} 2e Apo 
= $uw sin «(—) 2a f f IIpdpdé 
ho r 0 


= 12Ka*yw sin a(2a/ho)' since Il=g sin 8, 
where 


po 
K= -f p2(p)dp. 


From the solutions given in Table | we find 


for ss=1/200, K =0.057, 
so=0, K =0.053 (the approximate case). 


Clearly there is no point in keeping many 
decimals in a quantity such as K, so we may take 
as a round value K=1/20 and write the final 
result as 


P= 3a?pw sin a(2a/ho)!. (14) 


It is from this equation alone that we derive 
our criterion for the range of validity of the 
classical theory ; but it is convenient to postpone 
this until we have obtained the corresponding 
relation for the torque exerted by the viscous- 
frictional forces, and have applied both results to 
the particular orientation of spheres found in the 
four-ball instrument. These topics are dealt with 
in the next two sections. 


5. TORQUE ON THE STATIONARY SPHERE 


As was explained in the Introduction it was 
considered desirable to obtain the expression for 
the torque exerted on the stationary sphere in the 
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four-ball apparatus, and this section gives a 
derivation of the corresponding relation for the 
pair of spheres considered in the preceding sec- 
tions. It cannot be claimed that the result ob- 
tained here gives anything more than an indi- 
cation of the leading term of the required relation. 
For on the one hand we can take into account 
only the conditions in the region where the flow 
is converging, since we know nothing of the 
conditions in the diverging region; and on the 
other, the unknown dimension ro (of which the 
expression for the load was found to be effectively 
independent) appears explicitly in the final ex- 
pression. Nevertheless the result should show the 
essential features of the relation. 

The problem is to find the viscous-frictional 
torque on a curved surface immersed in a body of 
liquid whose motion is essentially three-dimen- 
sional ; this is necessarily rather complicated, and 
after various attempts to construct an approxi- 
mate solution it was decided that the only 
satisfactory method was to use the fundamental 
tensor-relations for the shear-stress in the liquid. 
Only the final expression for the torque is needed 
here, and it would serve no useful purpose to give 
at this stage the rather heavy algebra of the 
derivation; but as no explicit treatment of the 
general case was found in the standard textbooks 
the analysis is given in an Appendix. 

If fT is the torque about the axis of rotation 
exerted by the fluid on the stationary sphere, we 
find the approximate relation 


a 
P/33a* po = 1+ {loe. (—) _ 1 sin? a, 
ho 


and if we put log, (a//o) =2.3 logio (a/ho), 
[/$ra* yw = (4 — sin? a) 


a 
+2.3 sin? a logio (~). (15) 


0 


6. APPLICATION TO THE FOUR-BALL 
INSTRUMENT 


To derive from the foregoing results the ex- 
pressions for the load and torque for the four-ball 
instrument we have only to determine the angle 
a between the axis of the rotating sphere and the 
line joining the center of this sphere to the center 
of any one of the fixed spheres. 
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Fic. 4. Geometry of the 4-ball instrument. 


The geometry of the system is shown in Fig. 4. 
It is clear that the vertical axis of rotation must 
pass through the centroid G of the triangle 
formed by the centers Ci, Co, C3 of the fixed 
spheres. 

Then 


and so, from the triangle C,GC,, 





_ 2v3 a 
sin a=GC,/C\Cr=—a /2(a+h) = (4) ° 


at+ho 
Since hoa, we may take 
sin a=(4)!, cosa=(3)! (a=35° 16’). (16) 
The load P determined in §4 acts along the line 
of centers C,C,; by symmetry therefore the total 
load supported by the three spheres has zero 


horizontal component and the vertical compo- 
nent is 3P cos a. Thus the total load on the axis 
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is, by (14) and (16), 


6 2a\! 
W =3P cos a=-0'ns(—) (17) 


] ho 


The total torque on the axis is, by (15) and (16), 


2 a 
T=3r= tra uo| 0.17 log io (~)| (18) 
1 


q lo 


lf now we eliminate a/ho from these relations, 
we get the relation between T and W 


2 6v2 
T/$2a*pwo— = ——+1.54 logue| 1 —a'ne | 
15 5 
(19) 


showing that if experimental values of T and 
W are found and T/a*yw is plotted against 
log (W/a*uw), then the points should lie on a 
straight line. 

The sizes of the quantities involved are best 
shown by a numerical example. In applying the 
relations we must remember that all quantities 
are presumed to be given in absolute units, so 
that if a is measured in centimeters, w in radians 
per sec., and yin c.g.s. units (poises) then W will 
be given in dynes and T in dyne-centimeters. 

The usual size of instrument corresponds to 
spheres of 1-cm diameter, so a=}; for a rather 
viscous lubricant we may take yw=10 c.g.s. 
units at moderate temperatures. Let us take 
w=N r.p.m.=2rN/60 radians per sec. and 
ho/a =10-" (so that hyo =}3-10-* cm). 

A normal speed of rotation is 1500 r.p.m. and 
with this value we find (rounding off the nu- 
merical coefficients) 


W=500X 10"? dynes=0.51 X10" g wt., 
T =720n —125 dyne cm=0.74n —0.13 g cm. 


By taking different values for m we obtain a 
series of related values of load and torque corre- 
sponding to a series of film-thicknesses io; the 
results for the range 2<n<6 (i.e., for film- 
thicknesses ranging from 5X 10-* to 5X 10-7 cm) 
are given in Fig. 5. 

It is now clear that no great loads can be sup- 
ported by the apparatus if the conditions of 
film-lubrication are to be fulfilled. There is a 
lower limit to the thickness of the film, fixed by 
the size of the irregularities of the surfaces, for 
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which a value of the order of 10-4 cm has been 
suggested. With this value and even the rather 
viscous oil considered above, the load W cannot 
greatly exceed 50 g wt.; even if the speed were 
raised to 3000 r.p.m., the permissible load would 
rise only to about 100 g wt. Evidently a load of 
the order of a kilogram represents a virtual im- 
possibility for an instrument of this size under 
conditions of film-lubrication, unless it can be 
established that this regime can be maintained 
with film-thicknesses of the order of 10-* cm. It 
is quite certain that under the usual conditions 
of operation—with loads of 20 to 1000 kg—there 
can be no hydrodynamical lubrication; this is 
of course only what was expected, but the 
investigation has provided the criterion from 
which we can deduce the range of conditions in 
which the behavior of the apparatus obeys the 
laws of hydrodynamical lubrication. 


APPENDIX. TORQUE ON THE 
STATIONARY SPHERE 


The stress at any point in a viscous liquid has 
six independent components which form a sym- 
metrical tensor of the second rank. This tensor is 


V=—pl+uzT, 


where J is the unit tensor and T is the shear- 
stress tensor, and these can be conveniently 
written in “dyadic” form 
J =ii+jj+kk, 
; Ou dw dv 
T=2> —i+ZL{ —+— }Gk+k)), 
Ox Oy dz 
where ijk are as before the unit vectors along the 
coordinate axes.° 
The force on a surface of unit area whose 
normal has the direction of the unit vector n is 
itself a vector and is given by the scalar product 
of n with the stress tensor in the form 
F,,=V-n=—pn+uzT-n. 
Now if a force F is applied at a point P whose 
position vector with respect to a fixed point M 
is r, the moment of F about M is rAF; and if 
through M there passes a fixed axis whose direc- 
tion is the unit vector t then the torque about 


<= Weatherburn, Advanced Vector Analysis (G. Bell 
and Sons, London, 1928), p. 149. 
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Fic. 5. Numerical values for load (W) and torque (T) for 
varying minimum film-thickness (ho). 


this axis is t-rA F. Thus if for F we take the 
hydrodynamical force F,, we have the general 
result that the torque y about the axis t of the 
force exerted by the liquid on the unit area at 
the point r is 


y=t-rAF,=t-rA(—pn+uT-n), 


where it is understood that the pressure p and 
the stress-tensor T are evaluated at the point r. 

The resultant torque [ of the force exerted 
on the surface S in the liquid is then simply the 


integral of y 
r= {ws 


taken over the surface S. 

This result is perfectly general; in any par- 
ticular case we have only to substitute the 
appropriate values for t (given by the position 
and direction of the axis), r, n (given by the form 
and position of the surface), and p, T (given by 
the velocity distribution of the flow). 

We are interested in the torque, about the 
axis of rotation of the moving sphere, of the 
force exerted by the liquid on the fixed sphere; 
we therefore choose for t, r, n the following values 
(Fig. 2). 

MC, is the axis of rotation of the moving 
sphere and P;: the general point (x, y, —/) on the 
fixed sphere. 

Then 


i.€., 


t=sin ai+cos ak, 


r=MP,=[x+(a+hpo) tan a ji+yj—hk. 


The surface being spherical, the normal at any 
point has the same direction as the radius drawn 
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to that point; therefore, 


an =C.P,=xi+yj+Zk. 
Evaluating the triple product we find 


¥=(vF.4+-ZF,) sin a 
+(xF,—yF,) cos a+2hF, sin a, 


where F,, F,, F, are the components of F,,. 

In order to evaluate the tensor components, 
we need explicit forms for the pressure p and 
the velocity derivatives du/dx, etc. The Eqs. (2) 
and the continuity Eq. (3) together determine 
these latter; all except (dw/dx), (dw/dy) are 
given immediately, and for these we have first 
to determine w by integration of dw/dz with the 


condition w=0 when z= —h, and then to differ- 
entiate these with respect to x and y in turn. 
Putting z= —h in these expressions then gives 


the values required. Since 


aF ,,=an-¥ = —(xi+yj+Zk)p 
+ u(xi+yj+Zk) -T, 
= F,i+ F,j+ Fk, 
we can now express the force components as 
functions of x and y and hence also the torque 7; 


the algebra is straightforward but rather heavy, 
and the final result is 


(a) 
2Zh 
—¥ = (x?+-y?) cos? a+2(a+ho)x cos @ sin a 
pwd 


3 h 
+| (r+x" +(1 -—)y+2n2, sin? a 








(a) 
2h?/ ap ap 
-—(x—-y —) COS @ 
uw \ dy dx 
(b) 
h? ap ap dp 
- |»(=—+9—) 422°] sin a 
pwZ Ox oy oy 
(b) (b) 
4h® ap 4h‘Zy 
——— — sin a+ (V2p) sina (A) 
uw Oy a* ww ; 


éeye 


(the indications ‘“‘a’”’ and “‘b” will be explained 
later). 
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We have now to integrate this expression over 
a suitable area of the spherical surface; if we 
introduce the plane polar coordinates r, @ as in 
$2 we have for the element of area the expression 


dS =rdrd6/(1—r*/a?)!, 


and the limits for the integration are r<6<2z, 
O<r<ro. Alternatively we can use a pair of 
angular coordinates 0, ¢, where @ is defined as 
before and ¢ by the relation r =a sin ¢. The limits 
are now r<0<27, O<@< do, where 7o=a sin do, 
and the element of area is 


dS =a’ sin ¢déd¢. 


The first set of coordinates is more convenient 
for dealing with the terms in y which involve p 
explicitly, and the second for the other terms. 
We have also the relations 


Z=athy—h=(a?—r*)'=a cos ¢, 


r=p(2ahy)} 
2a\} 
p= Sw sin «(—) g(p) sin 8, 
ho 
the function g being given numerically (Table I). 
Examination of the terms of (A) shows that 
those marked (a) vanish identically on integra- 
tion and that those marked (b) lead to expressions 
which are of higher order in A/a than the re- 
mainder, and which may therefore be neglected 
since this is a small quantity and we want only 
the leading term in the expression. Thus we have 
to integrate only the terms 


~nandee, ——-+2}sin?a, ——sina. 
hZ pw Oy 


x?+-y? (— ) 2h ap 


The integration is straightforward and the re- 
sult is 


['/2a* pw = 1 —cos d9—log cos do 
+sin? a[2(1—cos do) 
+3(1—cos ¢o—log cos ¢o) 
+log (1+59—cos $0) /so cos do | 


po 
+0O(5so log so) —12 sin? af p’gdp, 
0 


where all logarithms are to the base e. 
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The leading terms of this expression can be 
written 


I'/}2a* pw = A + (log a/ho+B) sin? a, 
where 
A =1—cos ¢o—log cos do 


B=2(1—cos $9) +4(1—cos ¢9—log cos $0) 
po 
+ log [(1—cos ¢o) /cos $o | — 12f p’gdp. 
0 
From the values of g(p) in Table I we find 
po 
_ 12 f p’*gdp = —12X(—0.098) =1.2 approx., 
0 


but we cannot hope to determine the values of 
the terms involving cos ¢@) with any precision, 
since we have no precise knowledge of the value 
of @o. However we can obtain a rough estimate 
by the following argument. ¢o is defined by 
a sin ¢9=7o, where ro measures the extent of the 
region to which the film-theory may be applied. 
This dimension is determined by the thickness h 
of the film, and there is the relation 


h=at+hy—(a?—r*)'=a+hy—a cos ¢. 


If we can say that the theory can be applied 
for thicknesses up to a value h,,, then we can 
take for ¢o the value given by 


cos do = 1— (tm — ho) /a. 


We cannot say that this maximum permissible 
thickness has a definite value, since the value 
must depend upon the particular conditions pre- 
vailing; but we may not unreasonably suppose 
it to be considerably greater than the greatest 
value (10-* cm) of ho considered here, so that @o is 
effectively independent of ho, and cos do=1 
—h,/a. \f as a working hypothesis we assume 
that h,, is not less than about a millimeter, we 
have the figure 0.1 as an indication of the order 
of lower limit of h,,/a (since the radius a@ is of 
the order of unity) and 0.9 for the upper limit 
of cos do. 
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There is also the relation sin ¢)=70/a, and the 
analysis given in the preceding pages makes 
considerable use of the assumption that ro is 
considerably less than a (since the exact equation 
for the pressure has a singularity at r=a). We 
may suppose 7o/a<}, and hence it is suggested 
that cos o> 43(3)!'=0.866. 

These arguments are very specious, but they 
suggest that the order of the constants A, B will 
be given by taking a value in the region of 0.9 
for cos @o; if we tabulate the values for a few 
values of cos do we get 





COS go | A B 
0.95 | 0.1 —1.7 
90 | 0.2 —0.9 
85 | 0.3 —0.3 
80 | 0.4 +0.1 


so we shall not be seriously in error if we choose 
the round values 


A=}, B=-1. 


Thus we have the result that the torque is a 
linear function of log a/ho having approximately 
the form 


a 
P/32a* po = 1+| lor. (=) — 1 sin? a. 
0 


This indicates a very slow variation with a/hy 
and it is for this reason that the expression for 
y was evaluated in full; it is difficult to ensure 
that a “leading term’’ type of argument applied 
early in the calculation does not neglect terms 
which are in fact significant. 
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This is a general method of applying infra-red absorption spectroscopy to the analysis of 
multicomponent hydrocarbon mixtures in the liquid phase. The procedure is calibrated by 
measuring optical densities of synthetic standard samples. A constant thickness cell is used for 
both standard and unknown samples. Two procedures for converting optical densities to con- 
centrations are described. Examples are given for four- and five-component mixtures. Data are 
presented to show the reproducibility of repeated measurements of optical density on the same 
sample, and results are given to show the agreement between the infra-red analyses and the 


known composition of synthetic mixtures. 


INTRODUCTION 


GENERAL method of applying infra-red 

absorption spectroscopy to the analysis of 
multicomponent hydrocarbon mixtures in the 
liquid phase has been developed and has been 
applied to two particular mixtures, one con- 
taining four components and the other five com- 
ponents. It is an empirical method based upon 
analytical curves prepared from absorption 
measurements of pure samples of the com- 
ponents to be found in the unknowns and of 
binary mixtures of these components. 

A number of papers on the application of 
infra-red absorption spectroscopy to analytical 
problems have been published. Most of these 
have been concerned with the analysis of binary 
mixtures. However, two papers have described 
methods for the analysis of more complicated 
mixtures. Brattain, Rasmussen, and Cravath! 
have described a method for analyzing multi- 
component mixtures in the vapor phase. This 
method is most useful when the components are 
in the vapor phase under normal conditions of 
temperature and pressure, but can also be used 
with considerably more difficulty when the com- 
pounds are liquids under normal conditions. A 
method for analyzing a multicomponent mixture 
in the liquid phase has been described by Nielsen 


1R. Robert Brattain, R. S. Rasmussen, and A. M. 
Cravath, J. App. Phys. 14, 418 (1943). 
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and Smith? and the analysis of a three-component 
mixture has been discussed as an example. This 
method can be used only over concentration ranges 
for which deviations from Beer’s law are negligible. 

The procedures described in the present paper 
are used for multicomponent mixtures in the 
liquid phase and are applicable even though the 
concentrations of any of the components of the 
mixture vary from 0 to 100 percent. 

The absorption of monochromatic radiation as 
it traverses an absorbing liquid is expressed by 
Beer’s law. 

acx =logio (Jo/T), (1) 


where 7 is the intensity of the radiation trans- 
mitted by x centimeters of liquid containing a 
concentration c of the absorbing materials; J is 
the intensity of the radiation impinging on the 
liquid; and a is a constant of proportionality, 
known as the absorption coefficient. a is constant 
for any monochromatic radiation regardless of 
the values of c, x, or Jo. 

For a mixture of two or more non-associating 
components, all of which selectively absorb 
somewhat over the whole infra-red region, Beer’s 
law takes the following form for monochromatic 
radiation : 


I 
(iC, +4 2C2+a3¢3+--- +an¢n)x=logie—, (2) 


2 J. Rud Nielsen and Don C. Smith, Ind. Eng. Chem. 
Anal. Ed. 15, 609 (1943). 
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where a, is the absorption coefficient of material 
(1) and c; its concentration in the cell of thick- 
ness X; G2, @3°--d, and Co, C3---C, are the corre- 
sponding values for materials 2, 3-- -n. 

In actual practice, however, Eq. (1) and, 
therefore, Eq. (2) does not hold exactly. This is 
shown experimentally by the fact that a curve 
of log I)/I as a function of concentration with a 
constant thickness cell is not a straight line. 
This results from the fact that no spectrometer 
produces a monochromatic image at the re- 
ceiver. Instead, there is produced an image con- 
taining a wave-length band Ax. 

Equation (3) is obtained by combining Eqs. 
(1) and (2). 


Io 


, Io , fe o, 41 (3) 
og —=log —+log —+---+lo 
A Mines Ses *} 


n 


The terms 


ot wie —... 
og—, log—, ---log— 
£ I, g Ie £ I, 


are the values obtained for materials 1, 2---n 
at the chosen spectral position. In other words, 
the equation states that the optical density of a 
mixture is equal to the sum of the optical 


























Fic. 1. An assembled cell used to obtain the infra-red 
absorption of liquids. 
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Fic. 2. Unassembled parts of a cell used in measuring the 
infra-red absorption of liquids. 


densities of the individual components of the 
mixture. Equation (3) does not contain the ab- 
sorption coefficient and, therefore, is not subject 
to the same limitations as are (1) and (2). 


DEVELOPMENT OF ANALYTICAL PROCEDURE 


The problem in the work described here, as 
well as in the work of others, has been to de- 
velop a method of applying absorption spec- 
troscopy taking into account these limitations of 
Beer’s law. The fundamental problem of quan- 
titative analysis in absorption spectroscopy is the 
determination of the log J)/I values with suf- 
ficiently high accuracy, and the conversion of 
these values to percent concentration of the 
components in the mixture. ' 


Absorption Cells 


In the present work, optical densities are 
measured with a University of Michigan infra- 
red spectrophotometer* ‘ equipped with a sodium 
chloride prism. Two cells are used. The back- 
ground cell, consisting of a single salt plate 3” 
thick, is used to determine J» values. The absorp- 
tion cell, shown in Fig. 1 and Fig. 2, is a modi- 
fication of a cell described by Randall. The 
frame which holds the salt plates is designed par- 
ticularly for use with the University of Michigan 
instrument. 

Starting from the left, each piece in Fig. 2 
goes on top of the previous one. The salt plates 
of the absorption cell are }’’ thick. They are 
separated by a platinum shim, the thickness of 


3H. M. Randall and John Strong, Rev. Sci. Inst. 2, 585 
(1931). 

‘Robert A. Oetjen and H. M. Randall, unpublished 
thesis, University of Michigan, 1941. 

5H. M. Randall, Rev. Sci. Inst. 10, 195 (1939). 
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Fic. 3. A record of the interference fringe pattern obtained 
with a cell 0.152 mm thick. 


which determines the thickness of the absorbing 
layer of sample. For the work described here the 
shim is 0.152 mm thick. The positions of the salt 
plates and the shim are determined by the dowl 
pins in the cell supports. Two rubber washers 
separate the salt from the stainless steel housing. 
After the cell is assembled, a mercury seal is 
established by introducing mercury into the 
trough in the lower plate through a hole which 
extends through the upper salt plate, rubber 
washer, and support. The liquid to be examined 
is introduced through one of the two holes on 
the side of the upper plate. Rubber pads are 


_pressed over the holes to seal them without 


putting pressure on the contents of the cell. This 
method of sealing makes it possible to use the 
cell in an evacuated source box. The cell can be 
taken apart, cleaned, and reassembled in less 
than 10 minutes. However, it can be refilled in 
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less than two minutes by expelling the liquid 
with an air syringe and filling with an eye 
dropper. 

In order to use the present method success- 
fully, it is important that the thickness of the 
absorption cell remain constant. The two salt 
plates are ground flat with corundum in ethyl 
alcohol, and polished on a piece of white felt 
stretched over a piece of plate glass. Titanium 
oxide is used as the abrasive for polishing, and 
ethyl alcohol is used to wet the titanium oxide. 
The finished plates are examined with an optical 
flat and sodium light, and if they are not flat 
within four diffraction fringes (i.e., 0.0012 mm) 
the process is repeated. The plates are examined 
frequently to determine whether cleaning has 
caused deviation from flatness. 

The cells are supported in a tube directly in 
front of the salt window at the entrance slit of 
the spectrometer. The notch in the frame of the 
cell locates it in a given fixed position in front of 
the slit. The cells are handled by means of a 
knurled rod screwed into the upper plate. 

The thickness of the cell is determined by the 
interference fringe method as described by Smith 
and Miller.* (See Fig. 3.) Repeat measurements 
made several months apart show the thickness to 
be 0.152+0.0005 mm. 


Spectral Positions 


The approximate spectral positions at which 
optical densities are determined are selected from 
the records of the spectra of the pure components 
in the region from 3 to 13.5u. (See Fig. 4.) 

In order to obtain high resolution for these 
records the spectrograph is operated with narrow 
slits and slow rotation of the prism. The source 
of light is a. Nernst glower. The light beam is 
interrupted by means of a shutter having a 
period of 3.8 seconds, and the output of the 
thermopile is amplified by a Firestone amplifier.’ 
The recording galvanometer is five meters from 
the camera drum. The spectra are recorded in 
small sections with constant slits, and an over- 
lapping of these sections of spectrum occurs 
whenever the slit widths are changed. Table I 
gives the maximum theoretical resolution ob- 

6 Don C. Smith and Elmer C. Miller, J. Opt. Soc. Am. 


34, 130 (1944). 
7F. A. Firestone, Rev. Sci. Inst. 3, 163 (1932). 
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Fic. 4. The absorption spec- 
trum of 2,2 dimethylbutane in 
the region 2.5 to 13.5. 





tainable with the slit widths necessary to get 
half scale deflections with this amplifier. Ac- 
tually the resolution realized with the instrument 
is somewhat less than the theoretical value. 

In the analysis of an m component mixture it 
is necessary and sufficient to determine log Jo/J 
values of the mixtures at n different spectral 
positions. When a procedure is being developed 
for a given mixture, the spectra of the pure 
components are examined and m positions are 
selected in such a way that a different component 
contributes a maximum amount of absorption at 
each position. It is desirable but not necessary 
that at each wave-length position the absorption 
of one component be greater than that due to 
any other component; and further, that the com- 
ponent showing major absorption be different 
at each of the n wave-length positions. The 
spectral positions which yield most accurate 
analyses are those at which the absorption of one 
component is large compared to the sum of the 
absorptions of all other components (assuming 
all components are present in concentrations 
similar to those in an unknown mixture). In 
general, this condition is found at the position 
of maximum absorption of a band for one com- 
ponent. 

There are some other factors to be considered 
in the selection of these wave-length positions. 
For example, in order to reduce the time neces- 
sary to make an analysis, it may be desirable to 
select bands near each other. It may be desirable 
to select positions which are outside of the water 
vapor region. If one thickness of the absorption 
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cell is to be used, the selection of the wave-length 
positions may depend upon the intensities of the 
bands; however, as far as possible, thinner ab- 
sorption cells should be used when the bands are 
strong, and thicker ones should be used when 
the bands are weak. The selection of these 
positions is often complicated by the fact that 
several bands interfere with one another. This 
difficulty can be minimized by using the nar- 
rowest slits permitted by the sensitivity of the 
receiver and the available amplification of its 
output. 


Determination of log I,/I Values 


For the measurement of log Jo/J values, the 
thermopile in the spectrometer is connected 
directly to a galvanometer with a sensitivity of 
one-tenth microvolt per millimeter at one meter. 


TABLE I. Theoretical resolution of the spectrograph. 

















d du 
Actual Effective - B 

aN slit widths* slit widths dd dd 
Qu 0.04 mm 0.0097 208 640 
3 0.04 0.0117 257 440 
4 0.05 0.0128 314 560 
§ 0.06 0.0118 424 640 
6 0.07 0.0117 514 800 
7 0.08 0.0117 600 960 
8 0.09 0.0117 686 1080 
9 0.11 0.0127 711 1160 
10 0.13 0.0130 768 1240 
11 0.16 0.0140 787 1560 
12 0.20 0.0163 738 1720 
13 0.24 0.0183 711 1800 
14 0.35 0.0238 589 2080 
15 0.55 0.0339 443 2240 

* These slit widths produce half scale deflection. 
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TABLE II. Resolution for operating conditions. 


Spectral Actual slit Effective slit Actual Prism 
position widths width (yu) resolution resolution 
(pm) (mm) (dd) d/dd Bdn/dd 
8.23 0.55 0.071 116 1100 
8.55 0.60 0.075 114 1120 
8.88 0.70 0.092 96 1152 
10.50 0.90 0.085 123 1400 
11.28 1.10 0.096 117 1600 
9.13 0.7 0.080 114 1168 
10.31 0.9 0.088 117 1336 
11.89 1.2 0.100 119 1704 
12.30 1.3 0.116 107 1752 


The _ scale-to-galvanometer distance for the 
present work is five meters. With this system it 
is necessary to use wider slits and therefore to 
accept lower resolution than that realized with 
the amplifier, but the accuracy of measurements 
is much better. The actual resolutions for the 
slit widths used are shown in Table II. 

When the procedure for a given mixture is 
developed, the slit widths for each spectral 
position are set to give approximately 400-mm 
deflection when a non-absorbing material (back- 
ground cell) is in the light path. Of course, after 
the necessary calibrations have been made for 
the analysis of a particular mixture, the same 
slit widths must always be used for this analysis. 

After the proper slit widths have been selected, 
the spectral positions are re-examined under the 
conditions which would be used in making the 
analyses, that is with wide slits and very slow 
rotation of the prism. These records for the 
hexanes are shown in Fig. 5. The upper line on 
each of the records shows the deflection of the 
galvanometer for the background cell. The 
bottom lines in the figures show the zero posi- 
tions, and displacements along these lines in the 
vertical direction represent zero shift. The ab- 
sorption curves lying between the background 
curve and the zero lines are as designated. Cor- 
rect log Io/J values can be obtained only if (1) 
the deflection of the galvanometer is a known 
function of the energy striking the thermopile 
and (2) all energy reaching the thermopile lies 
in the selected wave-length interval. In the 
present case, all tests have shown the galva- 
nometer system to be linear. However, the 
energy reaching the thermopile actually includes 
some energy in addition to that from the wave- 
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length interval, AX. This wave-length interval is 
determined by the slit width and the region of 
the spectrum for which the spectrometer is set. 
The radiation in this wave-length interval which 
passes through the absorption cell and strikes 
the thermopile is the radiation which should be 
measured. Any other energy which causes a 
change in the galvanometer deflection (or other 
measuring instrument) when the shutter is 
opened is called false energy and its effect must 
be taken into account. 

Most of this false energy is probably due to 
scattering from various parts of the spectrometer 
of radiation in the region 1 to 4u which is the 
region of maximum intensity incident on the 
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Fic. 5. Portions of the spectra of the five hexanes. 


Curve Position 
no. Component of band 
il n hexane 11.28 
12 3 methylpentane 10.50 
13 2,3 dimethylbutane 8.88 
14 2 methylpentane 8.55 
15 2,2 dimethylbutane 8.23 
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entrance slit. This energy interferes mostly at 
long wave-lengths (6 to 15u) where the true 
energy is small and the slits are wide. When 
working above 7 most of the false energy can 
be cancelled out by using a mica shutter which 
transmits most of the radiation up to 5y but 
none above 7u. In this case the radiation below 
5u always passes into the spectrometer and is 
therefore included in the ‘“‘zero’’ reading. Only 
the radiation of wave-lengths greater than 7y is 
interrupted by the mica shutter. However, this 
method does not cancel out all the false energy; 
therefore, more refinement is necessary. A 
method of making partial correction for the re- 
maining false energy has been established by 
using a piece of Cellophane which has two regions 
of 100 percent absorption, one at 9.154 and one 
at 11.904, and obtaining the difference in the 
galvanometer deflections when the cells are in 
the light path with, first the Cellophane, and 
second the mica shutter. The amount of false 
energy when the background cell is in the light 
path is different from the amount when the 
liquid cell is in, but is practically the same for 
all liquids used in the liquid cell. It is propor- 
tional to Ip at any given wave-length and there- 
fore can be expressed as a percent of Jo. The 
false energy values at the two bands in the 
spectrum of the Cellophane film have been 
plotted on a graph showing wave-length versus 
false energy and a straight line has been drawn 
between the points. (See Fig. 6.) Curve (1) is 
for the background cell and curve (2) is for the 
liquid cell with any of the nine pure components 
in it. At every spectral position and given slit 
width, J» is made the same for all observations. 
Therefore, the false energy corrections can be 
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Fic. 6. False energy as a function of wave-length. 


converted to mm deflection and subtracted from 
the J and Jy values obtained, as shown in Table 
Ill. 

With the equipment used in this work there is 
some zero drift and the intensity of the light 
source is not perfectly controlled. For this reason 
the procedure in making a measurement at a 
given spectral position consists of making three 
readings as follows: (a) with the absorption cell 
alone in the light path, (6) with the cell plus a 
mica shutter in the light path, and (c) with the 
absorption cell again alone in the light path. 

The averaging of (a) and (c) in general cancels 
out the zero drift and the subtraction of this 
average from (b) gives (k). This number is 
proportional to the difference in the energy which 
falls on the thermopile with the two shutter 
positions. In other words, (k) is proportional to 
the energy which is interrupted by a mica 
shutter. A similar value (L) is obtained with the 
background cell. Then because the source may 
have changed its intensity during this interval 
of time, a value similar to (k), called (m), is 
again obtained. (k) and (m) are then averaged 
to give (m). This should, in general, cancel out 
any errors due to slow changes in the source. 
(n), therefore, is J plus the false energy F and 


TABLE III. False energy corrections for each spectral position and slit width. 











False energy 
Slit widths ion 
Principal absorber A(u) cm~ (mm.) Ie (mm) Io (mm) I (mm) 
2,3 dimethylbutene— 1 9.13 1096 0.70 320.0 0.5 0.5 
2,3 dimethylbutene—2 10.31 970 0.90 345.0 2.0 1.5 
2,2,3 trimethylbutene 11.89 841 1.30 390.0 3.5 3.0 
Methylpentenes 12.30 813 1.40 375.0 4.0 3.5 
2,2 dimethylbutane 8.23 1215 0.55 450.0 0.5 0.5 
2 methylpentane 8.55 1170 0.60 420.0 0.5 0.5 
2,3 dimethylbutane 8.88 1126 0.70 495.0 0.5 0.5 
2 methylpentane 10.50 952 0.90 430.0 2.0 2.0 
n hexane 11.28 886 1.10 430.0 3.0 2.5 
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Tas_e IV. Standard deviations of log Jo/J values for two 


samples. 
Sample 223 Sample 756 

0.549 0.549 0.153 0.154 

.550 554 .153 154 

.548 545 .153 .150 

.548 547 .152 151 

548 .548 .150 .150 

549 548 -152 .149 

.547 .544 .149 153 

551 .548 155 .148 

.550 .546 33 .147 

552 .548 .145 .150 

551 .552 .149 151 

547 .550 .150 .150 

549 .550 151 .152 

546 551 151 .149 

.548 .548 .149 151 

Average 548 151 

Max. dev. from 

the mean .004 .006 

Std. dev. .0022 .0023 


(L) is Jo plus the false energy Fo. These obser- 
vations are written in a data sheet as follows: 


Absorption cell 


Background cell 





Shutter Shutter Shutter Shutter Shutter Shutter 
open closed open open closed open 
(a) (b) (c) (d) (e) (f) 

(g) (h) (j) 


The values of Jo + Fo and J+ F are obtained from 
the above readings in the following manner: 


r (a) + (c)7 d)+ 
0 {[2F9]-0, @ [2M], 


r(g) + (J) 
(h) — ee l= cm), : 


ie = - - 


(L)=Iot+Fo, (n)=I+F. 


k)+ 
( (m) ‘eli 


’ 








F and Fo are the false energy values which are 
subtracted to determine J) and J (see Table III). 

It is essential in this method that changes in 
optical densities due to unavoidable changes in 
the two cells be known. It is, therefore, desirable 
to have a standard absorbing material which can 
be checked daily. This check is made with a 


* “pure” sample of 2,2,3 trimethylbutene at 11.89, 


at the beginning and end of each day. When the 
log Io/I deviates from the standard value 
(0.600) by more than 0.003, all values of log [o/J 
determined on that day are corrected by adding 
or subtracting the indicated amount. 
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Reproducibility in Determining log I,/I Values 


Thirty determinations of log J)/I values at one 
spectral position for two samples have been 
made over a period of three months. These values 
are listed in Table IV. 

The standard deviations for the two sets of 
data are essentially the same even though the 
amount of absorption differs. For this particular 
case a change of 0.002 in log J)/J represents a 
change of 0.3 percent in the concentration of the 
component sought. For a cell twice as thick it 
would represent 0.15 pereent, for one-half as 
thick, 0.6 percent. Likewise the deviation in 
terms of percent concentration would be more 
or less, depending on whether the principal 
absorber absorbs more or less at this wave-length. 


Analytical Curves 


In making infra-red analyses it is necessary to 
convert values of log J)/I to concentrations. The 
conversion may be made by means of analytical 
curves each of which shows log J)/J as a function 
of concentration of a particular compound. In 
the present work with liquids, the data for such 
a curve are obtained by measuring log Jo/J with 
binary mixtures in the absorption cell. 

The amount of absorbing material in the path 
is varied by using a cell of constant thickness and 
diluting the absorber with one of the compounds 
which has a high transmission at the wave-length 
in question. Three series of binary mixtures were 
made for the four-component mixture and four 
series were made for the hexane mixture. These 
solutions have concentrations as shown in Table 
V. It will be seen that for the four-component 
problem, one of the components of all the mix- 
tures is 2,3-dimethylbutene—2. This was se- 
lected because it absorbs strongly at only one 
of the spectral positions. At all the other positions 
it absorbs very little. In the case of the hexane 
series 2,2 dimethylbutane was selected as the 
common component. 

The preliminary analytical curves are obtained 
by plotting, for each spectral position,.a curve 
showing percent concentration of the material 
diluted, which is 100 percent in the ‘“‘a”’ solution 
(Table V), versus the log Io/I for each of the 
binary mixtures (a through h) in the series. For 
the four-component problem, this means three 
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TABLE V. Compositions of binary mixtures used for analytical curves. 











Mixture Amounts by percent volume 
oy Components a b c d e f g h 
Four component mixtures 

1 Methylpentenes 32 16 8 4 2 0 
2,3 dimethylbutene — 2 0 36 68 84 92 96 98 100 

2 2,3 dimethylbutene— 1 100 64 32 16 8 4 2 0 
2,3 dimethylbutene — 2 0 36 68 84 92 96 98 100 

3 2,2,3 trimethylbutene 100 64 32 16 8 4 2 0 
2,3 dimethylbutene — 2 0 36 68 84 92 96 98 100 

Hexane mixtures 

11 Normal hexane 100 60 30 15 7 4 0 
2,2 dimethylbutane 0 20 40 70 85 93 96 100 

12 3 methylpentane 100 80 60 40 20 10 5 0 
2,2 dimethylbutane 0 20 40 60 80 90 95 100 

13 2,3 dimethylbutane 100 80 60 40 20 10 5 0 
2,2 dimethylbutane 0 20 40 60 80 90 95 100 

14 2 methylpentane 100 80 60 40 20 10 5 0 
2,2 dimethylbutane 0 20 40 60 80 90 95 100 








curves at each of the four spectral positions, and 
for the hexanes, four curves at each of the five 
spectral positions. 

The final analytical curves are obtained from 
these preliminary curves by a procedure which 
will now be described. First, the final analytical 
curve of the component which has the least 
absorption at the spectral position being con- 
sidered is drawn as a straight line from the (0.0) 
point to the point at 100 percent concentration. 
This is illustrated by Fig. 7 which is a graph of 
the preliminary and final analytical curves at 
9.314 for the four-component mixture. At this 
spectral position 2,2,3 trimethylbutene (C, Fig. 7) 
is the least absorbing pure compound. The con- 
tribution of this component to the log Jo/J 
value of the binary mixtures of this compound 
with 2,3 dimethylbutene —2 is determined from 
this final curve for 2,2,3 trimethylbutene. This 
contribution (lower C curve, Fig. 7) is subtracted 
from the preliminary curve of the binary mixture 
of these two components (upper C curve, Fig. 7) 
at a convenient number of points. These dif- 
ferences are the contribution to the log Jo/J 
values of 2,3 dimethylbutene—2 in each binary 
mixture and, when plotted against the concen- 
tration of 2,3 dimethylbutene —2, represent the 
final analytical curve of this component. As 2,3 
dimethylbutene—2 is the common component 
in binary mixtures of each of the other com- 
pounds, the final analytical curves for the 
remaining two compounds can be obtained by a 
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series of subtractions similar to that described 
above. (Subtract difference between upper and 
lower C curves from upper A and D curves to 
obtain lower A and D curves). In this way an 
analytical curve for each of the components in 
the four-component mixture has been obtained 
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Fic. 7. Analytical curves at 9.314. Curves with circles are 
preliminary curves. 


. 2,3 dimethylbutene—1 
. 2,3 dimethylbutene—2 
. 2,2,3 trimethylbutene 
. methylpentene 
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at each of the four spectral positions (16 curves). 
For the hexane mixtures there are five curves at 
each of five spectral positions or a total of 25 
curves. 

After the analytical curves are obtained, it is 
possible to analyze mixtures which contain any 
or all of the components in any possible concen- 
tration. To do this, it is necessary only to obtain 
the log Jo/J values at all the spectral positions 
(nm, if an m component mixture) and to convert 
them to concentrations by use of one of the two 
methods described below. 


Two Methods of Converting Optical Densities 
to Concentrations 


One method of converting optical densities to 
concentrations is a method of graphical approxi- 
mation. In this case the analyses are arrived at 
by successive approximations based upon the 
final analytical curves. (See Fig. 7.) This method 
presupposes that the optical density of the 
mixture is equal to the sum of the optical den- 
sities of the component parts of that mixture as 
expressed in Eq. (3). 

In converting optical densities to concentra- 
tions by this method, the log Jo/J value at one 
spectral position is applied to the analytical 
curve of the material which has its major ab- 
sorption at that wave-length, and the first ap- 
proximation of percent concentration for that 
component is read from the curve. The log Jo/J 
value at a second spectral position is corrected 
by subtracting the contribution of the first com- 
ponent at this wave-length. The corrected 
log Io/I value is applied to the analytical curve 
of the material which has its major absorption 
at this position to obtain the first approximation 
for this component. The log Jo/J value at a third 
spectral position is corrected for the contribution 


‘of the first and second components and the cor- 


rected log Io/I value is applied to the appro- 
priate curve to obtain the first approximation 
for this third component. This procedure is con- 


- tinued until the first approximation is obtained 


for all components of the mixture. Successive 
approximations are obtained by repeating this 
process for each component, always using the 
latest approximate concentration of the other 
components. This is continued until two con- 
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secutive sets of approximations are the same. 
The last set of values is the analysis of the 
mixture. 

The analytical curves can be used graphically, 
as just described, or they can be used to obtain 
mathematical expressions to make the same ap- 
proximations. In this paper the latter will be 
called a method of applied corrections. The 
analytical curves, in general, are not straight 
lines (see Fig. 7). However, a straight line drawn 
between the points on the curve corresponding 
to concentration =0 and concentration = 50 per- 
cent does not deviate far from the analytical 
curve at any point. Each of these straight lines 
is represented by an equation such as (1). The 
slopes of the straight lines are the approximate 
absorption coefficients, so that the optical 
density (log Jo/J) of a mixture at a given 
spectral point can be written from Eq. (2) as: 


(log Io/I)x, = 61.61: +b262+b3¢+ re “bn.Cny (4) 


where );,, b2,---mn, are the slopes of the straight 
lines referred to above for components 1, 2---n 
at wave-length \;, and ¢, c2---c, are the con- 
centrations. There is one of these equations for 
each of the 1, 2, 3---m wave-lengths. 

For the four-component mixtures the equa- 
tions are: 


I 
(oe =) =0.01442c,+0.00270c: 
I Al 
+-0.00840c;+0.00499¢,, 
Io 
(og =) = 0.00193c,+0.00391c. 
he 


+0.00136c;+0.00601c,, (5) 


I 
(tog >) = 0.00155c;+0.00082c: 
I A3 
+4-0.00622c3+0.00757¢s, 
I 
(tog =) =0.00080c;+0.00125¢; 
4 


+0.00092c;+0.02792c,, 
where the X’s and c’s are as follows: 


A: =9.12u cc, is the concentration of 

2,3 dimethylbutene—1, 
Ae=10.314 Cc is the concentration of 

2,3 dimethylbutene — 2, 
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Fic. 8. Correction curves for four-component mixture 
(6) used with method of applied corrections. (See Fig. 7 for 


Io Io explanation of A, B, C, and D.) 
ct=— 14.48 (log =) —12.52 (108 =) 
Ii I, ; : ; . 
Equations (6) are used in connection with a 
wi7(1 Io are7lt Io series of correction curves which are obtained in 
F171 log he ne og h , the following manner: the differences between 


the log (Jo/J) values as obtained from the 

I I as . . 9 

c,= —0.2710 log — =11.69(log = analy tical curves and the straight lines men 
I; Ie tioned above are plotted against the concen- 

trations of each component at each wave-length. 


~3.214( log =) +39.19( log =). For the four-component mixture, these are shown 
T; I, in Fig. 8. 


TABLE VI. Results of analyses of synthetic mixtures of four components. 











2,3 dimethylbutene —1 2,3 dimethylbutene —2 2,2,3 trimethylbuten Methyipentenes 
Sample no. Actual Spec. Actual Spec. Actual Spec. Actual Spec. 

U-1 25.0% 25.2% 25.0% 23.2% 25.0% 25.7% 25.0% 24.9% 
U-2 30.0 30.4 20.0 18.5 30.0 31.4 20.0 19.9 
U-3 20.0 20.1 30.0 29.3 35.0 36.9 15.0 15.1 
U-5 33.3 33.2 33.3 33.9 33.3 33.6 0.0 0.0 
U-6 33.3 34.0 33.3 33.5 0.0 0.0 33.3 33.3 
U-7 0.0 0.0 33.3 34.0 33.3 34.3 33.3 33.6 
U-8 33.3 33.5 0.0 0.0 33.3 34.3 33.3 33.7 
O-1 25.0 24.6 25.0 23.9 25.0 24.4 25.0 25.2 
O-2 12.5 11.5 62.5 63.5 33.5 11.2 12.5 12.3 
O-3 42.0 42.8 18.0 17.3 20.0 19.0 20.0 19.4 
O-4 36.0 35.6 18.0 16.9 26.0 25.1 20.0 20.5 
O-5 33.3 33.6 33.3 31.9 33.3 32.6 0.0 0.0 
0-6 33.3 33.6 0.0 0.0 33.3 33.5 33.3 33.0 
O-8 33.3 33.1 33.3 33.0 0.0 0.3 33.3 33.0 
0-9 50.0 50.7 0.0 0.0 50.0 50.1 0.0 0.0 
O-10 0.0 0.0 0.0 1.1 50.0 49.7 50.0 49.9 
O-11 yf 1.9 3.0 3.0 90.0 88.8 4.5 4.9 
Av. dev. 0.36% 0.72% 0.77% 0.25% 
Max. dev. 1.0 1.8 1.9 0.6 
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TABLE VII. Results of analyses of synthetic mixtures of five components. 











2,2 dimethylbutane 2,3 dimethylbutane 2 methylpentane 3 methylpentane n hexane 

Sample no. Actual Spec. ‘Actual Spec. Actual Spec. Actual Spec. Actual Spec. 
H-1 20.0% 204% 20.0% 19.8% 20.0% . 2% 20.0% 19.6% 20.0% 19.8% 
H-2 10.0 10.4 10.0 10.1 10.0 9.7 10.0 9.6 60.0 59.2 
H-3 5.0 5.4 5.0 aa 5.0 5.0 55.0 54.0 30.0 28.4 
H-4 2.5 2.8 52.5 51.8 a0 aud 27.5 27.1 15.0 15.1 
H-5 1.2 1.6 26.3 25.9 1.2 50.3 13.8 13.1 7.5 7.8 
H-6 50.6 50.7 13.1 13.3 5.6 24.7 6.9 6.8 3.8 4.4 
Av. dev. 0.33% . 0.32% 0.54% 0.50% 0.60% 
Max. dev. 0.4 0.7 0.9 1.0 6 











To make an analysis the log (Jo/Z) values of 
the mixture are determined at each of the 
selected spectral positions and are substituted 
in Eqs. (6). The first approximate concentrations 
are obtained from these equations. The cor- 
rection to be applied to the experimentally de- 
termined log (Jo/J) at each wave-length position 
is determined by subtracting the sum of the 
values obtained from the correction. These new 
log (Jo/J) values are then used in Eq. (6) and 
second approximate concentrations are deter- 
mined. This process is carried on until the con- 


centrations are the same- for two consecutive 
approximations. The last set of values is the 
analysis. In the present work, three approxima- 
tions have been sufficient. 


DISCUSSION OF DATA 


In order to check the methods, synthetic 
mixtures have been made and analyzed. The 
results are shown in Table VI for the four-com- 
ponent mixture and in Table VII for the hexanes. 
It will be seen that in no case does the analysis 


obtained deviate from the actual composition 


TABLE VIII. Comparison of results by two methods. 












































2,3 dimethylbutene —1 2,3 dimethylbutene —2 2,2,3 trimethylbutene Methylpentenes 
Graphical Graphical Graphical Graphical 
Sample approxi- Applied approxi- Applied approxi- Applied approx- Applied 
no. mation correction mation correction mation correction mation correction 
0-3 20.1% 19.8% 29.3% 29.1% 36.9% 37.0% 15.1% 15.1% 
O-4 35.6 35.5 16.9 16.6 25.1 24.7 20.5 20.4 
O-5 33.7 33.6 32.2 31.9 32.5 32.6 0.0 0.0 
O-11 2.0 1.9 2.9 3.0 89.2 88.8 5.0 4.9 
TABLE IX. Data showing the effects of extraneous elements. 
2,3 dimethyl- 2,3 dimethyl 2,2,3 trimethyl- Methyl- 2 methyl Octene 
Sample butene —1 butene —2 butene pentene butene —2 cut 
As made 33.3% 33.3% 33.3% - 0.0% 0.0% 0.0% 
Analysis 33.2 33.9 33.6 0.0 Assumed Absent 
U-5 ' 
As made 32.5 32.5 32.5 0.0 0.0 
Analysis 33.4 35.8 28.1 2.4 Assumed Absent 
As made 0.0 33.3 33.3 33.3 0.0 0.0 
Analysis 0.0 34.0 34.3 33.6 Assumed Absent 
U-7 
As made 0.0 32.5 32.5 32.5 0.0 2.5 
Analysis 0.0 335 33.5 32.6 Assumed Absent 
As made 33.3 0.0 $3.3 33.3 0.0 0.0 
Analysis 33.5 0.0 34.3 33.7 Assumed Absent 
U-8 
As made 31.5 0.0 31.5 31.5 2:5 2.5 
Analysis 35.9 0.0 34.4 34.0 Assumed Absent 
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by more than 2.0 percent of the total mixture. 
The average error for all components in all 
mixtures shown in the two tables is 0.5 percent. 

All the synthetic mixtures have been analyzed 
by the graphical approximation method, and 
some have been analyzed by both methods. 
Results obtained by the two methods for four 
samples are compared in Table VIII. 

The complete analysis of a four-component 
mixture can be made in 45 minutes with the 
method of graphical approximations, and in 30 
minutes with the method of applied corrections. 

It should be pointed out that the presence in 
the unknown mixture of compounds for which 
analyses are not made will cause large or small 
errors depending upon whether the extraneous 
materials absorb strongly or weakly at any of 
the spectral positions. It has been found in the 
case of the four-component mixture that many 
of the unknowns submitted for analysis have 
contained two extraneous compounds. Synthetic 
mixtures including 2.5 percent of each of these 
components were prepared and analyzed, as- 
suming the extraneous components absent. The 
results of this study are shown in Table IX. It 
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will be noticed that these concentrations of the 
extraneous components cause errors up to 4.0 
percent in the analyses. If the approximate 
amount and effect of the extraneous components 
can be determined, and they frequently can, 
errors such as those shown in the example can 
be eliminated or at least minimized by making 
appropriate corrections. 


SUMMARY 


A method has been described for analyzing 
multicomponent mixtures of hydrocarbons in the 
liquid phase by means of their infra-red absorp- 
tion spectra. A constant thickness absorption 
cell is used, and analytical curves are established 
on the basis of synthetic binary mixtures of the 
different components in the multicomponent 
mixtures to be analyzed. Two procedures for 
using the analytical curves to convert optical 
densities to concentrations by means of suc- 
cessive approximations have been described. 
Four- and five-component mixtures are com- 
pletely analyzed in less than, 45 minutes per 
sample and the average error is 0.5 percent of 
the total sample. 
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Wide Range Static Strength Testing Apparatus for Glass Rods* 


T. C. Bakert AND F. W. PRESTON 
Preston Laboratories, Butler, Pennsylvania 
(Received November 1, 1945) 


By adapting an electro-dynamic loudspeaker, an apparatus was assembled to test the 
strength of glass rods in flexure. This apparatus will apply the stress in 2 to 3 milliseconds 
and permit static loading for times as short as 0.01 second. Using an additional long time 
loading apparatus, which simultaneously tested eighteen rods, it was possible to run static 
fatigue tests on glass rods from 1000 seconds to 86,400 seconds. The method of testing and 


errors due to the apparatus are discussed. 





INTRODUCTION 

HENEVER an attempt has been made to 

measure the strength of glass, it has been 
found that the rate of application of the stress 
greatly affects the results. If the load applied to 
the specimens is increased continuously, it is 
found that the specimens that are loaded slowly 
break at much lower stresses than those which 
are loaded rapidly. If the load is applied by incre- 
ments, with definite time intervals between each 
addition, the specimens which are loaded with 
long waits between steps break at much lower 
stresses than those loaded with short waits be- 
tween the steps. Moreover, a specimen under 
constant load may survive several hours before it 
breaks. This apparent weakening of glass with 
increased time of loading has herein been denomi- 
nated ‘‘fatigue.’’ Since engineers commonly use 
the term “‘fatigue’’ to mean failure under oscil- 
latory or repeated loads, we may denote the 
fatigue dealt with here as “‘static’’ fatigue, and 
the other kind as ‘““dynamic’”’ fatigue. 

Naturally it is important to know, for both 
theoretical and commercial reasons, what the 
real nature of fatigue is. The theoretical questions 
are: what is the strength of glass for very long 
loading periods, and does its instantaneous 
strength approach the value suggested by mo- 
lecular theory? The commercial questions are: 


what safety factors should be allowed for bottles 


containing pressure, and what impact stresses 
will a glass object stand? Although there had 
been considerable speculation on these questions, 


- the data were so limited as to make no generaliza- 


tion really sound. The object of this work was the 


* This apparatus was developed at the Preston Labora- 
tories in 1938 and 1939 for a project sponsored by the 
Hartford-Empire Company. 

t Present address: Libbey-Owens-Ford Glass Company, 
Plaskon Division, Toledo, Ohio. 
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extension of the data by several orders of 
magnitude to both shorter and longer times. 

Probably the most relevant data would show 
how long a specimen survives at a given stress. 
Unfortunately, to answer this question directly 
with experimental methods would be too tedious, 
but the inverse question, whether or not the 
specimen survives a given length of time at a 
predetermined stress, lends itself readily to ex- 
perimental technique. After sufficient data had 
been collected to answer this question, plotted 
curves would provide the answer to the original 
question. Accordingly, the method finally decided 
upon was to load the test specimen for a definite 
interval of time at a predetermined stress. If the 
specimen survived the test, it would be loaded at 
an appreciably higher stress and the test repeated. 

Round glass rods were chosen as the only satis- 
factory test specimens because of their availa- 
bility in several compositions, and because they 
have no cut or ground surfaces toconfuse theissue. 

When a constant load is to endure for only 0.01 
second, the problem of applying the load and 
removing it is far more difficult than when the 
load is to be held constant for hours, because for 
the brief application the load must be built up 
with great speed and without producing oscilla- 
tion. Therefore, two types of apparatus were 
designed: a short time loading apparatus that 
would apply the load for periods ranging from 
0.01 second to 100 seconds, and a long time 
loading apparatus that would apply the load for 
periods ranging from 1000 seconds to 24 hours 
(86,400 seconds). 


DESCRIPTION OF APPARATUS 
The Short Time Loading Apparatus 


The short time loading apparatus is essentially 
a dynamic motor consisting of a modified radio 
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loudspeaker. Instead of the usual voice coil, an 
armature or moving coil assembly was developed 
which would withstand the loads that had to be 
applied. 

A schematic diagram of the apparatus is shown 
in Fig. 1. A glass rod (1) is supported by two 
knife edges (2). A wire bail (3) is looped over the 
middle of the glass rod and is connected to two 
lugs (4) on the armature (5). An oversize field 
coil (6) which is actuated by a 12-volt storage 
battery (7) produces the magneto-motive force 
of a central field magnet (8). The magnetic 
circuit is carried from the central magnet through 
a shell (9) to the air gap (10), in which the 
armature (5) is freely suspended from the rod. 
When the key (11) is closed, the armature is 
actuated by the current flowing from the battery 
(12) of about 120 volts. The current is regulated 
by a resistance (13) and is measured by an 




















Fic. 1. Schematic diagram of the quick-loading apparatus. 


ammeter (14). Forces up to 30 lb. can be pro- 
duced momentarily by the armature. 


Details of the Dynamic Motor 


An eighteen-inch dynamic loudspeaker was 
used as the basis of the dynamic motor. The size 
of the field coil was practically doubled to provide 
a much stronger magneto-motive force. The 
magnetic air gap was operated at approximately 
10,000 gauss for the smaller armatures, and as 
low as 7,000 gauss for the larger ones. In Fig. 2 
the rod is shown in position to be tested. Several 
armatures were designed for different purposes. 
For the 0.01-, 0.1-, and 1-sec. loading durations 
the requirements were light weight and large 
damping coefficient to permit quick action with- 
out overtravel. The time interval was so short 
that the motor coil could carry large currents 
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Fic. 2. A rod in position to be tested. 


without overheating. However, for the 10-second 
and 100-second periods the chief requirement 
was that the coil should not overheat, necessi- 
tating a larger motor coil and greater weight, yet 
the time of application of the load could be much 
longer so that smaller damping coefficients and 
greater weight could be tolerated. The load could 
be satisfactorily applied by reducing the circuit 
resistance by hand, thus bringing the current to 
its correct value in approximately } second. 
Figure 3 shows a cross-sectional diagram of one 
side of the armature. An aluminum form serves 
as a frame and a spool. The wire bail to the rod is 
fastened at the hole (1) in the lug (2). A stiffening 
ring (3) prevents the form from collapsing due to 
the pull of the wire bail. A motor winding (4), 
which carries the current from the 120-volt 
battery, consists of a double layer of No. 28 
wire. The damping winding (5) is made of several 
turns of fine, uninsulated, aluminum wire. How- 
ever, the damping winding is insulated from the 
aluminum spool and from the motor winding by 








Fic. 3. Cross section of voice coil. 
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Fic. 4. Trip scale for calibrating the dynamic motor. 


a layer of insulating enamel. It was possible to 
thicken the spool and omit the damping winding 
in the design of some of the armatures. 

The operation of the dynamic motor depends 
upon the fact that a force is exerted on a wire 
carrying current in a magnetic field. This force is 
given by F=BII, where F is the force, B is the 
flux density in gauss, / is the length of conductor 
in the field, and J is the current flowing through 
the wire, in c.g.s. electromagnetic units. Thus the 
force is directly proportional to the current in the 
motor winding. As long as there are approxi- 
mately the same number of turns of wire in the 
magnetic field, the force will be substantially 
constant for small displacements of the armature. 


._In practice, the armature can be displaced over a 


range of 3°; inch while the force remains constant 
to +4 percent. The quick-loading armature 
designed for 0.01-second loading produces a force 


of 3 pounds per ampere. 


When the motor winding moves across the 
magnetic field, a back e.m.f. is developed. This 
back e.m.f. is given by E=Blv, where E is the 
electromotive force in c.g.s. electromagnetic 
units, / is the length of conductor in the magnetic 
field, B is the flux density, and v is the velocity. 
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Thus, when the coil moves, it is possible to de- 
termine its velocity by means of an oscillograph 
and so tell whether it is behaving as it should. 

The damping of the armature is due to the 
counter currents induced in the ring and damping 
winding by the motion of these conductors in the 
magnetic field. By definition, the coefficient of 
damping is k=F/v where k is the damping 
constant, v is the velocity, and F is the force pro- 
duced in opposition to the velocity. By combining 
this with the formulas F=B//, and E=Bly, it 
turns out that k= B°?]/E=B*P?/R. Here R is the 
electrical resistance of the conductor in the 
circuit. In the present case the damping constant 
k depends on three factors: the damping due to 
the aluminum spool, that due to the motor 
winding, and that due to the damping winding. 
Each of these factors can be computed from the 
formula k = B*/?/R. The coefficient of damping in 
the spool is the most difficult to determine be- 
cause the path of the current extends beyond the 
boundaries of the uniform magnetic field. The 
currents in the damping winding and in the 
motor winding, on the other hand, are confined 
within the part of the field where the geometry is 
relatively simple. The three separate coefficients 
have to be added together to obtain the total 
damping. When an armature is made, a good 
estimate of the actual value of the coefficient can 
be obtained by studying the voltages across the 
motor coil with the aid of an oscillograph, as will 
be explained later. Corrections can be made by 
adding to or subtracting from the damping 
winding enough wire to give the desired coefficient. 

The contribution of the motor winding to the 
damping constant is rather small if the circuit 
resistance stays fairly large. With a field strength 
of 10,000 gauss, and with 50 ohms in the circuit, 
the contribution is 4 lb./sec., and with 20 ohms 
in the circuit the contribution is 10 lb./sec. The 
circuit resistance was usually at least 20 ohms. 
This variation of 6 lb./sec. will change the total 
damping coefficient by only 3 or 4 percent. 

A preload is applied to the armature by passing 
a small current through the motor coil in order to 
take the slack out of the wire bail. This preload 
puts a maximum stress upon the specimen of 
around 1000 to 2000 Ib./sq. in. This is sufficiently 
below the breaking stress at any of the durations, 
and under any of the conditions where it is used, 
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that it causes no appreciable fatigue. This matter 
will be discussed more fully later. 

Several armatures were designed for quick 
loading, because none is effective for more than a 
small range of rod sizes. For example, to test the 
rods that vary in diameter from 0.210 inch to 
0.230 inch, a coil assembly was designed which 
has a damping coefficient of 168 Ib./sec. and 
which weighs 23 grams. 

For the 10- and 100-second loads a six-layer 
coil was designed which provides adequate loads 
at current values low enough to prevent over- 
heating. The additional weight and reduced 
damping of this coil offer no problem at the 
relatively slow loading rates appropriate for these 
long loading periods. 


Calibrating the Dynamic Motor 


In order to measure the force exerted by the 
armature, a trip scale is mounted with the center 
of one pan over the motor coil. Figure 4 shows the 
trip scale in position. The wire bail of the arma- 
ture is attached below the empty pan, and the 
opposite pan is loaded with a ten-pound weight. 
The empty pan is held down by hand, keeping 
the armature in its working range, while the 
current is increased until the coil just supports 
the ten-pound weight. Such determinations can 
be made to about } percent accuracy with an 
ammeter readable to about } percent of full scale. 
A better ammeter would be of no value as there is 
sufficient damping and friction in the dynamic 
motor to make finer readings uncertain. With a 
satisfactory ammeter only one calibration point 
is necessary, since the relation between the cur- 
rent through the motor coil and the force exerted 
by it is as linear as can be measured, if the field 
magnet is properly saturated. Of course, the 
current through the field coil must be kept 
constant or the calibration will change. 


The Timing Device 


The timing device which was designed to time 
the 0.01-second, 0.1-second, and 1.0-second inter- 
vals appears in Fig. 5. It consists of a 100-lb. 
pendulum mounted on a heavy frame and has a 
period of 2 seconds. The pendulum is released by 
a magnet and is caught at the end of its swing 
by a latch. At the start of the swing a limit switch 
is closed to energize the field magnet; at the end 
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Fic. 5. The timing device. 


of the swing another limit switch breaks the field 
coil circuit when the test cycle is complete. 

The 0.01i-second and 0.1-second intervals are 
controlled by means of a cam on the pendulum 
bob, which closes a pair of contacts in the middle 
of the swing. At the end of the interval another 
cam opens the contacts, and the resulting arc is 
blown out with a powerful jet of CO». The arc 
is greatly diminished in about 0.0001 second, but 
it takes about 0.001 second to blow it out. 
Without the blast of CO, it is nearly impossible 
to put out the arc at all because the contacts open 
only } inch and currents up to 6 amperes at 120 
volts d.c. have to be interrupted. The timing 
cams are calibrated by observing the oscillograph 
pattern and making the indicated adjustments. 

The 1.0-second interval is timed by having the 
pendulum close a switch in the middle of its 
forward swing, and open it on the return swing. 
Since the pendulum has a 2-second period, this 
interval is one second. 

The 10- and 100-second intervals are timed by 
a stopwatch. A slidewire rheostat shown in Fig. 2 
and indicated by (13) in Fig. 1 is operated by 
hand to apply and release the load. 
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Fic. 6. The electrical control circuit. 


Accessories to the Apparatus 


Before many tests had been made, it became 
evident that a great deal of accessory apparatus 
was necessary. While the current in the loading 
circuit is being adjusted preparatory to a test, a 
dummy resistance is substituted for the motor 
coil in order to protect the latter from over- 
heating. This resistance has to be switched in 
when the current is adjusted and out when a rod 
is tested. The control resistance becomes heated 
during the loading current adjustment, so that an 
a.c. current must be supplied to it between times 
to keep it warm and maintain its resistance con- 
stant. The CO, used to blow out the arcs of the 
timing contacts is turned on and off with an 
electric valve, because otherwise the flow would 
quickly dissipate the supply. 

' It was soon found that the manual switching 
arrangement was too complicated. Turning the 
wrong switch could easily cause an accident or 
spoil the test. For example, the pendulum bob 


~had to be released when testing but not when 


calibrating. An interlocking circuit was finally 
evolved which was relatively foolproof and which 
reduced the switching procedure to a minimum. 
This circuit is shown in Fig. 6, and is described 
below. 


166 





It has been found that a cathode-ray oscillo- 
graph is indispensable. The timing cams are 
adjusted with the oscillograph, and the trace of 
every test is watched on the oscillograph screen 
to be sure that everything goes correctly. By 
watching the voltage developed across the motor 
winding as indicated by the oscillograph it is 
possible to tell whether the load goes on smoothly, 
whether the armature sticks, whether the rod 
breaks, or whether the contacts arc excessively. 


Operation of the Apparatus 


Because of the fact that the rods to be tested 
vary in diameter, the experimental procedure is 
somewhat complicated. First it is necessary to 
construct a chart giving the relationship between 
the maximum fiber stress of the rod and the cur- 
rent through the motor winding for every rod 
diameter. Each rod is carefully measured for 
maximum and minimum diameter and then is 
oriented in the apparatus so that the maximum 
fiber stress will come at one end of the mini- 
mum diameter. 

Referring to Fig. 6 the sequence of operations 
is as follows: the knife switch (S5) is closed to 
energize the field magnet, and the manual switch 
(S3) is set-on “‘calibrate’’; the start button ($1) 
is pressed, and the: current indicated by the 
ammeter (A 2) is set to the proper value (as given 
by the stress vs. current table) by adjusting the 
control resistance (R1). Finally, to test the 
specimen, the switch ($3) is set on “test,”” and 
the start button (§1) is pushed again. At this 
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Fic. 7. Theoretical curve of load of 0.01-second duration 
for rods 0.230 and 0.210 inch in diameter. 
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point it may be helpful to trace on Fig. 6 the 
sequence of part of the resulting operations. 

Closing (S1) energizes the coil of relay (1) 
which causes an open contact (C1) of this relay to 
close and energize the coil of relay (2). Other 
contacts (C1) open to turn off the a.c. supply 
maintaining the temperature of the control re- 
sistance (R1) and close to make the moving coil 
circuit ready for closure of the timing contacts 
(TC). The coil of relay (2) being energized by 
closure of contacts (C1), the contact points (C2) 
open in the pendulum release circuit, allowing the 
pendulum to fall. Thus the pendulum bob is 
released, and by closure of the open contacts (C1) 
of relay (1) the preload controlled by (R2) is 
turned on, the CO, valve is opened, and the 
timing contacts (TC) become operative. The 
pendulum bob swings down and operates the 
timing contacts (JC). This energizes the moving 
coil (MC) of the armature from the 120-volt 
battery, and loads the specimen. At the end of its 
swing, the pendulum bob operates a limit switch 
(S4) which shuts off the apparatus while the 
pendulum is caught by a mechanical latch. In the 
meantime the operator watches the oscillograph 
to see that everything has functioned properly. 
If the rod has not broken, the pendulum is reset, 
the current adjusted for a higher stress, and the 
whole procedure repeated. 


Dynamics of Quick Loading 


In order to get some idea of the accuracy of the 
measurements it is necessary to study the loading 
mechanism as a vibrating system. Since most of 
the mass is in the armature and most of the 
compliance is in the rod, the problem may be 
treated as a vibrating system with a single degree 
of freedom. The differential equation becomes: 


d*y dy 
m—+k—+sy=F, 
dt? dt 


where y is the displacement of the center of the 
rod, m is the effective mass of the whole system, k 
is the damping constant, s is the stiffness of the 
rod at its center, and F is the suddenly applied 
force. The load F is applied at t=0, so the equa- 
tion is subject to the boundary conditions: 


y=0 and dy/dt=0 when (t=0. 
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Fic. 8. Theoretical armature velocity curve during applica- 
tion of the load. 


The solution of the equation becomes: 


F 


yro--- 


Ss scos¢ 





e*t/2m cos (wt—¢), 


which is plotted for two rod diameters in Fig. 7. 
Then 
dy F »w 





ek tl2m sin wt, 
dt s cos*@ 


which is plotted in Fig. 8. In these equations: 
s ke 4 
(- a s k? \3 
and o-(-- ) : 
m 4m? 


cos ¢ =———_—_-- 
(s/m)} 

As shown in Fig. 7, the displacement is mini- 
mum at t=0 and passes through a maximum at 
2.8 milliseconds or, in general, at t=7/w. 

Likewise, as seen in Fig. 8, the velocity passes 
through a maximum followed by a minimum, and 
these occur at t=(2r/2—@)/w for the maximum 
and t= (32/2—¢)/w for the minimum. 

These correspond to the maximum and mini- 
mum e.m.f.’s across the motor coil as indicated 
at D and G in Fig. 9, which is described later. 

The ratio between the overtravel, i.e., the 
excess displacement over that for the steady load, 
and the steady load displacement is e~* **" *. The 
ratio of the minimum velocity to the maximum 
velocity is also e~* “" *, Therefore, because the 
voltage induced in the moving coil is proportional 
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to its velocity, e~* **" * must equal the ratio of the 
departures from the steady load voltage at D 
and G in Fig. 9. By measuring this ratio in the 
actual tests, the amount of overtravel and the 
damping coefficient can be readily computed. 
It turns out that for 2 percent overtravel, 
k =1.20(ms)'. For a range of rod sizes to be tested 
with an armature of given mass and damping 
coefficient, the largest diameter is determined by 
the condition that there should be no more than 
2 percent overtravel. The smallest diameter is 
determined by not allowing the duration of the 
load to be decreased by more than 10 percent for 
the 0.01-second test. 

Originally it was intended to load the rods in 
one millisecond, which would require the use of 
rods greater than } inch in diameter, with a 
4-inch span between supports. Unfortunately, 
the rather surprising increase of strength of glass 
in the shorter time ranges completely overtaxed 
the equipment with the larger rods. It would 
have required a larger field magnet, much higher 
d.c. voltages, and a general redesign of the appa- 
ratus. So we have had to use rods 0.220+0.010 
inch in diameter with a 5-inch span, and to be 


content with loading the rods in 2 to 3 milliseconds 


instead of the one millisecond originally planned. 
However, as one can see from Fig. 7, the time the 
rods are under stresses between 80 percent and 98 
percent of the full loads is not too great. 

A typical armature designed for measuring 


168. 
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Fic. 9. Oscillograph trace pro-« 
duced during a _ 0.01-second 
test. 


rods of approximately 0.220-inch diameter has 
the following constants: 


m=0.051 lb., k=168 lb./sec. 


It turns out that this coil causes a momentary 
2 percent overtravel on rods of 0.230-inch diame- 
ter, and it can be used with rods of 0.210-inch 
diameter without decreasing the time of full 
loading by more than 10 percent. Figure 7 shows 
a theoretical loading curve for rods 0.230 inch 
and 0.210 inch in diameter. The constants of the 
system with the 0.230-inch diameter rod are: 


m= .060 Ib., k&=168 lb./sec., 
and 
s = 194,000 Ib. /sec.?. 


The mass constant is increased over the mass 
of the armature by the effective mass of the rod. 
The constants for the rod 0.210 inch in diame- 
ter become: 
m=.059 lb.; k=168 Ib./sec., 
and 
s = 135,000 Ib. /sec.?. 


In Fig. 8 is a theoretical curve of the velocity 
of the armature for a rod 0.230 inch in diameter. 
It turns out that this pattern is closely followed 
when the velocity is studied in the cathode-ray 
oscillograph. 


Oscillograph Patterns 


Since the velocity of the armature produces an 
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e.m.f. given by 


E=Blv, and since v=dy/dt, 
then 


Fw ; 
E=Bl e*tl2m sin wt, 
$s cos? @ 


and we should expect the voltage developed 
across the armature to have this e.m.f. as one 
component. The other component is, of course, 
the JR drop across the coil. Accordingly, when 
the oscillograph is connected to read the voltage 
across the armature, a pattern is given as shown 
in Fig. 9. The line ABCDIJKL is a trace of the 
oscillograph in a 0.01-second test. AM is the 
IR drop due to the preload current; HP is the 
IR drop due to the full load current. The voltage 
DN is the maximum e.m.f. due to the velocity. 
The load comes on at about F, and the overtravel 
occurs in the region between F and G. The time 
the load is on is between F and J. If a rod breaks, 
a pattern results as shown by the dotted line STU. 

It is thus quite easy to determine from the 
oscillograph pattern whether the load has been 
properly applied, how long the load was on, 
whether the rod broke, and whether or not any 
deleterious accident, such as the armature 
sticking, took place during the test. It is some- 
what difficult to estimate the amount of overtravel 
by increasing the scale, because the 7R drop is 
considerably larger than the e.m.f. due to velocity 
in the 2 percent overtravel region. It is found 
easier to use an auxiliary winding of a few turns 
of No. 40 wire around the outside of the armature 
when damping coefficients are being estimated by 
the overtravel. 


Accuracy of Tests 


The quick-loading device is able to apply a 
load in 2 to 3 milliseconds, maintain it for 10 
milliseconds, and then remove it in } millisecond. 
Even when a range of rod sizes is tested, the error 
in timing the load is less than +5 percent, and 
the overtravel is less than 2 percent. Since the 
overload is on for only about 10 percent of the 
time it is doubtful that a 2 percent overload 
would cause as much as 1 percent difference in 
the results, because of the rapid increase of 
strength of glass with shortening of duration of 
loading. An added factor in favor of accuracy is 
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Fic. 10. Cross section of the long time loading device. 


that, when a large number of rods are tested, 
several will not be overloaded, because they are 
on the low side of the diameter range. Therefore, 
for the 0.01-second test, the average result can 
be relied upon to at least 1 percent as far as the 
overloading is concerned. 

The overload will not be serious in the 0.1- 
second and 1-second tests, because it will occur 
for only about 1 percent and 0.1 percent of the 
total time under load. The manner of application 
of the 10- and 100-second loads makes over- 
loading in these tests impossible. Thus it is safe to 
say that 1 percent accuracy for the load and 
+5 percent accuracy for the timing of the load 
are available for all the short-time loading points. 


The Long Time Loading Device 


The problem in the long time loading device 
was to test several rods simultaneously so that 
the length of time and expense of the tests would 
be kept to a minimum. Figure 10 shows a cross- 
sectional diagram of the apparatus. Two knife 
edges (1) are fastened to stiff channels (2), and 
support a rod (3). A central knife edge (4) is 
attached to a container (5) holding sand. When 
the rod breaks, the container falls on the angles 
(7) which are supported on a different frame from 
the channels. In this way, when a rod breaks, the 
fall of the container does not create a shock that 
would seriously affect the other rods still on test. 
The apparatus is in three tiers, each simultane- 
ously testing six rods. It is thus possible to 
conduct tests on 18 rods at a time. The apparatus 
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is located in the basement where the vibration 
conditions are quite satisfactory, as no heavy 
machinery or railroads are in the vicinity. Tests 
with dial indicators have shown no serious 
vibration present. 

A table was constructed giving the load versus 
the extreme fiber stress for every rod diameter. In 
making a test, the diameter of the rod is meas- 
ured, and the load is determined from the table. 
The container is then filled with sand to make up 
the proper load and a special levering device is 
used to apply the load gently. In general, an 





accuracy of 1 percent could be attained in 
measuring the load and its duration. The 1,000-, 
10,000-, and 86,400-second tests were made with 
the long-time loading device. 


ACKNOWLEDGMENT 


We are greatly indebted to Major J. M. 
McCormick, formerly of the Preston Labora- 
tories, for his work and helpful suggestions in 
constructing the switching system whereby the 
labor in making tests on the short time loading 
apparatus was greatly reduced. 





Fatigue of Glass under Static Loads 
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Static loading tests were run on glass and porcelain rods #3 inch in diameter for times 
ranging in duration from 0.01 second to 24 hours, by using specially designed apparatus. It was 
found that glass can support for 0.01 second about 3 times the stress that would break it in 
24 hours. The effects are generally the same for all glassy materials. Porcelain showed the 
effect somewhat less than glass. It appeared that adsorbed moisture and gases reduced the 


strength. 


INTRODUCTION 


T has been recognized for many years that 
glass fatigues* under a perfectly steady load, 
breaking unexpectedly and without apparent 
cause under a stress it has long endured. There- 
fore, in measuring the mechanical strength of 
glass, the duration of the test must be specified as 
well as various other circumstances, some of 
which are discussed in this paper. 

Whereas most glassware gets broken by impact, 
that is, by loads whose duration is measured in 
microseconds and which vary widely even in so 
short a time, yet some glassware is subjected to 
constant loads of long duration—for instance, 
soda, champagne, and beer bottles, aquarium 
windows, gauge glasses, and so on. An under- 
standing of the relation between short time 
strength and long time strength is therefore 
important. 

* This use of the term “fatigue’’ appears to us perfectly 
legitimate, m3 it is not the same as the usual engi- 
neering usage. If necessary, our fatigue may be distin- 
guished as ‘‘static’’ fatigue, and the engineering usage as 
“dynamic” or “‘vibration’’ fatigue. 


t Present address: Libbey-Owens-F ord JGlass’ Company, 
Plaskon Division, Toledo, Ohio. 
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This relationship has already been worked out 
roughly over a substantial range of loading 
intervals by several experimenters, so that from a 
strictly practical point of view there is reason to 
believe we know enough to risk no hazards to the 
public. The interest, accordingly, shifts to the 
nature of ‘‘static”’ fatigue, and to such questions 
as whether the glass has any strength at all under 
an exceedingly long loading. If the strength after 
a great period is not zero, what is the value of it, 
and what is the meaning of it? On the other hand, 
we have reason to suspect that at loads of very 
short duration the strength may be very great. Is 
it theoretically infinite, or is it finite? If the 
latter, how great is it, and can any clear meaning 
be ascribed to it? Or is the question somewhat 
meaningless, in view of the inertia effects which 
make “steady” loads of very short duration 
inconceivable? 

At first we had hoped our investigations would 
shed some light on these theoretical questions, 
but as more evidence accumulated all our beliefs 
as to the nature of fatigue in glass tended to 
change, and we felt much less inclined to specu- 
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late on the strength of glass and the ultimate 
effects of fatigue until the concept of strength 
was more precisely defined and the nature of 
fatigue was more thoroughly understood. Our 
initial feeling was that fatigue was due to the 
thermal effects at the tip of a flaw. The highly 
stressed bonds at that point were supposed to 
break when their momentary energies grew 
greater, due to the random shuffling of the 
thermal energies between the atoms. If this were 
the case, the cause of fatigue would be the same 
for long and short durations of load, and so we 
could expect a correlation between time and 
breaking strength. One might then assume that, 
if we had data encompassing a few orders of 
magnitude, we could extrapolate to both longer 
and shorter times with reasonable safety. How- 
ever, as our work progressed it soon became ap- 
parent that atmospheric conditions were a factor 
in the fatigue of glass and that the phenomenon 
was so complicated that we could have little 
justification for extrapolation to either shorter or 
longer times, except as mere speculation. 

The data presented are therefore disap- 
pointingly restricted, from the viewpoint of de- 
termining ultimate ‘“‘strength” of glass. On the 
other hand, certain practical questions are 
answered as to the reliability of testing pro- 
cedures, and at the same time a sounder under- 
standing of the nature of the phenomenon has 
emerged, which may lead to a valuable line of 
research. 

These tests were designed to investigate the 
phenomenon of strength decay, or fatigue, of 
glass in its normal condition, that is, with the 
amount of surface damage ordinarily present in 
commercial glasses, rather than the fatigue in 
glass having specially prepared, pristine surfaces. 
Glass samples were used, therefore, which had 
light surface damage easily visible under a micro- 
scope, which made their strengths comparable to 
that of commercial glassware. It is quite possible 
that all the breaks started at these flaws, and so 
the flaws in glass may be an essential part of the 
phenomena that we have studied. Other effects 
noticed in glass fibers, or other glass objects 
having pristine surfaces, may bear no simple 
relationship to the results obtained in the pres- 
ence of large flaws. 

Since it is not possible to experiment for an 
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indefinitely long time and report the results to 
the present generation, and since there are, on 
the other hand, practical difficulties in experi- 
ments for an infinitesimal time, we are limited to 
a range of duration which is not extremely 
extensive in ‘powers of ten’ (orders of magni- 
tude). From this relatively restricted range, we 
may perhaps hazard good guesses as to what 
might happen at somewhat shorter and some- 
what longer periods of loading. The validity of 
our guesses depends on how accurately we can 
plot a curve of strength versus duration of load 
over the accessible range. 

In the work herein reported, this range has 
been considerably extended beyond that of previ- 
ous investigators. The extension is largely at the 
short duration end. This direction was chosen by 
one of us (T.C.B.) as offering more immediate 
promise than the other, because experimental 
conditions can be more readily kept constant for 
a brief time than for a long time. For instance, 
conditions such as temperature and humidity are 
important. Further, a lot of experimenting can be 
done in a matter of weeks or months at the short 
time end, while only a little can be done in many 
years at the other. 

The essence of the short time loading experi- 
ments is the use of “electromagnetic loads,” 
which can be very powerful while substantially 
free of inertia, so that the loads can be applied in 
0.003 second without significant overtravel or 
oscillation. The quick-loading apparatus was 
described in a previous paper! in which its details 
and accuracy were discussed. This apparatus was 
not used for experiments lasting minutes or hours, 
since gravity provides a more convenient steady 
force. For these longer periods of testing, the 
essence of the set-up is the simultaneous testing 
of many specimens under conditions as constant 
as possible in respect to temperature and 
humidity. 

The tests were made in flexure on 35-inch glass 
rods, the simplest method experimentally. The 
shortest practical duration was about 0.01 second 
with our apparatus. Theoretically, somewhat 
shorter times are practical in making tests in 
direct tension on fibers, but the experimental 


1T. C. Baker and F. W. Preston, ‘‘Wide range static 
strength testing apparatus for glass rods,”’ J. App. Phys. 17, 
162 (1946). 
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difficulties were believed to be greater and apt to 
introduce more error. The longest tests lasted 24 
hours. It is clear from the shape of the curves that 
little more would have been learned by extending 
the time to ten days, or even to one hundred or 
one thousand days. Furthermore, since each 
specimen must be tested at several successively 
greater loads, such tests would have taken many 
years tocomplete, even if experimental conditions 
could have been maintained constant. 

It is not unthinkable that the curve may yet be 
extended, experimentally, farther than we have 
gone, but it seems at present unlikely that much 
practical, or even theoretical, information could 
be acquired thereby. Therefore, in the belief that 
the subject has reached a measure of complete- 
ness, we present an account of some of our 
experiments. 


EXPERIMENTAL PROCEDURE 


Because the strength of glass seems to be a 
function of its previous history of exposure both 
to the atmosphere and to mechanical contact, 
special precautions had to be taken in order to 
keep constant the effects of handling, adsorbed 
moisture, etc., in measuring the variation of 
breaking stress with changing duration of the 
load. Ideally, the controlled experimental pro- 
cedure should start during the manufacture of 
the glass canes from which the rods were cut. 
Unfortunately, that was impossible. Our experi- 
mental procedure was constrained to begin after 
the glass reached the laboratory, which naturally 
left uncertain the previous history of the canes. 

The initial step was to select canes that met 
certain specifications as to diameter and taper, 
and cut them into rods six inches long. These rods 
were then distributed among six or eight samples, 
of 25 rods each, and their diameters were meas- 
ured and recorded. The rods were next set aside 
for a week or two so that any recent mechanical 
damage would tend to reach equilibrium, and the 
rods would have a chance also to come to 
equilibrium with the laboratory atmosphere. 
-Finally, the rods were broken in the strength- 
testing apparatus. 


Selection of Canes 


In order to get comparable test samples, glass 
canes were used, whenever possible, which had 
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been made at the same time. After being cut into 
four-foot lengths the canes were sent to the 
Hartford-Empire Company-to be annealed, where 
special precautions were taken to prevent extra 
damage during the annealing process. The 
annealed canes were then individually wrapped 
in tissue paper and sent to the Preston Labora- 
tories. There each cane for a test lot was ex- 
amined, and only those which met the following 
specifications were used. 


Specifications for Cane 


Taper per foot: -less than 0.005 inch 
Out of roundness: less than 0.005 inch 
Diameter: 0.220+0.010 inch 


Radius of curvature of axis: greater than 5 ft. 
Canes showing unusual abrasion were rejected. 


Each cane was cut into six-inch lengths, or 
rods, and these were distributed among either six 
or eight samples, depending on whether the cane 
was three feet or four feet long. The method of 
distribution was as follows: the first rod from the 
first cane was put into the first sample, the 
second rod was put into the second sample, and 
soon. The second cane was distributed by putting 
the first rod into the second sample, the second 
rod into the third sample, and so on, until the 
last rod was put into the first sample. The subse- 
quent canes followed out the progression. This 
method of distribution meant that each sample 
contained a rod from each cane, and each con- 
tained approximately the same number of rods 
cut from the ends and the middle of the canes. 
Thus 25 four-foot canes would make eight samples 
of 25 rods each, 25 three-foot canes would make 
six similar samples. 


Care of Test Samples 

As the canes were cut, the six-inch rods were 
put into wooden blocks having twenty-six num- 
bered holes bored 1 inch deep. Only twenty of 
these specimens were used for the tests, the other 
rods being saved for spares. For example, if an 
accident occurred in the test on one rod, the rods 
taken from the same cane were discarded from 
all the other samples, and a complete set of 
spares was substituted. Great care was used in 
handling the rods during and after cutting, be- 
cause new damage inflicted during subsequent 
handling could destroy the homogeneity of the 
specimens. 
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The ends of the rods were measured for maxi- 
mum and minimum diameter, and a scratch was 
made at one end to indicate the direction of the 
minimum diameter of the rod. The central 
minimum diameter of the rod was estimated to 
be the average of the minimum diameters at both 
ends, and was recorded on the data sheet, next to 
the specimen number. In testing each individual 
rod the apparatus was adjusted to produce the 
appropriate stresses in accordance with special 
tables showing, for each value of load, the appa- 
ratus adjustment in terms of rod diameter. 


The Apparatus 


The rods were tested in cross bending, with a 
five-inch span. The load was applied by a single 
central knife edge. The two types of apparatus 
for testing the rods have been described in the 
previous paper. Briefly, one consisted of a 
modified dynamic loudspeaker to apply the loads 
of short duration (0.01 sec. to 100 sec.), and the 
other was a multiple loading device designed to 
accommodate numerous rods for the 1,000-sec. to 
24-hour tests. 

The load was applied to the rods as quickly as 
possible, left on for a predetermined time, and 
then removed as quickly as possible. The load 
was applied in two or three milliseconds and re- 
moved in 4 millisecond in the 0.01-, 0.10-, and 
1.0-second tests. In the 10- and 100-second tests 
about } second was required to apply and to 
remove the load, and in the tests of 1000 seconds 
and over about four or five seconds were required. 
For the short time loading tests, tables were 
constructed giving the maximum stress in a rod 
as a function of the current through the armature 
of the dynamic loudspeaker, for all values of rod 
diameter. In making the tests, each rod was 
removed from the sample holder and placed in 
the testing apparatus with the minimum diame- 
ter in the direction of the load. This assured that 
the maximum fiber stress would be applied at the 
end of the minimum diameter. The rod was then 
tested at the minimum breaking load or at one 
loading interval below that. If the rod did not 
break, it was tested at an appreciably higher 
stress, the increment depending on the average 
strength. 

In the case of the long time loading apparatus 
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the procedure was very similar, except that 
instead of the moving coil used in the electro- 
magnetic loading a two-gallon can filled with 
sand was used to supply the load. 


The Loading Increment 


Glass ‘“‘fatigues’’ under constant load, and 
unless it recovers promptly upon the removal of 
that load (a matter that is not certain), the 
breaking of the glass at a given load might be 
affected adversely by any previous loadings, even 
though these might be of less intensity. Experi- 
mental evidence had been accumulated in the 
past, however, that when step-wise loading is 
used, previous, lower loads had very little or no 
influence if they were 10 percent or so below the 
breaking load. Thus if our increments of loading 
are 10 percent or more above previous loadings, 
we can ignore the effects of the latter. For good 
measure, it might be well to make the increments 
of the order of 20 percent. 

On the other hand, if we use large increments 
the statistical data will fall into very few class 
intervals and the average will be uncertain. 
However, it was determined that, when we had 
five class intervals, the error would probably be 
quite a bit less than that from other sources. 
Accordingly, the loading increment chosen was 
500 lb./sq. in. for samples whose averages were 
between 3000 and 6000 Ib./sq. in.; 1000 Ib./sq. in. 
for samples whose averages were between 6000 
and 12,000 lb./sq. in.; and 2,000 Ib./sq. in. for 
samples whose averages were over 12,000 Ib./sq. 
in. The strength was generally found to vary over 
a 2 to 1 range for any sample. Thus in a sample 
with 12,000 Ib. /sq. in. average strength some 
specimens would break at about 8000 Ib./sq. in. 
and some at about 16,000 Ib./sq. in. This resulted 
in five class intervals for the data, using steps of 
2000 Ib./sq. in. 

If arod survived a test at 14,000 Ib./sq. in. and 
broke at 16,000 Ib./sq. in. it would be recorded as 
having broken at 16,000 Ib./sq. in. Obviously, it 
could have broken at a load anywhere between 
14,000 and 16,000 lb./sq. in., and on the average 
would break at 15,000 Ib./sq. in. Therefore, it 
was necessary to apply the usual statistical cor- 
rection for the magnitude of the class interval. 
When the loading increment was 2000 Ib./sq. in. 
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the correction would be — 1000 Ib./sq. in. When 
the loading increment was 1,000 Ib./sq. in., then 
the correction would be —500 Ib./sq. in. 


Special Precautions 


Special precautions were necessary to keep the 
temperature and atmospheric conditions as con- 
stant as possible during the tests. The room 
temperature was kept as near 75°F as possible. 
The rods were placed near the apparatus previous 
to the test so that thermal strains would not be 
introduced when the rods were inserted in the 
apparatus. A few minutes before the test the rod 
was wet with water at room temperature, and 
water was applied to the rod whenever it became 
dry during the successive loadings. 





RESULTS OF VARIOUS STATIC LOADING TESTS 


The results herein reported were obtained from 
tests run on five glasses and three porcelains. 
These are presented in Table I, which gives, for 
each loading period, the average breaking strength 
of 20 to 25 specimens. The results for the various 
glasses are plotted in Figs. 1 and 2, and the results 
for the porcelains are plotted in Figs. 3 and 4. 
All results are corrected for the loading incre- 
ments, as described in the preceding section. 

Lot I consisted of Pyrex rods which had already 
been cut to length, since there were no long canes 
available at the time of the tests. In that respect, 
this constituted an exception to our usual pro- 
cedure. The rods were rather strongly disan- 
nealed, with a surface compression of 2000 to 
4000 Ib./sq. in., the average of several rods being 
about 3000 Ib./sq. in. Moreover, as these were 
our first tests and the effect of moisture was not 
appreciated at that time, the Pyrex rods were 
tested dry. Various difficulties, including the fact 
that no really reliable average for the surface 
compression could be taken, greatly reduced the 
accuracy of these results. 

Lot II was made up of the same type of disan- 
nealed Pyrex rods as used in lot | except that 
they were scratched circumferentially with a 
medium grade of emery. By this time the effects 
of humidity had been detected and the rods were 
tested wet. 

Lot III consisted of samples made from four- 
foot, annealed, machine-drawn, soda-lime canes 
which were unusually straight and round and 
made an excellent stock for selection of samples. 


TABLE I. Average breaking stress in pounds per square inch for various durations of load. 

















Ratio Ratio 
of 0.1- of 0.01- 
to 24- to 24- 
Test Wet 0.01 0.1 1.0 10 100 1000 86,400 hour hour 
Lot Material or Dry sec. sec. : sec. sec. sec. sec. test test 
I Disannealed Pyrex Dry 24,600 19,600 17,100 15,900 12,800 11,300 1.7 2.2 
Il Disannealed Wet 12,900 10,600 9000 8250 5750 5050 2.1 2.5 
scratched Pyrex 
Ill Annealed soda- Wet 19,900 15,900 10,500 7950 6450 2.4 3.3 
lime glass 
IV Annealed Wet 19,000 14,900 11,200 9400 7950 2.4 
lead glass 
V Fused silica Wet 23,880 19,800 15,660 14,280 11,740 2.0 
VI Porcelain Type A Dry 15,800 13,900 12,300 9700 9100 = 11.53 1.74 
VII Porcelain Type B Dry 15,400 13,500 11,900 9900 9400 9100 8900 =1.52 1.73 
VIII Porcelain Type C Dry 10,900 9400 8400 6900 6000 5900 1.59 1.85 
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Their composition was very close to that of bottle 
glass. Considerable care was taken in annealing 
and preparing comparable samples, and when 
certain check tests were run, the average strengths 
agreed to better than one percent, indicating 
excellent uniformity of the samples. These rods 
were tested wet. 

Lot IV was composed of glass rods containing 
approximately 50 percent lead’ by weight. The 
glass was unusually soft and easy to scratch, and 
the resulting damage on some rods made the 
tests somewhat erratic. These rods were an- 
nealed, and the tests were conducted with the 
rods wet. 

Lot V consisted of rods made from silica-glass 
canes three feet long. As the rods were apparently 
free of stress, no special annealing treatment was 
given them. All rods were tested wet. 

Lot VI was composed of porcelain rods* of the 
following composition, fired at cone No. 9. 

22 percent Kentucky Ball Clay 
25 percent Florida Kaolin 
37 percent Canadian Feldspar 


10 percent Pennsylvania Flint 
6 percent Water 


The rods were tested dry because the effect of 
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Fic. 2. III. Annealed soda-lime glass, wet; IV. Annealed 
lead glass, wet; V. Fused silica, wet. 


moisture on the strength of porcelain is con- 
siderably less than on the strength of glass. The 
possible errors caused by not testing wet are 
probably less than two or three percent. 

* Porcelain rods used in these experiments were supplied 


by George Aderhold of the Saxonburg Potteries, Saxonburg 
Pennsylvania. 
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_ Lot VII had the same composition as lot VI 
but the porcelain was fired at cone No. 10 instead 
of cone No. 9. These rods were also tested dry. 
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Fic. 4. VII. Porcelain, type B, dry; VIII. Porcelain, 


type C, dry. 


Lot VIII consisted of porcelain rods of the 
following composition, fired at cone No. 12: 


50 percent Canadian Feldspar 
10 percent Florida Kaolin 

30 percent Kentucky Ball Clay 
10 percent Water 


The rods were tested dry. 

Of the results shown in Table I, those for soda- 
lime glass, lead glass, and fused silica may be 
considered most reliable. In the case of the Pyrex 
glass the records are ambiguous regarding the 
amount of the correction that was applied for the 
disannealing present in the rods. The amount of 
surface stress in the porcelain rods due to 
disannealing could not be determined, so some 
error from this cause may be present in these 
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data also. The rods of lots I1I-V, however, were 
known to be free of any disannealing stress and 
these data, therefore, do not include the possi- 
bility of such an error. 


DISCUSSION OF ERRORS 


At the present time, an exhaustive treatment 
of the errors in the fatigue testing of glass is quite 
impossible, because of the entanglement of the 
variables, which cannot as yet be independently 
investigated. However, a partial list of the 
factors will be given so that the reader may be 
assisted in evaluating the results. 


1. Errors in Apparatus 


In a previous paper the errors involved in 
applying the load and the errors of timing the 
duration of load were discussed. In general it was 
felt that these errors were not greater than 
+5 percent. 


2. Errors in Measuring the Rod Diameters 


This error was not large because the diameters 
could be measured to 0.0005 inch, or } percent. 
Such an error would lead to only ? percent error 
in the applied stress, and since it was a statistical 
error it would practically disappear in the aver- 
age for twenty rods. 


3. Errors Due to Single Knife Edge 


Since the rods were measured in cross bending, 
with a 5-inch span and single knife edge loading, 
a certain amount of error existed when the rod 
did not break exactly underneath the knife edge. 
At a quarter-inch away from the central loading 
point the stress would be 10 percent less than at 
the center. This changing stress caused a peculiar 
type of error to enter the results, because the 
long time loads nearly always broke the rod 
exactly under the knife edge, whereas under the 
short time loads the rod often broke as much as 
} inch away from the central point. This would 
lead to a correction upwards in the short time 
strength, because the portion under highest 
stress was not broken, and therefore would have 
broken only at a still higher stress. Scratched rods 
nearly always broke under the knife edge, even in 
short time tests. The effect seems also to be a 
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function of the number and kind of flaws, It is 
estimated that this error could reduce the aver- 
age by as much as 5 percent in the short time 
loads. 


4. Error Due to Scattering of the Data 


The strength of the strongest rod in a sample 
might be twice as great as the weakest one. This 
scattering of the data resulted in a standard 
deviation for the individual rods that was about 
10 percent of the average of the sample. The 
standard deviation of the average of 20 to 25 rods 
would be about 2 to 2.5 percent of the average. 
It would appear’ from the foregoing that the 
average strengths of different samples could be 
expected to agree to only 5 percent, and differ- 
ences of as much as 8 percent would not be 
infrequent. On the other hand, actual check tests 
indicated that the procedure of preparing the 
samples had led to a strong correlation between 
them and that generally the average strength 
from sample to sample agreed to one percent or 
better. Unfortunately, insufficiency of data pre- 
vents an accurate computation of the correlation 
coefficient; and there is no assurance that the 
averages of 20 to 25 rods might not differ by as 
much as 5 percent due to purely statistical 
fluctuations. 


5. Errors Due to the Loading Increments 


The choice of the loading increments involved 
a compromise between two errors. One error 
would be that caused by the weakening of the rod 
from previous loading, and would become greater 
as the increment was made smaller. The other 
error would be due to an insufficient number of 
statistical class intervals, and would increase as 
the increments were made larger. This latter 
error is fully described in books on statistics? so 
that it need not be discussed here, except to state 
that it would probably be less than one percent 
in the average, when a loading increment of be- 
tween 10 and 20 percent of the average is used. 

The weakening effect of the previous loading 
has not been fully investigated. Some tests were 
made which indicated that a loading increment 


2 See, for example: G. U. Yule and M. G. Kendall, An 
Introduction to the Theory of Statistics (Charles Griffin and 
Company, Ltd., London, 1940) Chap. 6. 
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of 10 percent of the average or more, would cause 
less than a one percent decrease in the average 
strength. 


6. Errors Due to the Atmospheric Conditions 


Considerable data were obtained indicating 
that moisture and possibly gases affect the 
strength of glass. The specimens were tested wet 
after the phenomenon was discovered, but even 
then the presence of gases such as CO: possibly 
affected the results from one to five percent, the 
greater effect being for the longer load durations. 


7. Errors Due to Temperature 


Evidence has been obtained showing that tem- 
perature has a profound effect on both the 
strength and the fatigue of glass. It turns out that 
even small changes such as occur indoors will 
cause a change in the strength of glass amounting 
to one or two percent. 


8. Errors Due to Storage of Rods 


Since it was impossible to make all the tests at 
the same time, some of the rods were stored 
longer than others while the atmospheric condi- 
tions changed. There was some evidence that 
when rods were stored in a dry atmosphere for 
some time they did not have the same strength as 
rods stored in a damp atmosphere, even though 
they were tested wet in both cases. This resulted 
in another possible error of one or two percent in 
the final averages. 


9. Combination of Errors 


Considering these errors as a whole, it is esti- 
mated that the results cannot be considered 
accurate and controlled to more than +6 percent. 
Moreover, it must be appreciated that certain 
types of error will tend to distort or bias the data, 
making the short time strengths appear less than 
they should be in comparison with the long time 
strengths, or vice versa. Therefore, any generaliza- 
tions that depend upon an accuracy of better 
than 10 percent in the results, or series of results, 
should be treated with extreme caution. 


DISCUSSION OF RESULTS 


Since a primary objective of this work was to 
investigate the fatigue in commercial glassware 
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under normal conditions of temperature and 
humidity, the data are limited to normally 
scratched glass tested either wet or in an atmos- 
phere of high humidity. As the difference in the 
strength of glass between 30 percent humidity 
and completely wet is only about 5 percent, the 
results are quite comparable. The surface con- 
dition of the rods was similar to the normally 
scratched and abraded condition found on bottles 
and window glass. Scratches, chatter marks, and 
percussion cones could be easily seen under a 
50-power microscope. 


Amount of Fatigue in Glassy Materials 


The fact most apparent from the data is that 
all of the glassy materials tested exhibit the 
phenomenon of static fatigue to a marked degree. 
In the case of the soda-lime rods the 0.01-second 
strength is 3.1 times the 24-hour strength. The 
ratio of the 0.1-second strength to the 24-hour 
strength is between 1.7 and 2.4 for all the glasses 
tested. The corresponding ratios for porcelain 
ranged from 1.52 to 1.59. 


Nature of Fatigue Curves 

The plotted curves giving the relationship 
between the breaking strength and the logarithm 
of the duration of test show that there is a rapid 
decrease in strength up to a minute or so, and a 
very slow decrease from then on. The relationship 
between the strength and the logarithm of the 
duration of load is not linear, as has been previ- 
ously supposed.* At present there is some ques- 
tion as to whether or not the curve is truly a 
monotonic decay, with steadily decreasing slope 
and steadily increasing radius of curvature, as we 
go to longer times. There remains a possibility, 
within the limit of accuracy of the tests, that 
there may be several points of inflection, implying 
deviations from our drawn curves. 


Effect of the Condition of the Surface 
upon Fatigue 
The condition of the surface has a profound 
effect on fatigue as well as upon the strength of 
the glass. Scratched glass loses much less strength 
in a given time interval than does unscratched 
glass. However, the percent loss of strength is 
roughly the same. This may be seen by comparing 


3A. Smekal, Ergebr. d. exakt. Naturwiss. 15, 106-188 
(1936). 
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the tests on the scratched and unscratched Pyrex 
rods. The decrease for the scratched rods was 
7850 lb./sq. in. and the decrease for the un- 
scratched rods was 13,300. Yet the ratios of the 
. 0.01-second strength to the 24-hour strength were 
2.5 and 2.2, respectively. Unfortunately, this is 
all the data that are available for the comparison 
of scratched with unscratched rods, and the 
variable annealing stresses that existed in these 
rods as well as the different moisture conditions, 
tended to obscure the results. 

That badly damaged glass fatigued in the same 
ratio as less damaged glass was borne out by the 
fact that the percent standard deviation of a 
sample of each kind tested at 0.01 second was ap- 
proximately the same as that of a sample tested 
for 24 hours, although the strengths had changed 
by a factor of 3.1 to 1. Thus in the case of the 
soda-lime glass the standard deviation for the 
0.01-second test, in which the average was 8500 
lb./sq. in., was 800 lb./sq. in. In both cases the 
Pearson coefficient of variation was approxi- 
mately 10 percent. Thus the stronger rods lost 
more strength than the weaker ones, yet they 
lost it in the same ratio. On the other hand, the 
fact that the rods broke under the central knife 
edge in long duration tests, but tended to break 
at some distance from it during the short dura- 
tion tests, indicates that strong glass (i.e., glass 
with little surface damage) has a greater per- 
centage loss of strength than weak glass. For it 
may be argued that in the short time tests some 
flaws made the rods weaker than other flaws did 
and the rods would break at 10 percent less stress. 
If these flaws were only } inch away from the 
central knife edge, the rod would break there. 
Then if the minor and major flaws reduced the 
strength in the same ratio, we should get the 
breaks in about the same way on the long time 
tests. The fact that this does not happen indi- 
cates that the minor flaws must cause a greater 
fatigue rate in order to originate the breakage on 
the long time tests. 


Effect of Chemical Composition 


Because the loss of strength of glass is con- 
nected with surface damage, which cannot be 
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accurately standardized between different glasses, 
it is difficult to make accurate judgments con- 
cerning the effect of chemical composition upon 
fatigue. Apparently the short time to long time 
strength ratios are about the same for all glasses 
studied. The results on the disannealed Pyrex 
rods have to be discounted because of the 
impossibility of applying a good correction for 
annealing stresses. We are left, therefore, with 
the results for silica glass, lead glass, and soda- 
lime glass as being perhaps accurate enough to 
judge this point. The ratio between the 0.10- 
second test and the 24-hour test is 2.4 in the last 
two cases, so no distinguishable difference is 
evident. The ratio however is only 2.0 in the case 
of the silica glass, so there is some evidence that 
silica glass fatigues slightly less than soda-lime or 
lead glasses. 

In the case of the porcelains there is no question 
that they fatigue much less than the glass, as the 
ratio of the 0.10-second test to the 24-hour test 
ranges from 1.52 to 1.59. However, so little is 
known about the nature of fracture in porcelains, 
or aggregates in general, that it would be futile 
at this time to try to explain the fatigue in 
porcelain in terms of the fatigue phenomenon in 
glass, although they are probably closely related. 


SUMMARY OF RESULTS 


The results have corroborated previous investi- 
gations showing that the strength of glass depends 
on its surface damage, on the time it is under 
load, and upon atmospheric conditions. As the 
work progressed it became evident that these 
factors were interdependent in a complicated 
way. Nevertheless, to an accuracy of about 10 
percent we can be assured that glass has about 
three times the strength under our briefest loads 
that it has under long time loads; that the 
relative loss of strength with time is about the 
same for lightly damaged and badly damaged 
glass surfaces, although the absolute strength 
varies radically with the damage; and finally, 
that chemical composition has only a small effect, 
if any, on fatigue. 
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The Effect of Water on the Strength of Glass 


T. C. BaKerf AND F. W. PRESTON 
Preston Laboratories, Butler, Pennsylvania 


(Received November 1, 1945) 


The strength of glass was measured in various surrounding media. It was found that glass 
is 20 percent stronger when dry than when wet, and two to two and a half times as strong 
when baked in vacuum than when tested wet (on tests of ten seconds duration). Apparently, 
the chief cause of the loss of strength, as compared with that in vacuum, is moisture, but gases, 
especially CO2, seem to have some effect. It was found that fatigue of glass disappears when 


the glass is tested in a vacuum. 


Autoclaving tests showed that glass is rapidly etched by water at higher temperatures, and 
that even silica glass is considerably attacked. It appears that the attack of water on glassy 
surfaces produces gels, and this attack goes on inside the flaws, thus considerably complicating 


the problem of the strength of glass. 





INTRODUCTION 


HE previous two papers!” have shown that, 

under the conditions of the experiments 
therein described, there is a precipitate fall in 
the breaking strength of glass, including fused 
silica, and of porcelains in the first few milli- 
seconds of loading, and a slower but steady fall 
with the lapse of further time. They have also 
shown that moisture seems to have something 
to do with this, so that the breakage of glass can 
perhaps not be understood simply in terms of 
stress, strain, and elastic constants; rather does 
“strength,’’ as we usually measure it, involve 
the properties of the surrounding medium, per- 
haps a reaction with the surrounding medium, 
and so pass out of the field of engineering and 
even of physics into the domain of physical 
chemistry. 

There have, indeed, both prior to the experi- 
ment detailed below and also between the time 
of the experiments (c1937-39) and the present 
(1945), appeared several papers by other investi- 
gators on the effect of the surrounding medium. 
Milligan,* for instance, found that exposing 
freshly cracked glass to water vapor or alcohol 
vapor considerably reduced its strength. He 
found that moistening a freshly scratched piece 


t Present address: Libbey-Owens-Ford Glass Company, 
Plaskon Division, Toledo, Ohio. 

1T. C. Baker and F. W. Preston, ‘Wide range static 
strength testing apparatus for glass rods,”’ J. App. Phys. 
17, 162 (1946). 

2 T. C. Baker and F. W. Preston, “Fatigue of glass under 
static loads,”’ J. App. Phys. 17, 170 (1946). 

°L. H. Milligan, “Strength of glass containing cracks,” 
J. Soc. Glass Tech. 13, 351 (1929). 


VOLUME 17, MARCH, 1946 


of glass reduced its strength 20 percent or more, 
and immersing the freshly scratched glass in dry 
paraffin oil increased the strength about 20 
percent. Initial experiments of the authors indi- 
cated that changes of the relative humidity of 
the atmosphere would affect the results of 
strength measurements by 5 percent or more. 
It was thus evident that the surrounding medium 
played a very important part in tests for the 
strength of glass. . 

So it appeared desirable to investigate the 
effects of the surrounding medium on the strength 
of glass in greater detail than had been done 
heretofore. It was hoped that the studies would 
shed light on the causes of fatigue and possibly 
lead to a more complete understanding of the 
phenomenon of breakage. 

Accordingly, the glass specimens were tested 
after they had been put into various solutions, 
and after they had been dried in desiccators; 
they were even tested in vacuum. Unfortunately, 
the results of the tests were generally not clear 
cut, and often it was not possible to duplicate 
them, so the general picture became more com- 
plicated with every new experiment. At the time 
the work was stopped a great many loose ends 
remained, and our understanding of the strength 
of glass had not become sufficiently fixed so 
that we could feel sure of our ground. Neverthe- 
less, from a study of our data, in combination 
with that available in the literature, a general 
picture of the strength of glass began to emerge, 
indicating a muéh greater complexity than is 
usually supposed. 
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First of all, commercial glass has tiny flaws in 
its. surface which can easily be seen with a micro- 
scope. After the flaw is made, water and gases 
enter it and radically change the nature of the 
surfaces and alter the physical and chemical 
conditions at the tip of the crack. If the glass at 
the flaw is under tensile stress, the water vapor 
affects the bonds at the tip of the flaw, causing 
great weakening. If the flaw is allowed to stand, 
the water vapor forms thin gels on the surfaces 
and rounds out the tip of the flaw, causing an 
increase in strength. There is reason to believe 
that complete equilibrium was not established 
in any of our experiments. 

The “strength” of glass began to appear to 
have little meaning as a physical constant of the 
material. It might be more correct to say that 
the phenomenon of rupture is due to several 
causes, including tensile stresses at the surface of 
the glass, the time during which the tension is 
applied, the nature and depth of the flaws, and 
the penetration of atmospheric or other reagents 
into the inside of the flaws. It is obvious that a 
great stress in itself would be a sufficient condi- 
tion to break the glass, yet there exists evidence 
that atmospheric and weather conditions also 
would cause some flaws to extend to the point of 
rupture.‘ So it is doubtful that any one single 
cause is an essential element in the breakage. 

With such a complicated picture to investigate 
it becomes apparent that controlled tests are 
particularly difficult. Many of the factors cannot 
be truly controlled, such as the size and orienta- 
tion of a flaw, so it is not surprising that effects 
should occur that could not be duplicated. For 
this reason-much of the data in this paper has 
been included not to prove points finally but 
merely to indicate the nature of effects, so that 
future research on this problem will be on a 
sounder basis. 

The strength of ordinary commercial glass 
depends primarily on the size, depth, and orienta- 
tion of the tiny cracks and flaws on its surface. 
These flaws cause the strength to be reduced 
from a theoretical strength of several million 
pounds per square inch to one that varies be- 
tween 2000 Ib./sq. in. and 20,000 Ib./sq. in. in 


* See A. Silverman and L. G. Ghering, ‘‘Deterioration of 
old whiskey bottles,’’ Glass Ind. [9] 15, 339 (1942). 
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glassware whose surface has been damaged by 
abrasion and impact. The effects of the surround- 
ing medium upon the strength, therefore, can 
only be comparative. So, in our experiments, 
two or more samples of about 25 glass rods each 
were tested under different conditions and the 
ratio of their strengths determined. 

Any surrounding medium seems to lower the 
strength by a certain ratio below the vacuum 
strength, rather than to subtract a fixed amount 
from the vacuum strength. For example, if the 
strengths of two comparable samples were ob- 
tained, one in vacuum and the other in air, and 
the strength in vacuum was 20,000 Ib./sq. in. 
and in dry air 10,000 Ib./sq. in., the ratio would 
be 2.0. Two other somewhat stronger samples 
tested in vacuum and dry air might have 
strengths of 30,000 and 15,000 Ib./sq. in., and in 
this case also the ratio would be 2.0. This rule, 
that the surrounding medium affects the strength 
by a ratio, seems to be approximately true of all 
our studies on fatigue of glass as affected by 
the surrounding medium. 


APPARATUS 


The apparatus used to measure strengths of 
glass rods in the following experiments was the 
quick-loading apparatus described in a previous 
paper.' When it became necessary to measure 
the strength of glass in vacuum, the quick- 
loading apparatus was put into a vacuum 
chamber. Otherwise the rods were tested out in 
the open air after they had been treated accord- 


' ing to the various procedures described. 


The vacuum strength-testing apparatus in- 
cluded a vacuum chamber, consisting of a bell 
jar resting on a cylindrical steel shell, in which 
the entire quick-loading apparatus. was placed. 
A glass rod is supported by two knife edges five 
inches apart. A wire bale is looped over the glass 
rod and connects to the armature which hangs 
freely in a magnetic air gap. A field coil produces 
the magnetic field in the air gap. This field coil 
was covered with cement to prevent outgassing 
when the vacuum was drawn. On either side of 
the glass rod were added two heater coils con- 
sisting of chromel wire wrapped around porcelain 
rods. An easily removable furnace shield pre- 
vented radiation from the heater coils from 
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warming the bell jar too much while the rod was 
being baked out. A glass window at the end of 
the furnace shield permitted observation of the 
rod during the test. The shield had several 
openings to the main chamber so that there was 
no doubt of its having the same vacuum as 
that in the main chamber. The vacuum chamber 
is enclosed by a large bell jar resting on a steel 
cylinder which has a flat bottom welded to it. 
The vacuum chamber is evacuated through an 
outlet leading to a McLeod gage and a diffusion 
pump. 
PROCEDURE 


The procedure for testing the strength of glass 
rods was described in a previous paper.” 

Briefly, a table was constructed which gave 
the stress in terms of the current through the 
armature coil for every rod diameter. The diam- 
eter of each rod was measured before inserting 
it into the quick-loading apparatus. A minimum 
stress was put on the rod for a given length of 
time ; if the rod did not break, it was then stressed 
at a higher level for the same length of time, and 
so on until it broke. The stress was increased in 
steps of 2000 !b./sq. in. if the average strength 
was greater than 12,000 lb./sq. in., and in steps 
of 1000 Ib./sq. in. if the average strength was 
less than 12,000 Ib./sq. in. Preliminary experi- 
ments enabled the operator to tell which in- 
crement to use. The average breaking stress of 
each sample of rods was corrected, as explained 
previously,? by subtracting one-half the incre- 
ment from the average. 

The procedure in testing glass rods for strength 
under vacuum was to insert the rod in the quick- 
loading apparatus, put the furnace shield into 
position and, after placing the bell jar on top 
of the iron cylinder, start the mechanical vacuum 
pump. Because of the large volume of the test 
chamber, the mechanical pump was operated for 
a half-hour or so before it became feasible to 
start the diffusion pump. After an hour and a 
half the vacuum would be down to 1 micron or 
less, at which time the rod would be baked out at 
about 350°C for 10 minutes. After baking out, 
the rod was allowed to cool for one hour because 
it lost heat very slowly. Even at the end of an 
hour, when the rod was tested for strength, the 
temperature of the rod was about 80°C. 
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TABLE I. 








Temp. Time Dura- 











Test of bake of bake tion 
No. Pressure out out of test Average strength 
1 0.05 micron Hg 350°C 10min. 10sec. 24,700 Ib./sq. in. 
2 0.2 micron Hg 350°C 10min. 10sec. 23,800 Ib./sq. in. 
3 0.1 mm Hg 350°C 10min. 10sec. 18,900 Ib./sq. in. 
4 760 mm (wet) 10 sec. 12,100 Ib./sq. in. 
RESULTS 


Results of the Vacuum, Drying, and 
Autoclaving Tests 


The strength of glass rods was measured under 
various degrees of vacuum with different heating 
treatments. These tests were aimed at discover- 
ing the effect of vacuums, and the effect of heat- 
ing plus vacuums, upon the strength and fatigue 
of glass of various chemical compositions. In 
general, comparisons were made between the 
average strengths of twenty rods. In this way a 
series of twelve experiments was carried out on 
comparable samples to determine the effect of 
various treatments on “‘strength.”’ 


Experiment I—Variation of Vacuum, 
Soda-Lime Glass 


Soda-lime rods approximately two years old 
were made into four comparable samples of 
approximately twenty rods each. The rods had 
been stored part of the time in an unheated 
warehouse and had been alternately moist and 
dry, causing a slight amount of weathering or 
bloom to form on their surfaces. Three of the 
samples were tested under vacuum as described 
above. The fourth sample was tested wet under 
atmospheric pressure with no heating. The re- 
sults are shown in Table I. 

From these results it can be seen that the 
strength is greatest in the highest vacuum. The 
strength of the glass when tested in the 0.05- 
micron vacuum is greater than the strength 
when the glass is tested wet at atmospheric 
pressure by a factor of 2.04. The tests are limited 
to ten-second duration. It is probably significant 
that the strength at 0.05 micron is slightly greater 
than at 0.2 micron, but since the probable error 
of the results is around 250 Ib./sq. in. from 
statistical fluctuations alone, this conclusion will 
have to be treated with reserve. 
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TABLE Il. 











Temp. Time Dura- 
Test of bake of bake tion 
No. Pressure out out of test Average strength 
1 0.1 micron Hg 350°C 10min. 10sec. 44,300 lb./sq. in. 
2 0.1 mm Hg 350°C §=10 min. 10sec. 32,200 lb./sq. in. 
3 10 mm Hg 350°C) =610 mink 10sec. 29,200 Ib./sq. in. 
4 760 mm Hg (wet) 10 sec. 17,700 lb./sq. in. 








Experiment II—Variation of Vacuum, 
Pyrex Glass 


Four comparable samples of carefully handled 
Pyrex rods were tested by the same procedure 
as used for the soda-lime rods in experiment I. 
The results are shown in Table II. 

The strength is greatest in the highest vacuum 
and is less as the vacuum is reduced. In this 
experiment on Pyrex glass the strength in the 
0.1-micron vacuum is 2.50 times the strength 
wet in the atmosphere. 


Experiment III—Silica Glass 


This experiment was made on two comparable 
samples of silica glass rods. The procedure was 
the same as above. The results are shown in 
Table III. 

These tests show that for silica glass the ratio 
between the strength of rods baked out in 
vacuum and the strength of those tested wet is 
only 1.67 on the ten-second test. 











TABLE III. 
Temp. Time Dura- 
Test of bake of bake tion 
No. Pressure out out of test Average strength 
1 0.05 micron Hg 350°C 10min. 10sec. 27,900 lb./sq. in. 
2 760 mm Hg (wet) 350°C 10min. 10sec. 16,700 lb./sq. in. 








Experiment IV—Porcelain 


In order to see whether or not a similar effect 
‘is obtained with porcelain, strength tests were 
made on three comparable porcelain samples. 
For the first sample the procedure was to bake 
out the test specimens in a vacuum for 10 minutes 
- and allow them to cool for an hour to about 
80°C before testing. On the second sample the 
rods were dried for one week in a desiccator in 
the presence of phosphorous pentoxide, then 
taken out of the desiccator and quickly tested 
in the test chamber at atmospheric pressure in 
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the presence of fresh phosphorous pentoxide. 
The third sample was tested wet, at atmospheric 
pressure. The results are shown in Table IV. 

From these tests it can be seen that the ratio 
between the strength in vacuum and the strength 
wet is only about 1.5. The effect is apparently 
markedly less than on glass rods. 


Experiment V—Effect of Heat, Soda-Lime Glass 


There was some question as to whether or not 
heating had much effect upon the strength of 
glass rods under vacuum. This test was made to 
compare the strength of glass rods under vacuum 














TABLE IV. 
Temp. Time Dura- 
Test of bake of bake tion 
No. Atmosphere out out of test Average strength 
1 0.01 micron 350°C 10 min. 10 sec. 10,400 Ib./sq. in. 
vacuum 
2 Dried in 10 sec. 7400 Ib./sq. in. 
desiccator 
3 Wet 10 sec. 6950 Ib./sq. in. 








with and without baking out. Since from the 
time the vacuum pumps were started to the time 
of test on the baked-out rods was about four 
hours, the sample that was not baked out was 
kept under vacuum for that time also. The 
results (Table V) were obtained on three com- 
parable samples of 20 soda-lime rods. 

These tests clearly showed that the strength 
was greatest on the 10-second tests when the 
rods were baked out in vacuum. 


Experiment VI—Effect of Vacuum on Fatigue, 
Soda-Lime Glass 


In order to study the effect on fatigue of 
different levels of vacuum, four comparable 
samples of soda-lime glass rods similar to those 
in experiment I were used. These rods had been 
stored in an unheated warehouse and had visible 











TABLE V. 
Temp. Time Dura- 
Test of bake of bake tion 
No. Pressure out out of test Average strength 
1 0.1 micron Hg 350°C 10min. 10sec. 25,950 Ib./sq. in. 
2 0.1 micron Hg None 10 sec. 21,900 Ib./sq. in. 
3 760 mm Hg (wet) 10 sec. 12,900 Ib./sq. in. 
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evidence of weathering on their surfaces. The 
results are shown in Table VI. 

These results indicate that there is no fatigue 
when the rods are baked out at 0.1 micron of Hg 
vacuum, but there is considerable fatigue at 
0.1 mm of Hg vacuum. The fact that the 
strength of the rods on the 0.1-second tests was 
greater for the lower vacuum than for the higher 
vacuum seems hard to explain, especially since 
the results could not be duplicated in experiment 
VII which was performed on another batch of 
soda-lime rods of the same chemical composi- 
tion. However, when Test No. 3 was being con- 
ducted, it became evident as soon as a few rods 
had been tested that the data were going to be 
very different from that of Test No. 1, and the 
apparatus was carefully checked without dis- 
covering anything unusual. There was no doubt 
of the increased stress the rods were bearing; 
this was shown by the amount of bending they 
underwent before breaking. Moreover, check 


TABLE VI. 








Temp. Time : 
Test of bake of bake Duration 
No. Pressure out out of test Average strength 


0.1 micron Hg 300°C 10 min. 0.1 sec. 25,900 Ib./sq. in. 
0.1 micron Hg 300°C 10min. 10 sec. 25,900 Ib./sq. in. 
0.1 
0.1 





mm Hg 300°C =: 10 min. 0.1 sec. 32,200 Ib./sq. in. 
mm Hg 300°C §=610 min. 10 sec. 21,600 Ib./sq. in. 


Pane 








tests with a few remaining canes were run at 
conditions of Test 3, which confirmed the above 
results on this batch of rods. 


Experiment VII—Effects of Vacuum and Heating 
on Fatigue 


In order to investigate the peculiarities noticed 
in experiment VI, Tests 1 and 3, a series of tests 
was made on a new batch of soda-lime rods of 
the same chemical composition as those used in 
experiment VI. These rods were only a few 
months old and showed no signs of weathering 
at the time the tests were made. Various heat 
treatments were used to see if there was any 
indication of the effects found in experiment VI. 
The results are shown in Table VII. | 

These tests confirm the results of experiment 
VI in that there is practically no fatigue when 
the glass is baked out at high vacuum: The 
tests do show however that the increased strength 
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TABLE VII. 








_Test Nature of Duration 
No. treatment Pressure of test Average strength 





1 Heated 15 minutes 0.01 micron 10 sec. 25,520 Ib./sq. in. 
at 0.01 micron 
using P20s 


2 Heated 15 minutes 0.01 micron 10 sec. 25,330 lb./sq. 
at 0.01 micron, 
no desiccator 


n. 


3 Heated 15 minutes 0.01 micron 0.1 sec. 25,809 Ib./sq. in. 
at 0.1 mm; then 
again at 0.01 mi- 
cron 


4 Heated 15 minutes 0.1 mm 0.1 sec. 25,619 Ib./sq. in. 
at 0.1 mm, tested 


cool 


5 Heated 5 minutes 0.1 mm 0.1 sec. 25,050 Ib./sq. 
at 0.1 mm 


- 


n. 


6 Heated 10 minutes 0.1 mm 0.1 sec. 23,050 Ib./sq. in. 
at 0.1 mm; held 
1 hr.; heated 19 


minutes 


~ 


No heating 0.1 mm 0.1 sec. 25,250 Ib./sq. 


n. 








on Test 3 of experiment VI was not a repeatable 
phenomenon on the newer rods, and that the 
0.1-second strength at lower vacuum was prac- 
tically the same as that at higher vacuum. It is 
our opinion that weathering had formed gels in 
the cracks that had very much influenced the 
strength tests in experiment VI. 


Experiment VIII—Effect of Dried Air 


In order to get some idea as to whether or not 
gases other than water vapor affect the strength 
of glass, four comparable samples of soda-lime 
rods were prepared. In the first test the rods of 
sample 1 were baked out and tested in vacuum 
as before. In the second test the rods of sample 2 
were baked out at 0.01 micron and allowed to 
cool for 1 hour, then the air was allowed to enter 
after it had been dried by bubbling through two 
stages of concentrated sulphuric acid. After a 
period of half an hour, the rods were tested for 
strength. The third test, with sample 3, was the 


TABLE VIII. 








Temp. 
Test of bake Duration 
No. Pressure out Pressure oftest oftest Average strength 





1 0.01 micronHg 300°C 01 micron... 10{sec.’, | 25,500 Ib./sq. in 
2 O0lmicronHg 300°C 760mmafter 10sec. 22,600 lb./sq. in. 


$ hour 

3 0O0lmicronHg 300°C 760 mmafter 10sec. 22,000 lb./sq. in. 
2%hours ~— 

4 760 mm wet 10 sec. 12,500 Ib./sq. in 
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TABLE IX. 











Test Relative Duration 

No. humidity of test Average strength 
1 0 100 sec. 11,300 Ib./sq. in. 
2 30% 100 sec. 10,200 Ib./sq. in. 
3 100% 100 sec. 9,800 Ib./sq. in. 





same as the second except that the rods remained 
in the dry air atmosphere for 24 hours before 
testing. The rods of sample 4 received no heat 
treatment and were tested wet at atmospheric 
pressure. The results are shown in Table VIII. 

Apparently even dry air causes some reduction 
in the strength of glass, as can be seen from a 
comparison of Test 1 with Tests 2 and 3. How- 
ever, it is possible that slight amounts of moisture 
were able to come through the drying apparatus 
and of course there must have been some SO; 
and H.SO, gases present. 


Experiment [X—Effects of Dry Storage 


It was thought advisable to obtain compara- 
tive figures on the strength of “dry” glass rods 
as compared with glass under average atmos- 
pheric conditions and as compared also with 
thoroughly wet glass. The ‘‘dry’’ rods were dried 
in a desiccator for a week in the presence of 
P.O, before testing. The results are shown in 
Table IX. 

An increase of strength of about 15 percent 
above the strength when wet is indicated when 
rods are dried in a desiccator, and an increase 
of about 5 percent when the relative humidity is 
around 30 percent. However, this work was not 
too carefully controlled and is included merely 
to show approximately what can be expected. 


Experiment X—Effect of CO., Soda-Lime Glass 


A few tests had indicated that some factor 
other than water was a weakening agent. We 
suspected CO, because it is known that CO, in 
weak solution will dissolve silica more than will 
plain water. Four comparable samples were 
used. Each sample was put into a different liquid 
for two hours previous to testing and kept wet 
with the liquid during test. The liquids used 
were distilled water, distilled water carbonated 
at room temperature, distilled water with satu- 
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rated Ca(OH)s, and distilled water with normal 
NaOH. It was felt at this time that the differ- 
ences of strength observed were due more to 
differences of fatigue rates rather than to a 
fundamental lowering of the “strength.”” There- 
fore the rods were tested for 100 seconds duration 
in hopes that differences would be more apparent. 
The results are shown in Table X. 

There is some indication that the CO, is an 
additional weakening agent besides the water. 
If so, the rest of the data might be explained _ on 
the basis that the CO, that is normally in the 
air was kept away from the flaws. However, 
until more carefully controlled tests are available 
it will not be possible to be absolutely sure about 
the effect of CO». 


Experiment XI—Effect of Etching by 
Water Vapor 


Since the attack of water on glass appeared to 
be the chief cause of fatigue and the chief 
weakener of glass, some studies of the effect of 
water upon glass at an elevated temperature 
were tried. A high pressure autoclave allowing 
use of steam up to 1000 Ib./sq. in. was used to 
make the tests. It was found that under pro- 
longed treatment at pressures as high as 800 
lb./sq. in. all glasses, including silica glass, 
were attacked severely by water. The soda-lime 
glasses were attacked the most, Pyrex glasses 
next, and silica glass least of all. The attack was 
similar to etching with hydrofluoric acid. In the 
case of the soda-lime glasses a very thick jelly- 
like layer formed on the surface. Since the etching 
was similar to that of hydrofluoric acid, one would 
expect the rods to get stronger when autoclaved. 
Unscratched silica rods and rods scratched with 
coarse emery were autoclaved at 700 to 800 
lb./sq. in. for twenty-four hours and then re- 
moved and tested for strength (Table XI). 








TABLE X. 
Test Duration ‘ 
No. Liquid surrounding rod of test Average strength 
1 Distilled water 100 sec. 13,400 Ib./sq. in. 
2 Distilled water satu- 100sec. 13,150 lb./sq. in. 
rated with CO, 

3 Saturated solutionof 100sec. 13,950 lb./sq. in. 
Ca(OH)2 

4 Normal solution of 100sec. 13,750 lb./sq. in. 


NaOH 
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From these tests it is immediately evident 
that autoclaving can increase the strength of 
silica glass, probably by etching out the flaws. 

Tests similar to those on the silica rods were 
made on soda-lime glass rods. In this case the 
rods were autoclaved for four hours at 600 Ib./ 
sq. in. then removed and tested wet in the quick- 
loading apparatus. The results are shown in 
Table XII. 

These tests on soda-lime rods showed that 
their strength also could be increased by etching 
in an autoclave. 


Experiment XII 


Some further tests were made on the breaking 
strength of rods in an autoclave. A glass rod was 
supported by knife edges and loaded centrally 
with a heavy metal weight. The whole apparatus 
was put into the autoclave and the steam pres- 
sure was applied for about one hour. Then it 
was observed whether or not the rod had broken. 
When the rod did not break in the autoclave, it 
was next tested in the open and found to be 
very much stronger than comparable rods that 
had not been autoclaved. But such rods as did 
break in the autoclave broke at unusually low 
stresses of 2000 Ib./sq. in., and even less, com- 
pared with about 5000 Ib./sq. in. for the lowest 
breaking strength of comparable rods when 
tested for 24 hours in the open. It would appear 
therefore that glass can fatigue very much more 
in an autoclave than in ordinary atmospheric 
conditions. 


Discussion of Results of Vacuum Tests 


These experiments on the strength of glass in 
a moderately high vacuum show clearly that the 
strength of glass is much greater in vacuum than 
when the glass is wet and in the open, and that 
fatigue is greatly reduced or entirely eliminated 
by vacuum treatment. 

The increase in strength due to baking in a 














TABLE XI. 
; Time Steam Duration 
Test No. autoclaved pressure of test Average strength 
1-Scratched 24 hours 700-800 10 sec. 30,000 Ib./sq. in. 
2-Scratched 0 10 sec. 6000 Ib./sq. in. 
3-Unscratched 24 hours 700-800 10 sec. 54,000 Ib./sq. in. 
4-Unscratched 0 10 sec. 19,000 Ib./sq. in. 
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TABLE XII. 








re Duration ; ; 
Test of auto- Pressure of Duration 
No.  claving autoclave of test Average strength 





5 4hours 600 Ib./sq.in. 10sec. 27,400 lb./sq. in. 
6 0 10 sec. 13,000 Ib./sq. in. 








vacuum is shown by experiments I-IV. The ratio 
of the vacuum strength to the wet strength is 
2.5 for Pyrex, 2.0 for soda-lime glass, 1.7 for 
silica glass, and 1.5 for porcelain. Thus it would 
appear that this effect is dependent somewhat on 
chemical composition. The fact that it is rela- 
tively low for porcelain and silica glass (which 
are chemically resistant) would suggest that the 
effect is correlated with the chemical attack of 
moisture or gases at the tip of the crack. 

It is quite evident that glass fatigues very 
little, if at all, in a high vacuum. Tests 2 and 3 
of experiment VII show that the strength is 
25,800 Ib./sq. in. on the 0.1-second test and 
25,300 Ib./sq. in. on the 10-second test. Also, 
Tests 1 and 2 of experiment VI give the strength 
as 25,900 on both the 0.1-second and 10-second 
tests in high vacuums. The fact that there was 
very little fatigue at high vacuums was evident 
from observations with a cathode-ray oscillo- 
graph of the times required to break the rods. The 
rods generally broke as soon as the load came on 
or not at all; only occasionally was there definite 
evidence of delayed breakage. However, unless 
the rods were baked out at 350°C or more, at 
high vacuum, there was considerable evidence of 
fatigue. It would therefore appear that fatigue is 
caused primarily by the action of water vapor 
and adsorbed gases acting in the flaws. It is very 
doubtful that even heating the glass up to 
350°C or more removed all the water vapor 
from the flaws, since they were in the neighbor- 
hood of 0.005-inch deep in the glass used and 
probably only a few Angstrom units wide at the 
tips. Too little theoretical knowledge is available 
to estimate the length of time required to evacu- 
ate the water vapor under those circumstances. 
It is possible that the most sensitive test of the 
presence of water vapor in the tip of a flaw is the 
occurrence of fatigue in strength tests. 

There is considerable evidence that the 
“strength” of glass is a much more complicated 
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matter than is generally realized. Such evidence 
is apparent in’a comparison between the tests in 
experiments VI and VII. In both cases the rods 
were of the same composition. The rods in 
experiment VI had been stored for a few years, 
some of the time in an unheated warehouse, and 
showed slight weathering. These rods showed a 
strength of 32,200 lb./sq. in. on the 0.1-second 
test in a 0.1-mm vacuum, as compared with 
25,900 lb./sq. in. for rods baked out in a vacuum 
of 0.1 micron. Careful checks were made so 
there was no doubt of the phenomenon, yet it 
was unrepeatable in experiment VII on glass 
that was not so old and had been stored in a 
heated room most of the time. 

The only hypothesis to explain these rather 
erratic results seems to be that gels form in the 
flaws which strengthen the glass but have much 
greater fatigue rates than glass. When the glass 
is thoroughly baked out, the gels crack or lose 
their strength to an extent that makes them a 
negligible factor. 

The chemical attack of water on glass at 
elevated temperatures and pressures was shown 
in experiments XI and XII. It was somewhat 
surprising to see the rapidity of attack of water 
on silica glass at steam pressures around 700 to 
800 pounds per square inch. That the flaws were 
etched out and that the strength of the glass 
was increased was clearly demonstrated. It was 
surprising to find that the process of fatigue could 
be remarkably speeded up at these temperatures, 
so that glass fatigues to much less than its 24- 
hour strength in a matter of an hour or possibly 
only a few minutes. These results strengthen 
our feeling that the process of fatigue is essen- 
tially a chemical attack. 

The chemical attack upon glass by moisture is 
extensively discussed in the literature and sum- 

marized by Morey.’ That a chemical attack 
inside the flaws should decrease the strength of 
glass is not surprising, as one would expect the 
bonds under the greatest tension to be those 
most ready to rupture and combine to form new 
* compounds. This chemical attack will of course 
depend on two factors, the first being the ability 
of the water to penetrate to the tip of the flaw 
and the second being the chemical durability of 


5G. W. Morey, The Pro 


rties of Glass (Reinhold 
Publishing Corporation, New 


ork, 1938), Chap. V. 
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the glass. The fact that all cracks are pretty 
much alike dimensionally as well as the fact that 
all glasses, including silica, are attacked by water 
probably accounts for the similarity of strength 
between glasses at room temperature. That there 
are differences in the amount of fatigue in glasses 
of different compositions and that the ratio 
between the strength of glass in vacuum and the 
strength of glass wet is different for different 
glasses indicate that there is a partial correla- 
tion between chemical composition and fatigue 
of wet glass. It was found? that silica appears to 
fatigue less than soda-lime and Pyrex, and the 
Pyrex appears to fatigue even more than soda- 
lime if one takes into account the fact that the 
samples of Pyrex glass had an initial com- 
pression. Then of course in the vacuum work 
just cited the ratio of the vacuum strength to 
the wet strength was greatest for Pyrex, with 
soda-lime following and, finally, the ratio for 
silica glass was least of all. In general, Pyrex 
glass is considered to be more chemically durable 
than soda-lime, but this is not necessarily true 
when the attacking medium is basic in nature. 
In a flaw there is no reason for the alkali to be 
leached out. Since probably the attack in both 
Pyrex and soda-lime glass is basic in nature the 
correlation with chemical durability seems to 
hold. 

These tests on the strength of glass under 
different conditions of the surrounding medium 
give considerable insight into the ultimate ques- 
tions pertaining to the fatigue of glass under 
constant load. Specifically, it seems to be neces- 
sary to distinguish between glass in which the 
flaws are of recent origin and glass in which 
the flaws have existed for some time and have 
been exposed to the atmosphere. 

In the first case the glass has a definite strength 
under high vacuum no matter what the duration 
of the load (at least to the extent covered in 
our tests). When even slight amounts of moisture 
are present the glass loses strength with duration 
of load, although this loss is slight in moderate 
vacuum. As the humidity is increased,. the loss 
of strength with duration of load becomes greater, 
that is, as the absolute vapor pressure of the 
water is increased a condition arises where fatigue 
is very rapid. Thus it appears that when the 
duration of load is decreased the strength of 
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glass approaches its vacuum strength. This sug- 
gests that in extrapolating fatigue curves to 
shorter times it would be unjustified to go to 
durations which yield strengths greater than the 
vacuum strength. Moreover, we should expect 
the fatigue curves, now extended experimentally 
only to 0.0i-second duration, to level off as the 
duration of load is decreased to around 0.0001 
second. In all probability the extrapolation of 
the present fatigue curves to shorter times is 
not legitimate. 

In the second case, where the flaws have had 
time to weather and gels have formed inside 
them, the picture of the short time strength is 
too complicated to permit any definite statement. 
The short time strength probably is determined 
more by the strength of the gels than by the 
strength of the glass. Some glass may act like 
freshly scratched glass and some may show 
behavior characteristics similar to the weathered 
glass of experiment VI. The very short time 
strength can be equal to or greater than the 
vacuum strength, depending on the gels. 

In the middle range of durations of load (from 
0.01 second to perhaps several years) the decay 
of strength with time of wet glass at room tem- 
perature and atmospheric pressure may be repre- 
sented by the equation 


S= [Cy /log Cot J+ C3. 


This equation may be regarded as empirical, 
and other empirical formulae fit the results 
equally well. One such equation is discussed by 
Preston and Glathart in the fourth paper of 
this series, entitled ‘“‘The Fatigue Modulus of 
Glass.’’ 

Under ‘‘atmospheric”’ conditions, it is likely 
that there is very little difference in fatigue 
between glass with weathered flaws and glass 
with fresh flaws. But this may depend on the 
way we define ‘fatigue’; see the paper of 
Preston and Glathart mentioned above.® 

If however the conditions of the atmosphere 
are changed the picture is greatly altered. As the 
temperature and water vapor pressure are in- 
creased the fatigue is rapidly accelerated until 
at pressures of a few hundred pounds of steam 


°F. W. Preston and J. L. Glathart, J. App. Phys. 17, 189 
(1946). 
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a decline in strength takes place in a few minutes 
that might take hours or years at room tem- 
perature. However, if the glass is not broken the 
flaws will etch out in a few hours and the speci- 
men will become very much stronger than before. 
Ultimately, the attack of the steam will destroy 
the glass rods completely and so we may con- 
sider the strength to be zero for glass subjected 
to 1000 Ib./sq. in. steam pressure for a week 
or two. 

If a glass rod were under constant load and 
surrounded by water at room temperature, it is 
exceedingly probable that its strength would 
behave in a manner similar to that of glass 
subjected to elevated temperatures in an auto- 
clave. For fairly high stresses the tips of the flaws 
would tend to round out without much extension 
and the glass would eventually become stronger. 
Ultimately, the glass would dissolve and be 
destroyed by the water, but many thousands of 
years would be required in the case of chemically 
resistant glasses. 

It is a moot question as to what will happen 
to glass under load at room temperature in a 
dry atmosphere (relative humidity between 1 
percent and 50 percent). There is evidence of 
chemical attack on most glasses under these 
conditions. The glass will fatigue in about the 
same way as it does when it is wet. Therefore, 
it is likely that for many glasses the ultimate 
results will be the same. Under high stresses the 
flaws will extend until breakage occurs. Under 
low stresses of about 1000 to 2000 Ib./sq. in. the 
‘attack of the water will tend to round out the 
flaws and the glass will slowly get stronger. On 
the other hand, crystalline material may form, 
in the case of some of the less durable glasses, 
tending to force open the flaws and progressively 
extend them. In such a case the strength of the 
glass is zero for durations of load as short as a 
few years.! 


CONCLUSIONS 


In high vacuum, thoroughly baked out glass 
can be said to have a tensile strength fairly inde- 
pendent of the time of loading, at least for the 
shorter durations of load. If glass is exposed to 
moisture it rapidly loses strength under stress. 
In ten seconds the strength can decay to } or 2 


of the vacuum strength, and it may decay to } 
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the vacuum strength in 24 hours. This strength 
decay appears to be caused by a chemical attack 
in the flaws of the glass by water and possibly 
by some gases such as COs. 

It therefore appears that for very short dura- 
tions of load the strength of freshly scratched 
glass approaches the high vacuum strength. For 
moderate durations of load ranging from 0.01 
second to probably several days or so the 
strength of wet glass or glass under humid 
conditions approximately obeys the relationship 


S=[C,/log Cot ]+Cs, 


where S is breaking stress, ¢ is time, and C,, Co, 
and C; are arbitrary constants. For longer dura- 
tions of load, several years or more, the chemical 
composition is so important as to be the con- 
trolling factor in strength, but there is reason 
to believe that at low stresses the flaws can etch 
out and the glass actually become stronger, in 
the case of chemically resistant glasses, so that 
glasses that survive for years may survive for 
several millions of years. Inevitably of course the 
chemical attack of the atmosphere will disin- 
tegrate the glass to a point where the strength 
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is zero, in the case_of silica glass in probably 
billions of years. 
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The Fatigue Modulus of Glass 


J. L. GLatHaRT AND F. W. PRESTON 
Preston Laboratories, Butler, Pennsylvania 


(Received November 1, 1945) 


The decline in the breaking strength of glass with an increase in the length of time it is 
under stress has been noted by various experimenters. Based on the experimental results of 
T. C. Baker, an empirical relationship is obtained which indicates that the reciprocal of the 
breaking stress is a linear function of the logarithm of the duration of the stress. Certain 


implications of this relationship are discussed. 


I. INTRODUCTION 


ARIOUS experimenters have shown that a 
glass specimen may sustain a certain static 
stress momentarily but will ultimately fail under 
the same stress, with the mere passage of time. 
No complete theoretical basis for this effect has 
as yet been established, although Murgatroyd! 
has recently proposed a “‘viscous pocket”’ theory 
for the constitution of glass which, he states, is 
capable of accounting qualitatively for the ob- 
served decrease in strength with time. 

The object of the present paper is to present an 
empirical relation which fits the strength-time 
curves of all the glasses for which adequate 
experimental data are available. Some of the 
implications of this relation are discussed. 


Il. THE STRENGTH-TIME RELATION 


In another paper? T. C. Baker has given his 
experimental results of the fatigue effect in fused 
silica, three kinds of commercial glasses, and 
three porcelains. His results are summarized in 
Table I. 

The total range of durations of stress is seen to 
be from 0.01 second to 24 hours (86,400 seconds), 
a range of nearly ten million to one. This is a 
tremendous range, and it becomes virtually 
necessary to plot it on a logarithmic base. The 
various ranges of strengths, however, are rela- 
tively small, from 23,880 to 11,740 Ib./sq. in. in 
the case of fused silica, for instance, and from 
19,900 to 6,450 Ib./sq. in. in the case of soda-lime 
glass. It is natural to plot this on a linear or 
arithmetical scale. This method of plotting, f 
against log t, was the one adopted by Smekal* 

! J. B. Murgatroyd, J. Soc. Glass Tech. 28, 406 (1944). 

2 T. C. Baker and F. W. Preston, ‘Fatigue of glass under 


static loads,” J. App. Phys. 17, 170 (1946). 
3 A. Smekal, Ergeb. d. exakt. Naturwiss. 15, 106 (1936). 
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who, over a smaller range of durations, reported a 
linear relationship. This came about, as James 
Bailey pointed out later, because Smekal inad- 
vertently misplaced some of his points by a factor 
of ten. Had they been correctly plotted, his line 
would not have been straight. 

Baker adopted the same rather natural method 
of plotting, and obtained very definitely curved 
lines for all of his materials, the curvature being 
the more obvious because of his longer range of 
time intervals. Baker, with commendable cau- 
tion, refused to indulge in any speculation re- 
garding the strengths of his materials beyond the 
experimental time range, but warned that any 
point on any of his curves might be in error by as 
much as +10 percent. 

The Preston Laboratories have received vari- 
ous inquiries and suggestions regarding this 
matter of the fatigue of glass, especially since the 


TABLE I. Average breaking stress in pounds per square inch 
(and its reciprocal)* for various durations of load. 











0.01 0.1 1.0 10 100 1000 86,400 
Material sec. sec. sec. sec. sec. sec. sec. 
Disannealed 24,600 19,600 17,100 15,900 _ 12,800 11,300 
Pyrex (dry) (4.06) (5.10) (5.85) (6.30) — (7.81) (8.85) 
Scratched 12,900 10,600 9,000 8,250 — 5,750 5,050 
Pyrex (wet) (7.75) (9.43) (11.11) (12.12) _— (17.4) (19.8) 
Annealed 19,900 15,900 — 10,500 -- 7,950 6,450 
soda-lime (5.03) (6.29) — (9.53) _ (12.58) (15.50) 
glass (wet) 
Annealed lead — 19,000 14,900 11,200 _ 9,400 7,950 
glass (wet) —_— (5.26) (6.71) (8.93) — (10.63) (12.58) 
Fused silica _ 23,880 _- 19,800 15,660 14,280 11,740 
(wet) _ (4.19) _— (5.05) (6.39) (7.00) (8.52) 
Porcelain A 15,800 13,900 12,300 — 9,700 “= 9,100 
(dry) (6.34) (7.20) (8.13) — (10.32) — (11.00) 


Porcelain B 15,400 13,500 11,900 9,900 9,400 9,100 8,900 
(dry) (6.50) (7.41) (8.40) (10.10) (10.64) (11.00) (11.23) 














Porcelain C 10,900 9,400 8,400 6,900 _ 6,000 5,900 
(dry) (9.17) (10.64) (11.91) (14.50) _ (16.67) (16.95) 
* Reciprocals are given in parenthesis, and are to be multiplied by 107°. 
4]. Bailey, private communication. 
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publication of the Philadelphia lecture,> and 
the matter has recently been considered by 
Murgatroyd and Haward.* In view of this inter- 
est, it seems desirable to discuss Baker’s curves a 
little more fully with respect to whatever em- 
pirical information may be obtained from them. 

We have known for some time that if 1/f, the 
reciprocal of the breaking stress, is plotted 
against log ¢ a straight line can be drawn which, 
for every one of Baker’s materials, passes within 
the “‘possible error’ he assigned to his values. A 
note to this effect appeared in Nature.’ This 
relation has also been noted by Haward,* whose 
observations reached us after the note in Nature 
had been published. 

Figures 1 to 8 show the straight lines that 
result when Baker’s values are plotted in this 
way. Figure 9 is a similar plot of some previous 
data reported’ from this laboratory, which were 
obtained on plate glass laths of dimensions 
16’ 2""X1"’, considerably larger than Baker’s 
round rod specimens. Although they cover a 
much shorter time range, these data are included 
because they were cited by Haward.*® 

In the case of two materials, fused silica and 
soda-lime glass, this reciprocal relationship is so 
exact that nearly every point lies on the straight 
line. In the case of the Pyrex glass and the lead 
glass the fit is also quite good, but the points for 
the porcelains deviate considerably, although 

-even in these cases straight lines can be drawn 
within the errors assigned. The discussion pre- 
5 F. W. Preston, J. App. Phys. 13, 623 (1942). 
®R. N. Haward, Comment on Murgatroyd’s paper. See 
reference 1, p. 427. 


7F. W. Preston, Nature 156, 55 (July 14, 1945). 
*F. W. Preston, J. Am. Cer. Soc. 18, 220 (1935). 
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sented in this paper will apply more particularly 
to the four glasses, Pyrex, soda-lime, lead, and 
silica. 

Since the straight lines fit so well in these cases, 
it is natural to inquire as to their equations, and 
whether any theoretical basis can be found for 
them. We find that if there is any theoretical 
basis at all it is expressed in terms of concepts 
such as “‘weakness,”’ rather than strength, and a 
“fatigue modulus.”’ 
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It is more or less inevitable, also, that we 
should speculate on the extrapolation of the 
curves, to time intervals both longer and shorter 
than were attained in Baker’s experiments. There 
is a practical necessity for speculating on the 
longer times, and fortunately there is no apparent 
reason why the curves should not continue as 
straight lines to times far beyond anything of 
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technological interest. At the short time end 
things are more complicated. 


III. EXTRAPOLATIONS 


The various curves (straight lines) which repre- 
sent Baker’s results are expressed by equations 
of the form 


1/f=a+m log t, (1) 


where f is the breaking stress, ¢ is the time, and 
a and m are constants. Each line is thus defined 
by its slope d(1/f)/d(log t) and by the point (fo), 
where the line intersects the time axis. 

With such materials, once the fact that they 
plot to straight lines is conceded, the question 
arising in the mind of the experimenter—or 
reader—is apt to be: “how far can the curves be 
extrapolated at each end with any confidence?” 


(a) The Long Time End 


The equation above would make the strength 
zero at infinite time. This is psychologically 
satisfying, and, furthermore, is required by some 
theories of strength which have been recently 
developed. Otherwise stated, it means that even 
a slight stress will break the glass if only it lasts 
for a sufficiently long time. 

There seems to be nothing in the physics of 
long time loading that would cause us to deny the 
possibility that the straight-line law holds more 
or less indefinitely, provided temperature and 
humidity are constant and there is no physical 
change in the specimen other than changes which 
operate during our relatively brief loading inter- 
vals. The straight-line relation holds from 0.01 
second to nearly 10° seconds, a range of time 
durations of 10’ to 1. Our experiments end at a 
duration of one day, but there seems no reason 
why we should expect any departure from the 

‘linear relation if we extend the time range by a 
factor of 10*, that is, from one day to three years, 
or even by a factor of 10‘, that is, up to thirty 
years. This is much longer than most bottles last, 

~ and for most of their lives they are not under 
stress. There are not many examples of plate 
glass, or other kinds of glass, that are expected to 
last, under stress, for so great a length of time. 

If the curve were extended to the right, beyond 

one day, as far as we have determined it experi- 
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mentally to the left, that is, in the ratio of 10’ 
to 1, we arrive at a period of 30,000 years, or 
roughly the time that has elapsed since the 
maximum of the last (‘‘Wisconsin’’) glaciation, 
and several times as long as all recorded history. 

In Table II below we take the strength at the 
end of one day of stress as unity and calculate, 
for fused silica and soda-lime glass, the pro- 
portionate strengths at several other longer 
durations, on the assumption that this linear 
relationship holds indefinitely as we move to the 
right. 


(b) The Short Time End 


The strength-time equation would make the 
strength infinite for some very brief, but not 
infinitesimally short, duration (to). This point is 
obtained by extending the straight line of the 
graph to the left until it intersects the time axis. 
The values of to for the various materials tested 
are shown in Table III. 

This threshold time, fo, is obtained from the 
equation of the curve and its value depends, 
therefore, on the way the straight line is drawn 
among the experimental points. It is interesting 
to note, however, that for all those glasses which 
were tested wet, the value of fo is in the neighbor- 
hood of 1 microsecond, whereas its value for those 
materials which were tested dry is very much 
less. However, other factors than moisture are 
involved, of course. For example, porcelains A 
and B, which had identical compositions but 
different firing temperatures, have nearly the 
same value of to, whereas porcelain C, of different 
composition, has a t) value 20 to 40 times as 
large. 


TABLE II. Extrapolation at the long time end. 








Infi- 
2X10° nite 

















‘Duration 108 106 
Material iday iyear years years years time 
Fused silica 1.0 0.82 0.68 0.57 0.50 0 
Soda-lime glass 1.0 0.80 0.65 0.55 0.46 0 
TABLE III. The threshold time. 

Material to(sec.) Material to(sec.) 
Soda-lime glass 5.0X10-® Pyrex (dry) 0.0085 x10-* 
Lead glass 10.2X10-* Porcelain A 0.00004 x 10-6 
Fused silica 0.2X10-* Porcelain B 0.00008 x 10-¢ 
Pyrex (wet) 0.4X10-* Porcelain C 0.0016 «10-6 
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TABLE IV. Theoretical strengths at short durations. 








Duration 10-2 10-3 10-4 10-5 10-6 
Material, sec. sec. sec. sec. sec. 
Fused silica 28,200 37,000 50,700 80,700 196,000 


Soda-lime glass 19,900 28,800 51,300 232,000 © 


At present, to is merely a mathematical ab- 
straction. It could not be even remotely ap- 
proached experimentally with such a set-up as we 
were using. Yet on our logarithmic scale it is not 
very far beyond the short time end of the 
experimentally determined curve. In fact, for 
some of the materials it is much closer to the left- 
hand end of the experimental curve than this 
left-hand end is to the right-hand end, so the 
extrapolation, considered as a mathematical 
operation, is not too extravagant. 

Physically, however, Baker probably carried 
the curve toward the short time end as far as it 
can be carried unless some totally different and 
much smaller set-up can be devised. Possibly 
some arrangement using a silver plated glass fiber 
supported on a short span in a magnetic field, and 
deflected by the discharge of a condenser through 
the silver plating, might serve to reach shorter 
time intervals. 

Let us assume that by some means or other the 
curve could be extended, experimentally, to 
shorter times, and that it did, in fact, continue 
straight for some distance. The calculated 
strength values of fused silica and soda-lime glass 
at several shorter times are shown in Table IV. 

The behavior of glass under very high stresses 
cannot be expected to parallel its behavior under 
moderate stresses, because we cannot expect 
Hooke’s law to be obeyed at extreme strains. It 
seems not inconceivable that the straight-line 
law might apply fairly well until stresses of a 
million pounds per square inch are reached, with 
a resulting strain or extension of 10 percent in the 
length of the fiber. But it is scarcely likely that 
it could apply at stresses of ten million pounds 
per square inch, a figure equal to Young’s 
modulus as calculated for small strains. 

At the short time end the inertia of the ma- 
terial becomes an important factor and makes 
experimentation impractical, by any method as 
yet envisaged, much beyond the point to which 
Baker carried it. 

Thus, for more reasons than one, the shape, or 
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even the physical reality, of the curves beyond 
Baker’s short time values must remain doubtful 
for the present, although the theoretical impli- 
cations may make us wish to find means of 
extending them in that direction. This might be 
accomplished by totally different techniques, 
perhaps by impact testing, as hinted by Baker. 


IV. THE CONCEPT OF “STRENGTH” OF GLASS 


From the empirical equation 
feto/f =log, (t/to), (2) 


which can be deduced from Baker’s experimental 
results, it can be seen that the strength f is not 
a constant of the material but that it depends 
upon a variable (time, ¢) and on two constants of 
the material (or at least of the specimen). One of 
these is fo, which is the maximum duration of 
stress for which the strength is theoretically infi- 
nite, and the other is feto, which is the strength of 
the glass for a duration of stress ety, where 
e=2.718, the base of natural logarithms. 

We may tentatively call ft the ‘threshold 
duration.”’ It is the least interval of time in which 
a tensile stress, however large, can rupture the 
material (or the specimen). In our experiments on 
glass it is of the order of a microsecond, and in 
this time a crack can be propagated no more than 
60,000 X 10~* inches = ;'g inch, and the stress pulse 
no more than 4 inch. These dimensions seem to 
bear no relation to the ultimate structure of the 
material, but they are comparable with the 
shortest dimension of the specimen (0.22 inch), 
or with the shortest dimensions of an effective 
crack. They may even correspond somewhat 
closely with the radius of the smallest ‘mirror’ 
(initiation crack) that can be produced under 
infinite stress. This in turn may be, and probably 
is, a function of the size of the specimen, and not 
a true constant of the material. Obviously, the 
mirror on a fine fiber has to be smaller than 
ve inch. 

If to is a property of the specimen rather than 
the material, then so is efo, and likewise, pre- 
sumably, is the strength at this time interval, 
or feto. 

Equation (2) above merely expresses an 
equality of ratios; it describes the “fatigue,” but 
it says nothing about strength itself. The strength 
at one duration is expressed only in terms of the 
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strength at some other duration, not in terms of 
the physical constants of the material, or the 
dimensions of the specimen: Until to and feto can 
be so described, the equation is limited in exactly 
the same way as the experiments upon which it is 
based. 

By differentiation of (2) we get 


ds 1 
fuw—(- )- ‘ (3) 
dt\ f t 


This eliminates ¢) from the picture but does not 
simplify things very much. The reciprocal of the 
breaking stress, that is, 1/f, may be called the 
“weakness” of the specimen; it is the area that 
must be assigned to carry unit load, if the 
specimen is to last out the time interval ¢. If this 
is denoted by w, then 


1 1 
=w and 


f feto 
dw dt = Weto(1 t). (5) 


=Weio. (4) 


So 


The rate of increase of weakness, dw/dt, we may 
call “fatigue.’’ If so, Eq. (5) shows that fatigue is 
a property which varies inversely with the 
duration of stress. 

Equation (2) may now be written 


t to =e” Wety, (6) 
and we then have, from Eq. (5), 
dt to 


—_= -—e" Wety, ( 
dw  Weto 


~~ 
— 


t is the “lastingness,”’ or endurance, of the speci- 
men and is here expressed in terms of w, the area 
assigned to carry unit load. 

We may write, finally 


1 dt 1 
— a> @ <-> S85 Tele. (8) 
t dw Wety 


This says that the relative gain in endurance (as 
compared with the previous total endurance), 
due to unit increase in assigned area, is a constant 
of the material, or at least of the specimen. 

In other words, if we want to increase the 
endurance by 1 percent we must add a constant 
number of square inches of cross-sectional area to 
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TABLE V. Fatigue modulus. 


Modulus Modulus 

Material in.?/Ib. Material in.?, Ib. 
Pyrex (dry) 3.0 1076 Fused silica 3.2 1076 
Pyrex (wet) 1.6 10-6 Porcelain A 3.3 10-6 
Soda-lime glass 6.5X 10-6 Porcelain B 2.4X 10-6 
Lead glass 5.7X 10-6 Porcelain C 5.9 10-6 


carry unit load. The area, dw, to be added is 
inversely proportional to the “discriminant 
strength,” feto, since, from Eq. (8), fetoXdw is 
constant for a given value of dt/t. 

If, instead of talking about a 1 percent incre- 
ment in the life of the specimen, we define 
‘fatigue’ as that number of square inches that 
must be added to the cross-sectional area in order 
to double the life of the specimen, then 6¢=¢ in 
Eq. (8) and 


W = Wety=1/ fety. (9) 


In this case, 6w is the fatigue in square inches per 
pound of load. 

For the various glasses tested the éw values in 
in.?/Ib. are given in Table V. 

It should be understood that éw is to be de- 
duced from Eq. (8), not ascertained by de- 
termining experimentally what area must be 
added to a specimen to increase its life from, say, 
one minute to two minutes. dw is a kind of 
“modulus,” in the sense that Young’s modulus E 
is ‘“‘that stress which would double the length of 
the specimen.”” In practice, we do not try to 
double the length of the specimen. Rather, we use 
a small increment of length, assume Hooke'’s law, 
and deduce E from the equation 


(1/L)(dL/df) =1/E. 
This compares with our Eq. (8) 
(1/t)(dt/dw) = 1/ wet. 

Subject to this understanding, it seems we can 
say : Fatigue is that number of square inches that 
must be added to the cross-sectional area assigned 
to carry unit load (one pound) in order to double 
the life of the specimen. 


V. PHYSICAL CONSIDERATIONS 


The above discussion consists merely in the 
manipulation of the mathematical symbols in an 
attempt to interpret the straight-line relationship 
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derived from Baker’s curves in terms of concepts 
which may find further application in this field. 

There are other considerations however, some 
of which have been outlined. Another is the 
inherent unlikeliness that the lines representing 
the soda-lime glass and the silica glass actually 
cross one another to the left of the observed 
points, which would mean that the soda-lime 
glass becomes stronger than fused silica for 
stresses of sufficiently short duration. It seems 
more likely that, if the curves have any physical 
reality when extrapolated far to the left, they do 
not remain straight, but presumably flatten 
rapidly and become asymptotic to the time axis. 

In all probability, this means in turn that our 
experimental straight-line relationship ‘‘in the 
visible’ is itself only an approximation, even 
though it seems an extraordinarily good one in at 
least some of the instances. 

Haward’s theory® leads to an equation 


S(e?0/s —e) = kt, 


where S is the applied stress, t is the time required 
for the glass to break under this stress, and Po 
and & are constants. For t=0, this yields S= Po, 
that is, the strength is not infinite, as would be 
required if our straight line held indefinitely. 
On the other hand, Murgatroyd’s viscous 


pocket theory! implies a relation identical with 
our Eq. (1). It is obvious that much more 
theoretical work needs to be done on the matter 
of the meaning of the ‘‘strength”’ of glass and on 
the relation of strength to time. 


CONCLUSIONS 


The ordinary concepts of “strength” do not 
apply to glass, at any rate not to glass as we 
know it, exposed to the earth’s atmosphere and 
to atmospheric moisture. Under such normal 
conditions the “life” or “‘endurance’”’ of a stressed 
glass specimen depends on the cross-sectional 
area which we assign to carry unit load. The 
greater the area we assign, the longer the life of 
the specimen, the life, or endurance, varying 
exponentially with the assigned area. 

This permits the calculation of a certain 
constant, which we have called the “fatigue 
modulus,’ which describes the behavior of the 
specimen under the particular conditions of the 
experiment. 

However, in all probability, this modulus is not 
so much a property of the substance of the 
specimen as it is a function of the surrounding 
atmosphere or, more accurately, of the conditions 
obtaining at the surface of the specimen. 
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This paper deals with the heat input to hollow cylindrical metal cores subjected to induction 
heating. The rigorous expression for the heat input, valid for the whole frequency range, is 
developed and is presented by Eq. (27) of this paper. It is shown in an example how the equa- 
tion is used for numerical computations. Existing approximate formulas for the low frequency 
and the high frequency range are checked by means of the rigorous equation; the maximum 
deviation between the rigorous and approximate equations is about 10 percent. 


1. THE PROBLEM 


S a result of the growing number of practical 
applications of induction heating of metals, 
many papers have been published that present 
the theory of induction heating and its usefulness 
in solving practical problems. One of the chief 
applications of the theory involves the calcula- 
tion of the heat input to a metal core surrounded 
by a coil carrying alternating current. In the 


’ VOLUME 17, MARCH, 1946 


case of a solid core rigorous equations are usually 
given, from which are deduced the simplified 
equations that hold with good approximation, 
if the frequency is either higher or lower than 
certain critical limits. 

When hollow cores are discussed, the rigorous 
equations referring to that case are not given. 
It is occasionally pointed out how these may be 
obtained, but the final equation is not presented. 
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Instead, only the approximate equations for 
frequencies above and below certain limits are 
given. This type of treatment can be found, for 
instance, in papers by Greig,’ Wall,? and Baker,’ 
as well as in a book by Warren.‘ 

While the approximate equations just referred 
to seem to be sufficient for practical purposes, 
it appeared to the author that a presentation of 
the rigorous equations would be valuable for the 
following reasons: 

(1) The practical importance of the problem 
certainly justifies knowledge of the rigorous 
solution. 

(2) A reliable check on the validity of ap- 
proximate solutions can best be made by means 
of the rigorous solution. 

(3) It might well happen that certain variants 
of the hollow-core problem will be found for 
which the approximate solution is not applicable, 
in which case the knowledge of the rigorous 
solution would prove valuable. 

The complete solution for the heat input to 
hollow cores in induction heating processes is 
presented in the following. First, a development 
of the equation is given, and next a numerical 
application is presented. 


2. SOLUTION OF THE DIFFERENTIAL EQUATION 


For the following analysis, the unit length of a 
hollow cylindrical metal core that is located 
coaxially within a coil is considered. An annulus 
of the metal of radius r and thickness dr is now 
considered. If the magnetic field parallel to the 
axis of the core is H (instantaneous value), then 
the total flux @ within the radius r is 


6-20 f plirdr, (1) 
0 


where yp, the magnetic permeability, drops at the 
inner radius of the core from its value for the 
metal to unity. The electromotive force E along 
the circumference of the annulus is (with ¢ = time) 


(2) 
The current density G along the circumference of 


1 James Greig, World Power 25, 246 (1936). 

2 T. F. Wall, Engineer 173, 209 (1942). 

3R. M. Baker, Trans. A.I.E.E. 63, 273 (1944). 

*A. G. Warren, Mathematics Applied to Electrical Engi- 
neering (Chapman & Hall, London, 1941). 


E=—0/0t. 
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the annulus is 
G=E/2zxpr, 


(3) 


where p=electrical resistivity. Finally, the gra- 
dient of the magnetic field intensity H is 


dH /ar = —4nG. (4) 


These four equations contain four variables: 
H, ¢, E, and G as functions of r and t. Elimination 
of three of the four permits the setting up of the 
differential equation for any one of the four. 
For H, one has: 

8H 16H 4eu 6H 


—+- — =—_ —.. 


; (5) 
ér*> or “Or p 6t 


Since the alternating current in the coil is 
sinusoidal, H can be replaced by He?**/‘, where H 
now stands for the complex amplitude and 
f=frequency of the alternating current. Simi- 
larly, 0H/dt can be replaced by 2mifHe?**!'. 
Equation (5) then takes the form 


@H idil 
—+- —-—-1k*H=0, (6) 

dr? +r dr 
where . ‘ 
k?=8ryf/p. (7) 


In the following H stands for r.m.s.-value 
instead of amplitude; the form of Eq. (6) is 
obviously not affected by this change. The solu- 
tion of Eq. (6) can be found from the solution 
of the differential equation of the general form: 

d*y idy 





— +me"y =0 
dx*- x dx 


(m and n=constants), whose solution is: 
¥y =Z(miei"!*x), 


-where Zo is the sum of two Bessel functions of 

the first and third kind, both of zero order and 
each multiplied by an arbitrary complex con- 
stant. Substituting —k? for m and 2/2 for n, the 
solution of Eq. (6) becomes 


H=Zo(ibkir), (8) 


in which Zo, as stated before, is the sum of Jo 
and H,™, Bessel functions of the first and third 
kind, each multiplied by a constant. It is con- 
venient to introduce the so-called modified Bessel 


5 E. Jahnke and F. Emde, Tables of Functions (Dover 
Publications, New York, 1943). 
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functions according to the definitions 


Io(x) = Jo(ix) 


and 


cal 
Ko(x) = — (ix), 


so that Eq. (8) becomes 
H=A1,(ikr)+BKo(iikr), (9) 


where A and B are the complex constants of 
integration. 

By using Eq. (9), the other essential variables, 
namely the current density G and the electro- 
motive force E can be computed. Since both are 
sinusoidal and have the same frequency f as H, 
the form of Eqs. (3) and (4) is not changed if G 
and E from now on stand for the complex 
r.m.s. values, instead of instantaneous ones. 
From Eq. (4): 


URA VRB 


G=- Io (kr) -——Ky'(t}kr), (10) 
dn 





4r 


where the prime indicates differentiation after 
the argument. Further, from Eq. (3): 


iMkrpA BkrpB 
E-= He kr) —— Ke’ (@er). (11) 


Eqs. (9)-(11) give the essential variables as 
functions of the radius r within the metal part 
of the hollow core. 


3. COMPUTATION OF THE CONSTANTS 
OF INTEGRATION 


The next step is the evaluation of the con- 
stants of integration A and B from the boundary 
conditions. The two boundary conditions are: 

(1) For r=a (a=radius of the outer surface of 
the hollow cylinder), H assumes the value Ho 
which is the field inside the coil that is unaffected 
by the metal core. 

(2) For radii from zero to 6 (b=inner radius 
of the hollow cylinder) the field is constant and 
equal to Hy, the value that obtains from Eq. (9) 
for r=b. The first condition yields the equation: 


Hy) =AlI)(ttka) + BK o(itka). (12) 


The second condition can be expressed in the 
following manner. First, Hy is written: 


H,=AI)(i4kb) + BK o(i*kb). (13) 
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The total flux within the air space, ¢, is then 
b’rH,. Further : 0¢,/dt =i22°fb’H,. But this must 
be equal to — E,, where Ey, the e.m.f. at r=), is 
given by Eq. (11), if 7 is replaced by b. Hence, 
the second boundary condition yields the equa- 
tion: 


k 
A [1a — (V2) 2n'poru(ie)| 


k 
+B| Keb) . (Vi) 2a'fOK ib) =0. (14) 


The constants A and B can be evaluated in 
terms of Ho, a, and b from Eqs. (12) and (14). 
At this stage of the calculation, it is convenient 
to replace the Bessel functions by their real and 
imaginary components. The following notations 
are customary : 

I (tix) = ber x +i bei x 
and 
Ko(t!x) =ker x+i kei x. 


The corresponding differential forms are: 
Io! (t'x) =i- ber’ x+7 bet’ x 
Ko’ (i'x) =i-4 ker’ x+i} kei’ x. 


In order to obtain A and B in convenient 
forms, the following notations are introduced : 


u.= ber ka, 
v= bei ka, 
w=ker ka, 
z= kei ka, 


p=2 bei’ kb—kb ber kb, (15) 


q= —2 ber’ kb—kb bei kb, 
s=2 ket’ kb—kb ker kb, 
t= —2 ker’ kb—kb kei kb. 


Evaluation of the constants A and B can now 
be carried out. The result is: 





(s+1t) 
A =Hy>——_——_ — (16) 
(u+1v)(s+it) —(p+ig)(w+iz) 
anc 
—(p+iq) 
B=H,———_- ors ; (17) 








(u-+iv)(s+it) —(p+ig)(w+is) 
Both complex constants can be written as: 
A=A,+iA; (18) 
B=B,+1B;, (19) 


and 
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where the four real constants are given by Eqs. 
(16) and (17); numerical computation of these 
four constants is carried out according to stand- 
ard methods. 

It is easy to verify that for b=0, the solid 
core, B, becomes zero. The result for H in Eq. (9) 
is then given by the J» function alone.' 

It may be added that the present treatment 
holds not only for the case of a hollow metallic 
cylinder, but also for the general case of a solid 
core surrounded by a solid cylinder of another 
material. ¢ is then not simply b’xH,, but a more 
involved function of 6 and the constants of the 
solid core, as known from the equations valid for 
the solid cylinder. Since ¢@» is still linear in Hb, 
Eq. (14) will have the same general form, and 
only the second terms in the expressions (15) 
will be correspondingly more involved. 


4. CALCULATION OF THE HEAT INPUT 


The heat input into the hollow cylinder is 
calculated from the current density—Eq. (10)— 
in conjunction with Eqs. (18) and (19). The heat 
input per unit length W is given by its value 
_G|*p for unit volume, integrated over the cross- 
sectional area of the cylinder: 


W=2r0 { |G | *rdr. (20) 
b 


In order to solve this integral, it may be con- 
veniently written in the form 


W= ano [ GGrdr, (21) 
b 


in which G is the conjugate complex value of G. 
This integral can be solved by making use of 
the following general integral relation in which 
- Y, and Z, denote two different Bessel functions 
of pth order, both the sum of two functions of 
the first and third kind: 


f Y,(ax)Z,(Bx)xdx 





Bx Y,(ax)Zp_1(Bx) —ax Y,_1(ax)Zp(Bx) 
—~ a’? — 8? . 
Here a anc ° are constants. 
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Y ,(ax) is now identified with G, Z,(8x) with G, 
and x with r of Eq. (21). The four functions 
appearing on the right-hand side of the above 
equation, as well as a and 8, can then be easily 
evaluated, as follows. 

Concerning Y,(ax), it is given by Eq. (10), or, 
recalling the general relation 


Zo (x) = —Z (x) 
kA URB 
Y,(ax) = es +——K ,(i'kr). 


T T 


by 


(10a) 


From this it follows that p=1 and a=iik. 

Expressing Y,(ax) in the form of tabulated 
functions, Eq. (10) is converted, by using the 
form of Io’ as given in Section 3, into: 


kA 
Y,(ax) = —— kr+ bei’ kr) 
T 


kB 
——(ker’ kr+i kei’ kr). (22) 
4a 


Concerning Z,(8x), or G, it should be recalled 
that the conjugate value of any function is given 
by replacing all 2’s occurring in it by —7. Hence, 
from (10a) it follows that B=(—7)!k and from 
(22): 


kA 
Z (Bx) = ———(ber’ kr —i bei’ kr) 


us 


kB 
——(ker’' kr—i kei’ kr), (23) 
_ An 


where A and B are given by Eqs. (18) and (19) 
with —17 instead of i. 

As to Yp-1(ax), Eq. (10a) is used in replacing 
I, by Ip and K, by Ko. It then follows that 


kA 
Y,~1(ax) =——(i! ber kr — (—i)} bei kr) 


T 


kB 
+—(i! ker kr —(—i)! kei kr) 


T 


(24) 
and _ 
kA 
Zy-1(Bx) = 7 —i)* ber kr —i' bei kr | 
Tv 


kB 
+—[(—i)! ker kr—i' kei kr]. (25) 


us 
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Equations (22) to (25) represent the four 
functions to be substituted into the numerator 
of the right-hand side of the general integral 
given above. The denominator is given by 


a? — 6? = 2k? 


In carrying out the substitution indicated, the 
four constants A,, A;, B,, and B; are found to 





occur in four combinations as follows: 


Af+AP?=n, 

B/f+B;7= i, 
(26) 

A,B,+A.Bi=un, 

A,B;-—A ;B,= V. 


With these quantities, the heat input W becomes: 


p a 
W=— Rr (ber ber’ +-bei bet’) +-krd(ker ker’ +-kei kei’) +kru(ber ker’ +-bei kei’ 


Sr >| 


+ ber’ ker+ bei’ kei) +krv(bei ker’ — ber kei’+-bei’ ker—ber’ kei)| (27) 


where, for sake of simplicity, the argument kr 
after each function is omitted. 

It can be seen that, for the solid core for which 
B=0, only the first of the four terms within the 
bracket remains, giving the known solution for 
this case. 

5. NUMERICAL CALCULATION OF A SPECIAL 

CASE OF HOLLOW CORE 

In order to show the use of Eq. (27), a numeri- 
cal example treated by R. M. Baker* will be 
discussed. That author calculated the heat input 
by means of two approximate equations and a 
check of these, as explained in Section 1 of the 
present paper, can now be carried out by using 
the rigorous solution as given. It will be shown in 
the course of this section that the approxi- 
mate equations are in good agreement with the 
rigorous one. 

For low frequencies, characterized by the condi- 
tion that the depth of current penetration is larger 
than the wall thickness, Baker uses the formula 


61.8X 10-"5f2a°t,.pH,? 
V = ’ 
p?+15.5 K 10—'*f%at,,? 


where W, the heat input per unit length, is given 
in watts/cm and ¢,, the wall thickness, is meas- 
ured in cm, p is measured in ohms cm, and 
u=1 (non-magnetic metal). 

For high frequencies, characterized by the 
condition that the depth of current penetration 
is smaller than the wall.thickness, the formula 


W =7.92 X 10-‘a(fp)'Ho? 
is used. At the frequency where the two calcu- 


lated curves cross they are smoothed into one. 
In order to check the two approximate equa- 
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tions, one of the cases calculated by Baker was 
recalculated by means of the rigorous Eq. (27) 
developed in Section 4 and valid for the whole 
frequency range. This special case refers to a non- 
magnetic metal (u=1), an outer diameter of the 
cylinder of 10 cm (a=5), and a wall thickness of 
1 millimeter (6=4.9). p is taken as 7X 10~* ohm 
cm (=7X10* e.m.u.). Since the heat input is 
proportional to H,? according to all the equations 
used, W/H,? is calculated. A frequency range 
from 300 to 50,000 is covered. 

Before giving the results of the complete calcu- 
lation, the computation will be shown in detail 
for one particular frequency. A frequency very 
close to the crossing point of the two approximate 
curves is selected, i.e., 17,900. From Eq. (7) we 
then have k=14.2, which gives ka=71.20 and 
kb=69.78. 

The first step is the evaluation of the sixteen 
Bessel functions involved for ka and kb as argu- 
ments. For arguments up to 10, the functions 
have been tabulated previously.5 They are some- 
times found under the real and imaginary parts 
(Re and Im) of the respective J and H™ func- 
tions. It is useful to remember ‘the following 
identities : " 


‘ o 
ber x = Re Jo(t'x), kerx= ‘‘: Im Hy (ix) ; 


T 
bei x = —Im Jy(i'x), kei x= 3 Re Ho (ax) ; 


us 
ber’ x= —Ret!J,(i'x), ker’ x=— Im Hy (ix) ; 
Po 


TT 
bet’ x =Im iJ ,(i'x), kei’ —. Re tH, (ax). 
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For arguments larger than 10, tables do not 
exist, but there are asymptotic formulas that 
permit the computation. The eight formulas are 
listed below. Let 


F=e7!*/(29x)! F’ =(x/(2x))'e7-2/” 


and 
g=x/(2)!+x/8 p=x/(2)!—x/8. 
Then: 


ber x = F{.(1+0.088/x) cos 
+ (0.088 /x+0.070/x") sin y ], 


bei x = F((1+0.088/x) sin y 
— (0.088 /x+-0.070/x*) cos y |, 


ber’ x = F((1—0.265/x) cos ¢ 
— (0.265/x+0.117/x*) sin ¢ ], 


bei’ x = F[(1—0.265/x) sin ¢ 
— (0.265 /x+0.117/x?) cos ¢], 


ker x = F’[ (1—0.088/x) cos ¢ 
— (—0.088/x+0.070/x*) sin ¢ ], 


kei x = F’[ —(1—0.088/x) sin ¢ 
— (—0.088/x+0.070/x?) cos ¢ ], 


ker’ x= FL —(1+0.265/x) cos p 
+ (0.265 /x—0.117/x*) sin y], 


kei’ x = F’[ (1+0.265/x) sin y 
+ (0.265 /x—0.117/x?) cos yp}. 


Using these formulas, the following data are 
obtained for the particular frequency mentioned 
above: 


ka=71.20 

ker = 0.01800 X 107?! 

ket = —0.00920 x 10-*! 
ker’ = —0.01936 x 10-! 
kei’ = —0.006165 x 10-*! 


kb =69.78 
ber = 0.03220 X 10! ker = 0.04800 X 10-*! 


bei = —0.1246 x 10”! kei = 0.02828 x 10-*! 
ber’ = 0.1106 x 107! ker’ = —0.01429 x 10-*! 
bei’ = —0.06443 X10"! kei’ = —0.05414 K 10-'. 


ber = 0.3033 X 10”! 

bet = —0.1070 X 10?! 
ber’ = 0.3068 X 10”! 
bei’ = 0.1587 < 10?! 


The next step is the evaluation of the eight 
quantities given by Eq. (15), by means of which 
A and B from Eqs. (16) and (17) are calculated. 
The four components shown in Eqs. (18) and (19) 
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Fic. 1. Heat input as function of the frequency for a 
special case of a hollow cylinder. Solid line calculated from 
Eq. (27); dotted line calculated by R. M. Baker from 
approximate formulas. 


are then in the present case: 
A,=2.823 X10-"' Ao, 
A; =0.8252 X10-"' Ay, 
B,= — 1.609 X 107'FZo, 
B;=0.8668 X 10777». 


These data permit evaluation of the four 
quantities shown in Eq. (26), and in this manner 
all magnitudes occurring in Eq. (27) are known. 
In the present case the result is 


W = 1.43 X10“? 


in watts/cm, if H» is measured in e.m.u. 

The calculation was carried out for five selected 
frequencies and the result of the complete calcu- 
lation, representing the heat input as a function 
of frequency, is shown in Fig. 1. The solid line is 
calculated from Eq. (27), the dotted line is 
Baker’s curve. As can be seen, the two methods 
are strictly identical below the frequency for 
which the depth of penetration is equal to the 
wall thickness ; for higher frequencies the approxi- 
mate formula gives values slightly higher than 
those from the rigorous equation. For still higher 
frequencies the two curves merge again; it can be 
shown that for depths of penetration markedly 
less than the wall thickness Eq. (27) degenerates 
to the second of the approximate formulas given. 
The highest frequency for which computation 
was made is 4000 kc (not shown in Fig. 1); W, 
according to both equations, is 20.9 watts/cm. 
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The Tensile Properties of Several Special Steels and Certain Other 
Materials under Pressure 


P. W. BRIDGMAN 
The Physics Laboratories, Harvard University, Cambridge, Massachusetts 


(Received November 16, 1945) 


INTRODUCTION 


N a series of papers! the effect of hydrostatic 

pressure on the tensile properties of a large 
number of different steels has been investigated. 
The pressure effects on all these steels fell ap- 
proximately into a one-dimensional pattern, the 
behavior being determined by a single parameter, 
the hardness or the absolute magnitude of the 
yield point at atmospheric pressure. In spite of 
the great ostensible variety of the steels, they 
were, nevertheless, restricted to a range of only a 
few percent in any of the conventional alloying 
elements, and to a hardness range below 50 
Rockwell C. In the present investigation two 
different series of steels were used falling outside 
the previous limits. The first of these comprised 
four stainless steels running up to 18 percent Cr 
and 10 percent Ni; the second was a series of 
five different hardnesses of a carbon steel, 
running from 48 to 63 Rockwell C. 

In anticipation of the results, it was found 
that two of the stainless steels fall markedly 
outside the previous pattern, but the others fell 
approximately within it. 

In addition to the measurements on the steels, 
qualitative measurements were made on several 
other materials. 

In the discussion, the new results of this paper 
and previous unclassified results are surveyed 
together. 


APPARATUS AND METHODS 


The apparatus and method for making the 
tensile tests have been fully described in the 
previous communications; only a brief charac- 
terization will be made here. The tensile speci- 
mens, of an over-all length of 1 inch and a 


1 Bight Reports to the Watertown Arsenal by P. W. 
Bridgman, from March, 1943 to December 21, 1944, 
numbered No. WAL 111/7, No. WAL 111/7- 1, etc., 
through No. WAL 111/7-7. See also P. W. Bridgman, Rev. 
Mod. Phys. 17, 3-14 (1945). There are also several reports 
by P. W. Bridgman to the N.D.R.C., on which the con- 
fidential classification had not yet been lifted at the time 
of writing. 
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diameter varying about a mean of 0.1 inch, de- 
pending on the absolute strength, are completely 
surrounded by a fluid in which pressures up to 
30,000 kg/cm? may be produced and maintained 
during stretching. The specimen is stretched by 
a forked yoke driven by the piston with which 
pressure is produced. The tension on the speci- 
men is measured within the pressure chamber by 
the change of resistance of a “grid,” in this way 
avoiding errors from the friction of the piston. 
Simultaneous measurements were made of hydro- 
static pressure, tension on the specimen, and 
elongation of the specimen, for progressively 
increasing elongations. The elongation can be 
determined with adequate accuracy only in the 
range of plastic flow, after reaching the maximum 
load. The parameters singled out for determina- 
tion from the complete data and presented in 
the following are: maximum load, load at 
fracture, reduction of area at fracture, and 
character of the fracture—all as functions of the 
pressure. In many experiments stretching was 
not pushed to fracture; for these the final loads 
were determined instead of the load at fracture. 

In calculating the stresses from the data the 
change of geometrical configuration due to 
necking is so large that it was necessary to make 
allowance for the distribution of stress across the 
neck. This was done by methods which have been 
outlined in several papers.” The stress system at 
the neck of a tension specimen can be analyzed 
into the sum of three stress systems. The first 
of these is the hydrostatic pressure in the ambient 
liquid ; the second is a tension along the axis of 
the specimen, constant across the section, called 
the ‘“‘flow’’ stress in the following; the third is a 
variable hydrostatic tension (three equal com- 
ponents of principal stress) varying from zero at 
the outer surface of the neck to a maximum on 
the axis according to a formula given in detail in 

2 P. W. Bridgman, Trans. A.S.M.E. 32, 553-574 oe 
Metals Technology, Technical Publication No. 1792 of 


Am. Inst. Min. Met. Eng. (Dec. 1944) with supplement in 
April, 1945; Wat. Ars. Report, No. WAL 111/7-4. 
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the references. It follows that the ‘‘flow’”’ stress 
is not equal to the average tension across the 
section (total tensile load divided by the actual 
area of the neck), but is less than it. The factor 
by which it is less will be called the “correction 
factor.” In the following work the correction 
factor varied through an extreme range from 1.0 
to 0.7. 


THE MEASUREMENTS 
1. Stainless Steels 


Measurements were made on four different 
varieties of stainless steel. For these I am in- 
debted to Dr. A. L. Feild of the Rustless Iron 
and Steel Corporation. The designations, com- 
positions, and heat treatments are given in 
Table I. 

The material was supplied in }-inch rods. 
These had been drawn to approximately the final 
size, then annealed and pickled. After this they 
were given a final cold drawing pass with further 
reduction in cross-sectional area of 5 percent. 
The final hardnesses on the Rockwell B scale 
were 87, 93, 95, and 93, respectively, for H-26, 
H-50, H—29, and H-88. 

Seven runs were made on each of these four 
steels at pressures ranging from atmospheric to 
25,000 kg/cm’. In four of these runs, including 
the run at atmospheric pressure, the test was 
carried to fracture, and in the three others it 
was stopped short of fracture. The purpose of 
the latter runs was to obtain points on the strain 
hardening curve comparatively free from the 
error introduced by the disturbing local flow 
which is very likely to occur in the moment of 
fracture. 

These steels fall into two groups, of markedly 
different properties. At atmospheric pressure 
H-26 and H-50 break in the fashion normal for 
- a steel of high ductility, with a reduction of area 
of about 78 percent. Necking starts, as is normal, 
at a strain of about 0.2. Under pressure these two 
steels show an increase of ductility so large that 

















TABLE I, 
3 Mn P Ss Si Cr Ni Mo Cu 
H-26 0.094 0.36 0.023 0.022 0.35 12.26 0.46 0.50 ND 
H-50 0.067 0.50 0.020 0.030 0.51 17.49 041 ND ND 
H-29 0.058 0.70 0.030 0.013 0.85 18.51 895 ND 0.20 
H-88 0.051 0.59 0.030 0.020 0.47 18.31 10.27 ND 0.20 
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Fic. 1. The flow stress as a function of strain (natural) 
for four stainless steels. An f on a number of points indi- 
cates the occurrence of fracture. 


the maximum pressure was not pushed beyond 
19,000 kg/cm’, the reduction of area at this 
pressure being so extreme that error is introduced 
by deviations from geometrical regularity be- 
cause of finite grain size. The two other steels, 
H-29 and H-88, break with the same reduction 
of area at atmospheric pressure, but the character 
of the stress-strain curve is quite different, since 
the elongation is much greater and necking does 
not start until a strain of 0.65 is reached. Under 
pressure the increase of ductility is not nearly 
as large as for the two other steels, and pressure 
could be raised to 25,000 without enough re- 
duction of area to introduce error from geo- 
metrical irregularity. 

Because of the much greater reduction of area 
before necking starts, the contour of a tension 
specimen of H-29 or H-88 has a different shape 
from the other two or from the other steels inves- 
tigated previously. Since the shape of the con- 
tour determines the correction for stress dis- 
tribution across the neck, the shape of the 
contour, that is, the radius of curvature of the 
contour at the neck, was measured for all speci- 
mens of all four steels. The curve of a/R (radius 
of curvature of the contour at the neck divided 
by the radius of the cross section at the neck) 
against strain was normal for H-26 and H-S0, 
but for H-29 and H-88 it lies below the normal 
curve in the initial stages, cutting the axis at a 
strain of 0.65 instead of at 0.2, and approaches 
coincidence with the normal curve at a strain 
of 3.5. The corrections for the stress distribution 
in H-29 and H-88 were determined by utilizing 
these special measurements of a/R. 

The results are reproduced graphically. In 
Fig. 1 is shown the flow stress in terms of the 
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STRAIN AT FRACTURE 


Fic. 2. The hydrostatic pressure at fracture as a function 
of the strain at fracture for four stainless steels. 


strain; the concealed variable is the hydrostatic 
pressure, which in general is greater for the larger 
strains. An f written against some of the points 
indicates that the specimen was fractured during 
the test. In general the scattering of the points 
is somewhat greater for the fractures; this is due 
to distortion in the process of fracture, which 
may be considerable at the largest strains. Within 
experimental error the relation between flow 
stress and strain is linear in the present range. 
More accurate measurements under especially 
favorable conditions have shown that the curves 
are slightly concave toward the axis of strain. 

In Fig. 2 is shown the pressure at fracture 
against the strain at fracture. In accordance with 
the great majority of previous results, the rela- 
tion is assumed to be linear, and is so shown. 
The strain at fracture can be determined with 
some accuracy, however, and instances have 
already been found in which the relation is 
almost certainly not linear within experimental 
error. It is therefore not at all impossible that the 
relation for H-88 should be represented by a 
curve concave toward the pressure axis. 

In Fig. 3 is shown the flow stress at the 
maximum load (that is, at the point where 
necking starts) as a function of the pressure of 
pulling. For all four steels this flow stress 
increases linearly with pressure, by some 3 
percent per 10,000 kg/cm? for H-26 and H-50, 
and by some 13 percent for the same pressure 
increment for H-29 and H-88. The former rate, 
3 percent, is somewhat less than for most steels 
previously investigated ; the latter rate, 13 per- 
cent, is distinctly higher. Except for a factor 
constant with respect to pressure, the ‘‘flow 
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stress at maximum load” is the same as the 
technical ‘‘tensile strength.’”’” The factor to 
reduce from flow stress to tensile strength is 
0.83 for the first two steels and 0.52 for the 
second two; this factor is determined by the 
strain at which necking starts, the ‘tensile 
strength” being calculated on the original cross 
section. 

The fracture of all four stainless steels roughly ° 
follows the same pattern previously found. At 
atmospheric pressure the fracture is of the con- 
ventional “‘cup-cone”’ type. With increasing pres- 
sure the tensile part of the fracture (bottom of 
the cup) becomes progressively less in extent 
relative to the total cross section, and eventually 
disappears, the fracture becoming all shear above 
a certain pressure. In the lower part of the 
domain of completely shearing fracture the sur- 
face of shear is conical; at higher pressures this 
multiplies into a double cone, and at still higher 
pressures the shearing fracture shifts to a single 
distorted plane. These features are all present in 
the stainless steels. In Fig. 4 is shown the ratio 
of the area of the tensile part of the fracture to 
the total area of the fracture as a function of the 
hydrostatic pressure prevailing at the moment 
of fracture. The fracture becomes all shear at a 
lower pressure for the two steels H-26 and H-50 
than for the two others, H-29 and H-88. The 
intermediate episode of shearing fracture on 
single or multiple conical surfaces is confined to 
a comparatively narrow pressure range, and was 
definitely observed for only one specimen, most 
of the shearing fractures being on a single dis- 
torted plane. The general character of the 
fracture was always more irregular and coarsely 
granular than for ordinary carbon steels. 
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Fic. 3. The flow stress at the maximum load (i.e., at the 
beginning of ce meng as a function of the prevailing 
hydrostatic pressure for four stainless steels. 
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2. Steels of Hardness Greater than Rockwell 
C 48 


I am greatly indebted to Captain John H. 
Hollomon, who had the specimens prepared at 
the Watertown Arsenal under his direction. The 
specimens were all made from the same bar of 
high carbon steel. They were heat-treated by 
holding at 1525°F for 30 minutes and then 
quenching into oil at room temperature. They 
were then drawn in batches of ten to five dif- 
ferent temperatures so chosen as to give median 
hardnesses of 48, 51, 56, 60, and 63. 

These experiments followed a number of 
unsuccessful attempts to demonstrate an effect 
of hydrostatic pressure on the ductility of tool 
steels of hardnesses greater than 60. The pre- 
liminary attempts had encountered brittle failure 
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Fic. 4. The ratio of the area of the tensile part of the 
break (bottom of the cup) to the total area as a function 
of the pressure at which the break occurs for the four 
stainless steels. The two groups H-26, H-50 and H-29, H-88 
are so nearly equal that a single line is drawn to represent 
the average behavior of the two. 
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at the shoulders of the test specimen because of 
stress concentration. After a number of trials the 
shape of the specimen shown in Fig. 5 was 
arrived at. The specimens were rough finished 
at the Arsenal to the shape shown and heat- 
treated. They were finished to final dimensions 
at Harvard by grinding; | am greatly indebted 
to the skill of Mr. Charles Chase for this opera- 
tion. The minimum diameter of the finished 
specimen was varied somewhat with the absolute 
strength of the steel; it was so chosen as to 
demand approximately the same maximum load 
for all heat treatments. The hardness of each 
specimen was determined on the flats ground on 
the ends as indicated. From each batch of ten, 
the four or five were selected whose hardness was 
closest to the median. Since the primary object 
of these measurements was to determine how 
much ductility was imparted by pressure, these 
specimens were all pulled to fracture at pressures 
ranging from atmospheric to a maximum of 
28,000 kg/cm*. The strain at fracture at at- 
mospheric pressure as a function of hardness is 
shown in Fig. 6. 

The same sort of measurements were made on 
these as on the stainless steels. 

In Fig. 7 is shown the flow stress at fracture 
as a function of the strain at fracture. The error 
from distortion in the act of fracture is much 
less for these hard steels than for the softer 
steels, and Fig. 7 therefore gives essentially the 
strain hardening curves. The points are some- 
what more scattered than for the softer speci- 
mens. Again the relation is linear within the 
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Fic. 5. The special shape of the tensile specimens of the hard steels. Diameters given provide for finish grinding. 
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Fic. 6. The strain at fracture at atmospheric pressure of 
the hard steels as a function of the Rockwell C hardness. 


limits of error. Two atmospheric points for 
hardness 63 lay distinctly out of line with the 
others and are not plotted; it is probable that 
there was some surviving stress concentration 
effect here. Only a single experiment under 
pressure was made on the 63 group, at 24,000 
kg/cm*. A slight but measurable ductility was 
found, the strain at fracture being 0.06. In view 
of the very small ductility at this pressure it did 
not seem worth while to try for measurements 
on this hardness at lower pressures. 

In Fig. 8 is shown the hydrostatic pressure at 
fracture as a function of the strain at fracture. 
Again the relation seems to be linear within 
experimental error. The lines drawn in Fig. 8 
are the best lines through the points of each 
group. Doubtless if a larger number of measure- 
ments had been made, lines would have been 
found more uniformly spaced in terms of the 
hardness. A vertical line in Fig. 8 means zero 
effect of pressure on ductility. The rapid decrease 
of pressure effect on ductility with increasing 
hardness in this hardness range is noteworthy. 
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Fic. 7. The fiow stress at fracture as a function of the 
strain at fracture for the hard steels. The figures on the 
lines indicate the average Rockwell C hardness of the 
group. 
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In Fig. 9 is shown the flow stress at maximum 
load as a function of the hydrostatic pressure at 
fracture. Since the hydrostatic pressure at 
fracture was only slightly higher than the hydro- 
static pressure at the time of maximum load, 
this is practically a graph of flow stress at 
maximum load against the prevailing hydro- 
static pressure. The factor for reducing from 
flow stress to tensile strength is 0.82 for all the 
specimens of Fig. 9. As in the case of Fig. 3, the 
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Fic. 8. The hydrostatic pressure at fracture as a func- 
tion of the strain at fracture for the hard steels. The 
numbers on the lines indicate the average Rockwell C 
hardness of the group. 


flow stress at maximum load increases somewhat 
with increasing pressure. The increase is greater 
on a percentage basis for the harder heat treat- 
ments. 

In Fig. 10 is shown the ratio of the area of the 
tensile part of the fracture to the total fractured 
area as a function of the pressure. At atmospheric 
pressure the tensile part of the fracture occupies 
a progressively greater part of the total area as 
hardness increases, the fracture of the steel of 
hardness 63 being practically all on a single fine 
grained tensile plane except for a very narrow 
rim of shear around the edge. As usual, the tensile 
break disappears at high pressure. The pressure 
at which it disappears is practically the same for 
all the hardness groups except the hardest, which 
differs abruptly from the others. The multiple 
shearing fracture does not often appear in the 
harder steels, but at the lower end of the hardness 
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range is more common. The most common type 
of shearing break is on a single distorted plane; 
this is always highly burnished. In a couple of 
instances there were also fissures on longitudinal 
planes. 


3. Tensile Properties under Pressure of 
Other Materials 


In addition to the tests on the two series of 
steels just described, other sorts of material were 
examined, chiefly from a qualitative point of 
view and primarily to determine whether duc- 
tility is imparted to all sorts of material by 
pressure. 

A number of experiments were made on gray 
cast iron. This was the ordinary gray cast iron 
of commerce, fine grained and without obvious 
imperfections. Several experiments were first 
made in which the specimen was given the con- 
ventional shape of most of the previous speci- 
mens. These tests all gave negative results, the 
specimen fracturing under pressure with a clean 
tensile break and no evidence of plastic flow. 
Two specimens were then tried with the special 
shape shown in Fig. 5 to avoid any effects of 
stress concentration. These were both pulled at 
29,000 kg/cm*. One was exposed directly to the 
action of the pressure transmitting liquid, tso- 
pentane, whereas the other was copper plated in 
order to prevent penetration by the liquid into 
any pores that there might be. Both experiments 
gave the same completely negative result, clean 
fracture perpendicular to the axis of the speci- 
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Fic. 9. The “average longitudinal stress,”’ which in this 
case is the same as the flow stress, at the maximum load 
as a function of the hydrostatic pressure for the hard 
steels. The numbers on the lines indicate the average 
Rockwell C hardness of the group. 
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Fic. 10. The ratio of the area of the tensile part of the 
break (bottom of the cup) to the total area as a function of 
the pressure under which the break occurs for the hard 
steels. The numbers on the lines indicate the hardness 
group. 


men, with no measurable reduction of diameter. 
It is perhaps natural to ascribe the negative 
results to the microscopic grains of pure graphite 
known to be present in this material. 

For unity of presentation, certain results with 
other materials which have already been de- 
scribed elsewhere will be referred to again here.* 
Tension tests under pressure on specimens of the 
conventional shape were made on aluminum, 
copper, bronze, and brass. The ductility of all 
these is markedly increased by pressure. In 
general it is not possible to push the strains to 
such extreme values without loss of geometrical 
regularity with these materials as with the steels. 
The cross section is likely to become dominated 
by a single crystal grain, so that the section 
becomes rectangular or square, with the result 
that the stress at the section becomes inhomo- 
geneous and the determination of the stress- 
strain relation is thwarted. No serious attempt 
was made to find whether this effect was due to 
abnormally large grains in the virgin material or 
whether large grains were developed by the proc- 
ess of straining under pressure; my general im- 
pression in the case of aluminum and copper was 
that it was the latter. 

Qualitative measurements were also made on 
a number of other materials that could not be 
readily fashioned into the conventional tensile 
specimens. These materials were made into 
slender rods, approximately 1 cm long and 0.8 
mm in diameter, supported near the two ends 


3 P. W. Bridgman, Rev. Mod. Phys. 17, 3-14 (1945). 
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and subjected to a push at the center by a blunt 
wedge. Stops so limited the travel of the wedge 
that ductile materials were bent to a maximum 
angle of 135°. The experiment consisted in ex- 
posing the rod to the bending action of the wedge 
under hydrostatic pressure and observing whether 
there was any measurable difference between the 
amount of bending in the neighborhood of the 
fracture (if there was a fracture) as compared 
with atmospheric pressure. The apparatus was 
first tried with rods of a 1.50 C tool steel, 
hardened glass hard and drawn to various tem- 
peratures. At 29,000 kg/cm? the undrawn speci- 
men broke with no appreciable bending, but all 
the others received the whole 135° bend without 
fracture. Of these the specimen drawn at 200°C 
broke at atmospheric pressure with no appreci- 
able bending and the one drawn at 275° broke 
after bending through a slight angle. Repetition 
of the same experiment at 25,000 kg/cm? showed 
a rather rapid falling off of the pressure effect, 
the specimen tempered at 200° now breaking 
after bending only 40°. 

At 29,000 gray cast iron broke with no appre- 
ciable bending, thus checking the result by the 
other method. 

Other materials tested by the same method 
were: antimony, tellurium, bismuth (all in the 
single crystal condition), borax glass, hard 
rubber, and limestone. A rod of bismuth was 
recovered after an exposure at 15,000 per- 
manently bent without fracture to an angle of 
80°; whether the difference between 80° and 
135° was due to elastic recovery was not inves- 
tigated. Antimony, arsenic, and tellurium might 
be expected to show some ductility in analogy 
with bismuth, but loading of these three ma- 
terials at 29,000 resulted in clean transverse 
fractures with no trace of any bending before 
the occurrence of fracture. Hard rubber was 
recovered after exposure to similar conditions 
unfractured and with a set of 10°; this set 
entirely disappeared and the rod became straight 
again on standing 12 hours at room temperature 
at atmospheric pressure. The elastic limit of 
hard rubber is thus obviously markedly raised 
by pressure. The borax glass was from the same 
source as that which Birch and Dane‘ had used 


4E. B. Dane and F. Birch, J. App. Phys. 9, 669-674 
(1938). 
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in measuring the effect of pressure on viscosity 
at elevated temperatures. This glass deforms 
at atmospheric pressure with measurable ve- 
locities at temperatures as low as 300°; it might 
therefore be expected to show some measureable 
deformability under pressure at room tempera- 
ture. No evidence of any bending before fracture 
was found, however, at a pressure of 29,000. 
Doubtless time effects are of importance here; 
the load was applied at such a rate as to reach 
the maximum deformation in a time of the order 
of a minute. 

Limestone was of special interest because it 
had already been extensively investigated in the 
experiments of Griggs and of Balsley.5 They 
found marked plasticity under compressive 
stresses in combination with hydrostatic pres- 
sure and also in simple tension, but much less 
ductility in bending. The pressures that could be 
reached in the present experiments were higher 
than those attained by Griggs and by Balsley, 
and the external conditions were somewhat more 
uniform, so that a positive result might perhaps 
be expected. The fracture was, however, com- 
pletely brittle at pressures up to 29,000 kg/cm? 
when the test was performed in the regular way. 
It is known from the experiments of Birch and 
Bancroft® and others that the behavior of rocks 
under pressure is in many cases qualitatively 
altered if the pressure transmitting fluid is 
prevented from having access to the interior of 
the specimen by jacketing it with a suitable 
covering. The experiments were therefore re- 
peated, encasing the limestone in a heavy sheath 
of lead or copper or steel, and bending under 
29,000. The lead and copper jackets were prob- 
ably completely successful; the specimens were 
markedly bent, and there were no obvious flaws, 
but the tests were not entirely conclusive because 
the specimens were broken on attempting to re- 
move the jackets at atmospheric pressure. In the 
case of the steel jacket, both the jacket and the 
specimen broke brittlely with no obvious perma- 
nent deformation. The conclusion would there- 
fore seem to be that the limestone is capable of 
receiving permanent deformation in tension 


5 D. T. Griggs, J. Geol. 44, 544-571 (1936); Am. Mineral. 
23, 28-33 (1938); Bull. Geol. Soc. Am. 51, 1002-1022 
(1940). J. R. Balsley, Trans. Am. Geophys. Union, Part II, 
519 (1941). 


6 F. Birch and D. Bancroft, J. Geol. 46, 113-141 (1938). 
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without fracture at pressures as high as 29,000. 
The extension of this conclusion to ordinary 
limestone is complicated by the fact that there 
is a polymorphic transition in the neighborhood 
of 25,000 so that the plasticity demonstrated 
above is plasticity of another modification than 
that existing under normal conditions. 

In the endeavor to acquire further under- 
standing of the effect of pressure on the me- 
chanical properties of glass hard steel, another 
type of purely qualitative experiment was per- 
formed which has been previously made for 
other materials. If a thin ring of hard rubber or 
glass is slipped over a tightly fitting solid steel 
core and then exposed to hydrostatic pressure in 
a medium in which ring and core are both totally 
immersed, the ring will be split, just as if it had 
been expanded by driving a tupp into it. The 
simplest explanation is that the 
because of its much lower cubic compressibility, 
forces the ring to be stretched beyond the dimen- 
sions that it would normally assume under the 
hydrostatic pressure. This suggests that, for 
brittle materials, the normal zero from which to 
figure the extension which will produce brittle 
tensile fracture is the configuration which the 
material would assume under the mean equiva- 
lent hydrostatic stress of the total stress system. 

The corresponding experiment was made with 
a 1.50 C tool steel by fitting a thin ring of the 
steel in the glass hard condition over a Carboloy 
core. The two were ground together with a 
slight taper, and were pushed snugly together, 
so that the steel ring was under a slight initial 
tension. The ring with core was then exposed to 
a pressure of 29,000 kg/cm’. On releasing 
pressure, not only was the steel ring unbroken, 
but it was still a tight fit on the core, showing 
that it had received no appreciable plastic set. 
If all the elastic constants of Carboloy, as well 
as the cubic compressibility, are assumed to be 
about one-third the corresponding constants of 
steel, which has never been proved by direct 
experiment but which is not unplausible, it may 
be calculated that at the maximum pressure the 
steel ring was strained by about 2.5 X 10-* beyond 
the dimensions which it would normally have 
had under a pure hydrostatic pressure of 29,000. 
This stretch is far beyond the elastic limit or 
fracture point at atmospheric pressure, and 


steel core, 
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Fic. 11. The flow stress as a function of strain (strain- 
hardening curves) for all the steels of this paper and 
various 1045 steels previously described in Watertown 
Arsenal reports. 


would demand a tension of approximately 
750,000 p.s.i. So large an effect of hydrostatic 
pressure on the elastic limit has not been found 
previously for the softer steels. It may be, 
therefore, although this conclusion should be 
confirmed by other methods, that the com- 
paratively small effect of pressure on the duc- 
tility of very hard steels is compensated by a 
large effect on the elastic limit. 

Special attention was given to an investigation 
of the effect of pressure on the tensile properties 
of Carboloy. The interest in this lies in the fact 
that Carboloy containers, immersed in a medium 
carrying a high pressure, successfully withstand 
an internal pressure of more than 100,000 
kg/cm?.? The interiors of such containers are 
subject to large tensile strains and stresses, so 
that a large effect of pressure is indicated, Car- 
boloy under normal conditions being brittle and 
with comparatively little strength in tension. 
The method used was the bending of slender 
rods, as before. No measurable plasticity in 
bending was imparted to a rod of grade 999 
Carboloy by a pressure of 29,000 kg/cm?. A 
special grid of greater sensitivity was then con- 
structed and with this it was found that the 
bending force required to produce fracture at 


7P. W. Bridgman, J. App. Phys. 12, 461-469 (1941). 
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Fic. 12. Collected results for the steels of this paper and 
of a previous Arsenal report for the pressure at fracture as 
a function of the strain at fracture. 


27,000 was at least threefold, and perhaps con- 
siderably more, greater than at atmospheric 
pressure. This very notable increase of strength 
must be an important factor in the functioning 
of the high pressure containers. The behavior is 
qualitatively the same as indicated in the last 
paragraph for hard steel. 


DISCUSSION 


In Fig. 11 are collected various results of the 
previous and the present investigations for flow 
stress versus strain for different grades of steel 
(that is, the strain hardening curves). The five 
lines designated by A2, A5, A6, A7, and A8 are 
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Fic. 13. Various significant stresses as a function of the 
strain at fracture for stainless steel H-29. 
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taken from the Watertown Arsenal reports, and 
refer to various heat treatments of a 1045 steel. 
The two groups of similar stainless steels are 
each averaged to give a single line. The lines of 
Fig. 11 form approximately a fan shaped family 
of non-intersecting lines, the slope of the line 
being higher the higher the absolute value of the 
stress, with the exception of the stainless steels 
H-29 and H-88. For these two steels the slope 
of the strain hardening line is higher than would 
be expected from the tensile strength at atmos- 
pheric pressure. That is, these two steels harden 
at an abnormally rapid rate. 

Figure 12 shows a similar collection for dif- 
ferent steels of the curves for pressure at fracture 
as a function of the strain at fracture, that is, the 
effect of pressure on ductility. The two undesig- 
nated lines are for the Watertown steels, which 
fall into two groups with respect to the property 
shown in Fig. 12. Again all the steels of this 
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Fic. 14. Various significant stresses as a function of the 
strain at fracture for stainless steel H-50. 


paper and two of the stainless steels fall into the 
same family of curves as the steels previously 
investigated, the increase of ductility under a 
given increment of hydrostatic pressure being 
greater the lower the absolute tensile strength. 
Again the two stainless steels H-29 and H-88 are 
exceptions. Their increase of ductility under 
hydrostatic pressure is abnormally low. Thus for 
these two steels a high rate of strain hardening 
with increase of strain goes with a low increase of 
ductility with increasing pressure. 

The data of this paper bear on the question of 
the best criterion of rupture in the present range 
of conditions. In the fifth Arsenal report it was 
found that for the steels of those reports the 
mean hydrostatic tension (one-third the sum of 
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the three principal stresses, abbreviated in the 
following as S) was approximately constant at 
fracture. In Figs. 13-16 data are shown with 
respect to fracture for four representative steels 
of this paper, that is, stainless steels H-50 and 
H-29 and the groups of hardness 48 and 60. In 
these figures the lines of flow stress, F, and of 
hydrostatic pressure at the instant of fracture 
taken as positive, P, were first drawn as func- 
tions of strain. Next, from the line of flow stress 
was drawn the curve of “hydrostatic tension,” 
HT, on the axis of the specimen. This hydro- 
static tension is the second of the three com- 
ponents of the complete stress system described 
above. It was computed from the flow stress by 
multiplying by a somewhat complicated function 
of a/R which has been given in detail in the 
papers describing the corrections for stress dis- 
tribution at the neck. This function of a/R may 
be expressed as one definite function of strain 
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Fic. 15. Various significant stresses as a function of the 
strain at fracture for the specimens in the 48 hardness 
group. 


for all normal steels, and as another function of 
strain for the two abnormal stainless steels. In 
both cases the conversion from a function of a/R 
to a function of strain was made with the use of 
empirical data. The function converting flow 
stress into “‘hydrostatic tension” is highly non- 
linear in the strain. It is a first empirical result 
that in the case of all steels, normal and abnormal 
alike, the curve for HT versus strain obtained 
with the use of this highly non-linear factor is 
linear within limits of error. It is a second em- 
pirical result that in all cases the line for HT 
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Fic. 16. Various significant stresses as a function of the 
strain at fracture for the specimens in the 60 hardness 
group. 


rises approximately parallel to the line for F. 
Next, the two lines for F and HT were combined 
to give the line for F+HT. This is the axial stress 
at fracture that would be found on the axis of 
the specimen if there were no _ hydrostatic 
pressure. Inspection of the figures shows that in 
all cases it increases rapidly with strain. If the 
hydrostatic pressure is subtracted from F+HT 
the actual stress Z, on the axis at the instant of 
fracture will be obtained. This line is also shown 
in the figures. Finally, the line for the mean 
hydrostatic tension (S=}F+HT-—P) was con- 
structed. 

This method of representing the results is 
somewhat different from any that I have em- 
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Fic. 17. Various significant stresses as a function of the 
strain at fracture for the steel A-5 of the fifth Arsenal 
report. 
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Fic. 18. Various significant stresses as a function of the 
strain at fracture for the steel A-7 of the fifth Arsenal 
report. 


ployed before. Two of the steels previously 
examined in the Arsenal reports are therefore 
analyzed in the same way, and the results shown 
for A5 in Fig. 17 and for A7 in Fig. 18. These 
figures correspond exactly to the four figures 
just described for the steels of this paper. 

In all six figures, 13 to 18, it is to be noted 
that the lines have physical significance for the 
conditions of the present experiments only for 
strains at which fracture occurs, that is, for 
strains to the right of the point where the line 
for P crosses the zero axis. To the left of this 
point fracture is formally indicated when the 
specimen is exposed to a hydrostatic tension 
instead of a pressure. What the actual physical 
behavior would be under such a tension is for 
the present conjectural. 

A valid criterion for fracture will be repre- 
sented by a horizontal line in the figures. Of the 
various lines in the figures it is obvious at a 
glance that the lines for Z, and S come nearest 
to meeting this requirement. It also appears that 
S does not stand up as well as a criterion for the 
steels of this paper as it did previously for the 
two Arsenal steels A5 and A7. In this connection 
too much weight should not be given to the 
apparently large failure for the steel of hardness 
60, since the entire experimental range of strain 
was only 0.0 to 1.0 and the experimental error 
is therefore higher than for the others. The 
criterion that Z, is constant shows up as well on 
the average for the four steels of this paper as the 
criterion S. It is distinctly poorer for A5 and A7, 
however, so that on the whole it is not very 
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good. The physical significance of the criterion 
of constant Z, is that rupture occurs at a definite 
fiber stress, irrespective of the other stress com- 
ponents, and therefore irrespective of the com- 
ponent of elastic stretch in the same direction. 
This criterion has been used with some success 
by other investigators as a criterion of fracture 
in a limited range of conditions. However, in the 
wider range of conditions of this paper it would 
seem to have little physical significance, par- 
ticularly when it is considered that under pres- 
sure the fracture loses its tensile character and 
becomes entirely shearing. The criterion defi- 
nitely fails for such phenomena as the “‘pinching- 
off”’ effect. The physical significance of the 
criterion that S be constant is that fracture 
occurs when the average volume distension 
reaches a critical value. There would seem to be 
some physical plausibility in such a criterion, 
but the data of this paper, as well as considera- 
tions which I have already advanced in other 
papers, show that it cannot be the sole deter- 
mining factor in fracture. However, it seems to 
me that this must contain an important aspect 
of the mechanism responsible for fracture. 

Of the various criteria which have been sug- 
gested for fracture, it is obvious that no criterion 
in terms of strain only can possibly be valid 
under present conditions in view of the enormous 
strains that can be reached under high pressures 
without fracture. Of the simpler criteria in terms 
of stress, it is obvious that the maximum shearing 
stress criterion cannot be valid. The maximum 
shearing stress is simply }F, and inspection of 
Figs. 13-18 shows that F and therefore }F in all 
cases rises rapidly with strain, so that this 
criterion cannot be general. A physical argument 
could be made for the condition that F+HT 
should be constant at fracture, since this has 
the physical significance of maximum tensile 
fiber stress calculated from a state of uniform 
hydrostatic pressure as the zero. This criterion 
however, evidently fails by a wider margin than 
any other indicated by the figures. Or, since the 
total component of hydrostatic stress on the 
axis, where fracture begins, is 717—P, a similar 
argument could be made for the reasonableness 
of the criterion that F should be constant at 


8 P. W. Bridgman, Phil. Mag. 63-80 (July, 1912). 
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fracture. This criterion is also evidently wide of 
the mark. 

It is probably worth while calling attention 
to the fact that for all the steels the value of S 
at the lowest strain at which fracture occurs is 
approximately the same, namely 10,000 kg/cm’. 


SUMMARY 


The tensile properties of the two stainless 
steels with low nickel content fall into the same 
pattern as the carbon or low alloy steels pre- 
viously investigated. Two stainless steels with 9 
and 10 percent nickel content fall outside the 
previous pattern in that necking starts at an 
abnormally high elongation, they strain harden 
as an abnormally rapid rate, and the increase of 
ductility under pressure is abnormally low. The 
carbon steels with Rockwell C hardness in the 
range 48 to 63 all fall in the previous pattern, 
the rate of strain hardening being greater and 
the effect of pressure on ductility being less the 
greater the hardness. Even the steel of hardness 
63 shows a measurable ductility under a hydro- 
static pressure of 24,000 kg/cm’. 

A search is made. for a possible criterion of 
rupture by plotting various relevant stresses 
against the strain at fracture. No criterion 
emerges good under all conditions. The constancy 
of the mean hydrostatic tension (one-third the 
sum of the three principal stress components) 
remains the best criterion over the whole range 
of conditions, but in certain circumstances it 
may show considerable variation, and its margin 
of superiority over the criterion that the total 
fiber stress in the direction of fracture, Z., be 
constant at fracture, is not impressive. The 
maximum shearing stress criterion for fracture 
fails by a wide margin. No criterion of fracture 
in terms of strain components only can possibly 
be valid. Certain other plausible criteria on the 
stress components are also found not to hold. 

In addition to the generalizations obtained 
from previous work, namely that the strain 
hardening curves are linear in the strain and that 
the ductility is a linear function of the hydro- 
static pressure, two new empirical results emerge 
from the present method of plotting. These are: 
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first, that the component of hydrostatic tension 
on the axis of the specimen arising from necking 
is a linear function of the strain at fracture and, 
second, that the line of hydrostatic tension 
versus strain is approximately parallel to the line 
of flow stress versus strain. 

With regard to other metals than steel, 
aluminum, copper, bronze, and brass show an 
even greater effect of pressure on ductility than 
do the mild steels. These four metals were tested 
in the form of conventional tension specimens. 
The tensile properties of a number of other 
materials were tested by bending slender rods 
under pressure. Results were obtained in this 
way on hardened tool steels in various tempers 
in line with those obtained with the conventional 
tension specimens. Antimony, arsenic, tellurium, 
borax glass, cast iron, limestone (uncovered), and 
Carboloy showed no measurable plastic defor- 
mation before breaking by bending under hydro- 
static pressure. Bismuth shows large plastic 
effects. Hard rubber acquires great plasticity and 
also shows prominent time effects. Limestone, 
when covered, will take permanent plastic 
deformation in bending under hydrostatic pres- 
sure, but the interpretation of the results is 
obscured by a polymorphic transition which 
takes place under pressure. Although Carboloy 
showed no plasticity in bending, it showed a 
bending strength under 29,000 kg/cm? more than 
three times greater than that at atmospheric. 
Another type of experiment pointed to a similar 
very large increase of the elastic distortion which 
glass hard steel will support without fracture in 
a medium under hydrostatic pressure. In general 
comment on the bending experiments, it would 
seem that because of the complicated nature of 
the conditions the results are not susceptible of 
as straightforward an interpretation as are the 
conventional tension tests. This is pointed by the 
fact that carboloy will show plasticity in simple 
compression under pressure although it shows no 
plasticity in bending. The same thing is indi- 
cated by the marked difference of plasticity of 
limestone in simple tension and in bending. 

I am much indebted to Mr. L. H. Abbot for 
making the readings and manipulating the pres- 
sure apparatus. 
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Letter to the Editor 


GEorRGE G. BRUCK 
157 South Harrison Street, East Orange, New Jersey 
December 22, 1945 


N 1937 Woodyard published in the Proceedings of the 

I.R.E. a paper in which he described the use of an 
oscillator as frequency discriminator. Recently the growing 
interest in circuit arrangements permitting frequency dis- 
crimination without amplitude sensitivity has brought 
about such developments ‘as the ‘‘Ratio Detector.” 

The circuit described here is essentially an oscillator 
which falls in step with an incoming signal. The differential 
plate-current of the two halves of the tube is the dis- 
criminator output (Fig. 1). 

This circuit lends itself to many other applications like 
f.i. frequency modulation, frequency stabilization, fre- 
quency modulation phonograph pick-up, etc. 

A more complete paper on the theory and experimental 
results with this circuit will be published shortly. 
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Fic. 1. F-m discriminator. 





Here and There 


New Appointments 








Philip Nolan, who has spent the past three and a half 
years in fire control optics and infra-red research and de- 
velopment in the Navy’s Bureau of Ordnance, recently 
joined the staff of the Farrand Optical Company, Inc., of 
New York as its chief physicist. 


School and Society reports that G. W. Fox, professor of 
physics, lowa State College of Agriculture and Mechanic 
Arts (Ames), has been appointed technical adviser to 
General Douglas MacArthur. Dr. Fox will spend about 
three months in Japan making recommendations as ‘‘to 
which scientific activities should be encouraged . . . and 
which, if any, discouraged.” 


After more than three years of war research with the 
Office of Scientific Research and Development and the 
U.S. Navy Bureau of Ordnance, George M. L. Sommerman 
has become associate professor of electrical engineering at 
The Technological Institute of Northwestern University. 


Three new appointments were recently announced by the 
National Bureau of Standards. Hugh L. Dryden is now 
Assistant Director of the Bureau, Ferdinand G. Brickwedde 
is chief of the Heat and Power Division, and Robert Simha 
has become a member of the Division of Organic and 
Fibrous Materials with special responsibility for the coordi- 
nation of fundamental research on the molecular properties 
of high polymers. 
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Necrology 


Dr. Parke Benjamin Fraim, associate professor of physics 
at the Polytechnic Institute of Brooklyn, died January 13 
at the age of 59. 


Professor Harry Bateman, mathematician at California 
Institute of Technology, died in Utah on January 21 while 
on the way to New York, where he was to have been 
honored for outstanding work with Albert Einstein and 
other leading scientists. Professor Bateman was to be 
made a Fellow of the Institute of Aeronautical Sciences 
and receive a medal for his contributions to aeronautics. 
He had held full professorships in three fields—mathe- 
matics, physics, and acronautics—and was the author of 
many standard higher mathematics texts. At his death 
he was 63 years of age. 


Physicist Honored in Chicago 

Dr. Glenn T. Seaborg, co-discoverer of plutonium and 
elements 95 and 96, who is ‘‘on loan” to the Metallurgical 
Laboratory of the University of Chicago, has been named 
the outstanding young man of Chicago for 1945 by the 
Chicago Junior Association of Commerce. Dr. Seaborg has 
been on leave from the University of California for the 
past four years. He discovered plutonium in 1940. The 
announcement of his discovery, with three co-workers, of 
elements 95 and 96 was made last November. Dr. Seaborg, 
who is 33, has been in charge of the research at the Uni- 
versity of Chicago for separating plutonium in amounts and 
purity required for war purposes. He has also been pri- 
marily responsible for the chemical separation procedures 
which were used in connection with the manufacture of plu- 
tonium at Clinton, Tennessee, and Hanford, Washington. 
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Technical Translations 


The R. T. P.'s listed below are available from the 
Durand Reprinting Committee in care of California Insti- 
tute of Technology, Pasadena 4, California. 

These were selected for reprinting from the British 
Ministry of Aircraft Production R. T. P. List No. 122. 


2478 Flow Tests on Steel Tubes under Combined Tension and 

Hohenemker Torsion. (Z.A.M.M., Vol. Il, No. 1, February, 1931, 
pp. 15-19.) (5 pages) 

2480 Measurement of the Characteristics of Frequency- 


Stablein, W. Modulated Oscillations. (Elektrische Nachrichten 
Technik., Vol. 4, April, 1943, No. 20, pp. 102-111.) 
(12 pages) 

The Testing of Graduations of Theodolites and Uni 
versal Instruments. (Ztschr. f. Instrumentenkunde. 
Vol. 45, February, 1925, pp. 70-84.) (15 pages) 


2481 
Heuvelink, H. 


2492 On the Stability of Elastic Shells Pt. I. (Z.A.M.M.. 
Deuker, E. A. Vol. 23, April, 1943, No. 2, pp. 81-100.) (26 pages) 
2493 The Stability of Elastic Shells Pt. 11. (Z.A.M.M., Vol. 
Deuker, FE. A. 23, June, 1943, No. 3, pp. 169-179.) (13 pages) 

2494 Oscillations in Gas Columns. (Forschung, Vol. 6, No. 6, 
Pischinger, A. November—December, 1935, pp. 273-280.) (12 pages) 
2495 The Induction and Discharge Process and Limitations of 


Schmidt, F. R. Poppet Valve Operation in a Four-Stroke Aero Engine. 
(Luftfahrtforschung, Vol. 16, No. 5, 20.5.1939, pp. 
251-275.) (45 pages) 

Rivet Load Distribution in Lap Joints of Constant Cross- 
Section under Tension. (Heinkel Werke, Rostock L. F. 
F., Vol. 15, No. 4, 20.1.38, pp. 41-47.) (14 pages) 
Measurement of Direction Velocity and Pressure in a 
Three Dimensional Current. (Luftfahrtforschung, Vol. 
19, No. 9, 20.10.42, pp. 326-330.) (10 pages) 
Physical-Chemical Effects of Supersonic Waves. (Archiv 
f. d. Eisenhuttenwesen, Vol. 12, No. 4, October, 1938, 
pp. 185-193.) (17 pages) 


2497 
Volkersen, O. 


2498 
Jegorow, G. 


2499 
Hiedemann, E. 


2500 Graphical Determination of Wall Temperature during 

Lutz, O. Heat Transfer. (Forschung, Vol. 6, No. 5, 1935, pp. 
240-244.) (7 pages) 

2501 The Theory of One-Sided Web Stiffeners. (Der Stahlbau. 

Chwalla, E. Suppl. to Bautechnik, Vol. 10, No. 10, May 7, 1937, 

Novak, A. pp. 73-76.) (10 pages) 


Bulletin of the Atomic Scientists of Chicago 


The Atomic Scientists of Chicago have published several 
issues of the first volume of their Bulletin, an eight-page 








printed leaflet. Thus far articles have appeared under the 


following headings: 


Technical Feasibility of Atomic Energy Controls 

General Groves on the Future of the Atomic Bomb Projects 
The Federation of Atomic Scientists 

Denver Conference on Atomic Energy 

UNO Faces the Problem of Atomic Energy Control 
Structure of the UNO 

British Public Opinion and the Atomic Bomb Menace 

The Atomic Scientists Back the McMahon Bill 

The Atomic Scientists and the Federation of American Scientists 
The Distribution of Uranium in Nature 

UNO Creates Atomic Energy Commission 


In addition, the-Bulletin always carries several columns 
under the heading ‘‘Congressional News.” 

The Bulletin has been featuting a series of articles on 
“Control Problems.”’ Thus far it has included the following: 
general survey of the control problem, control of uranium 
mining, and control of a gaseous diffusion plant. The 
Bulletin is published at 1126 East 59th Street, Chicago 37, 
Illinois. 


New Firm of Consulting Physicists 


Paul Rosenberg Associates, consulting physicists, have 
announced the opening of offices in New York City. Work 
will be done in radar, ultrasonics, television, optics, nu- 
clear physics, and general physics. The head of the firm, 
Paul Rosenberg, has recently completed four years’ work 
at the Radiation Laboratory of the National Defense 
Research Committee at the Massachusetts Institute of 
Technology. 


Dr. Meitner Arrives 


Dr. Lisa Meitner, who is credited with establishing the 
mathematical data and equations in the breakdown of the 
uranium atom, arrived in New York from Stockholm on 
January 25. She came to the United States to join the 
faculty of the Catholic University of America in Wash- 
ington. 





New Booklets — 











_ | Rheology Bulletin for December 1945, publication of The 
Society of Rheology, has the following table of contents: 


Bibliography of Books on Viscosity, Viscometry and Rheology 

Society of Rheology Annual Meeting Abstracts 

Rheology Progress Abstracts (covering apparatus and methods, cellu- 
lose, clays, fibers, foams, general, ink, paints, petroleum products, 
plastics, polymers, proteins, theory) 






General Radio Company, 275 Massachusetts Avenue, 
Boston 39, Massachusetts, has issued a 68-page Postwar 
Supplement to its Catalog K, now in the fourth edition. 
The Supplement brings all editions completely up to date. 
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Westinghouse Newsfront for January 1946, published by 
Westinghouse Electric Corporation, Pittsburgh 30, Penn- 
sylvania, features articles on devices used in radar, re- 
search into the causes of rust and corrosion, industrial use 
of high frequency “‘cooking,”’ and the value of electricity 
in increasing production, efficiency, and comfort on Ameri- 
can farms. 


Jarrell-Ash Company of 165 Newbury Street, Boston 16, 
Massachusetts, has issued a 16-page catalog of spectro- 
graphic equipment for research in spectroscopy and for 
spectrochemical analysis in metallurgy, mineralogy, agron- 
omy, chemistry, medicine, and other fields. 
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